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Abstract 

Lasers  offer  tremendous  advantages  over  RF  communication  systems  in  terms 
of  bandwidth  and  security  due  to  their  ultra-high  frequency  and  narrow  spatial 
beamwidth.  Unfortunately,  atmospheric  turbulence  significantly  increases  the  re¬ 
ceived  power  variation  and  bit  error  rate  (BER)  in  free-space  optical  communication 
(FSOC)  systems.  Further,  airborne  optical  communication  systems  require  special 
considerations  in  size,  complexity,  power,  and  weight. 

If  two  or  more  laser  beams  are  sufficiently  separated  so  that  their  turbulence 
effects  are  uncorrelated  (i.e.  anisoplanatic),  they  can  effectively  “average  out”  tur¬ 
bulence  effects.  This  requisite  separation  distance  is  derived  for  multiple  geometries, 
turbulence  conditions,  and  optical  properties.  In  most  cases  and  geometries,  the 
angles  ordered  from  largest  to  smallest  are:  phase  uncorrelated  angle  (equivalent 
to  the  tilt  uncorrelated  angle  and  phase  anisoplanatic  angle),  tilt  isoplanatic  angle, 
phase  isoplanatic  angle,  scintillation  uncorrelated  angle  (or  scintillation  anisoplanatic 
angle),  and  scintillation  isoplanatic  angle  (9^ind  >  Ota  >  0o  >  0Xind  >  0Xo).  Con¬ 
ventional  adaptive  optics  (AO)  systems  only  correct  for  phase  and  cannot  correct  for 
strong  scintillation,  while  multiple-transmitter  systems  use  several  transmission  paths 
to  “average  out”  effects  of  the  strong  scintillation  by  incoherently  summing  up  the 
beams  in  the  receiver. 

Since  all  three  airborne  geometries  (air-to-air,  air-to-ground,  and  ground-to- 
air)  are  studied,  a  comparison  of  multiple-beam  airborne  laser  communication  system 
performance  is  presented  for  the  first  time.  Wave  optics  simulations  show  that  a 
combination  of  transmitter  diversity,  receiver  and  transmitter  trackers,  and  adaptive 
thresholding  can  significantly  reduce  BER  in  an  air-to-air  FSOC  system  by  over  10,000 
times.  As  demonstrated  in  this  work,  two  transmitters  alone  separated  by  only  31  cm 
(100  km  path  length,  1.55  /mi  wavelength,  4  km  in  altitude)  provide  a  significant 
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BER  improvement  over  one  transmitter,  especially  for  the  strong  turbulence  regime 
where  the  required  SNR  for  a  fixed  BER  is  reduced  by  9  dB.  Including  the  track¬ 
ing  and  adaptive  thresholding  techniques,  resulted  in  a  13  dB  overall  improvement. 
Two  beams  also  reduce  the  fade  length,  suggesting  even  greater  improvement  can  be 
obtained  when  interleaving  and  forward  error  correction  coding  is  implemented. 
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I.  Introduction 

1.1  Motivation  and  challenges 

Laser  communications  offer  tremendous  advantages  over  radio  frequency  (RF)  in 
bandwidth  and  security  due  to  the  ultra-high  frequencies  and  point-to-point  na¬ 
ture  of  laser  propagation.  In  addition,  optical  transmitters  and  receivers  are  much 
smaller  and  lighter  than  their  RF  counterparts  and  operate  at  much  lower  power  lev¬ 
els.  Current  airborne  sensors  are  collecting  data  at  an  ever-increasing  rate.  With 
the  advent  of  hyperspectral  imaging  systems  and  other  sensors,  this  trend  will  con¬ 
tinue  as  two-dimensional  data  are  replaced  by  three-dimensional  data  cubes  at  finer 
resolutions.  Current  RF  communication  systems  cannot  keep  up  with  this  trend.  In 
addition,  free-space  optical  communication  (FSOC)  systems  could  provide  covert,  dif¬ 
ficult  to  jam  or  intercept,  high-speed,  broadband  connectivity  to  airborne  (especially 
low-observable)  platforms  [86]. 

However,  airborne  laser  communications  are  severely  affected  by  clouds,  dust, 
and  atmospheric  turbulence,  causing  deep,  long  fades  at  the  receiver.  Ultimately  a 
hybrid  communication  system  which  includes  RF  communication  is  necessary,  since 
clouds,  fog,  or  dust  occasionally  obstructs  the  path  for  laser  communication.  Even 
when  the  channel  is  clear,  the  same  atmospheric  turbulence  effects  that  limit  the  res¬ 
olution  of  optical  systems  and  make  the  stars  twinkle  can  severely  reduce  the  received 
power,  causing  long  deep  fades.  This  atmospheric  turbulence  in  the  propagation  path 
causes  the  laser  beam  to  wander,  spread  and  break  up.  These  effects  can  cause  the 
received  signal  power  to  drop  below  the  receiver’s  threshold  for  milliseconds  at  a  time. 
For  a  10  Gbit/s  binary  FSOC  system,  a  millisecond  fade  means  at  least  5  million  bit 
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errors.  Since  the  turbulence  of  an  air-to-air  link  extends  along  the  entire  path  and 
causes  long,  deep  fades,  simply  turning  up  the  power  would  not  be  effective.  In  addi¬ 
tion,  airborne  FSOC  systems  require  special  considerations  in  size,  complexity,  power, 
and  weight. 

There  are  essentially  two  different  ways  to  improve  this  condition:  increase  the 
signal  diversity  to  average  out  the  effects  or  compensate  for  the  turbulence  condi¬ 
tions  in  real  time.  The  first  approach  is  implemented  in  this  research  by  designing  a 
multiple-transmitter  system  where  the  uncorrelated  effects  of  each  path  are  averaged 
at  the  receiver. 

1.2  Spatial  diversity  and  redundancy 

In  the  night  sky,  objects  like  the  moon  and  the  planets  do  not  twinkle  like  the 
stars.  The  uncorrelated  amplitude  fluctuations  of  the  optical  wavefronts  from  these 
extended  incoherent  objects  are  received  by  our  eyes  and  these  scintillation  effects  are 
averaged  out.  Multiple  laser  transmitters  incoherent  with  each  other  and  adequately 
separated  act  in  much  the  same  way. 

Atmospheric  turbulence  causes  random  fluctuations  of  a  wavefront’s  phase  and 
amplitude  in  space  and  time  due  to  the  winds  and  platform  velocities.  Two  beams 
traveling  along  the  same  path  experience  correlated  effects,  but  as  the  separation 
distance  increases  these  effects  decorrelate.  If  these  beams  are  separated  so  that  the 
effects  are  relatively  uncorrelated,  each  path  experiences  fades  at  different  moments  in 
time,  allowing  a  multi-path  system  to  “average  out”  the  fades.  As  with  a  system  where 
redundancy  is  built  in  to  bypass  a  failure  of  a  particular  system,  a  multiple  transmitter 
system  provides  multiple  paths  for  the  information  to  the  reach  the  receiver.  This 
research  analytically  determines  the  requisite  separation  distances  for  uncorrelated 
phase  and  amplitude  effects  and  verifies  the  results  and  improvements  in  simulation 
for  an  air-to-air  FSOC  system. 


2 


1.3  Adaptive  system  for  temporal  signal  fluctuations 

Once  the  system  receives  the  signal,  a  binary  digital  receiver  must  determine 
if  a  ‘1’  or  ‘0’  was  sent  by  measuring  the  signal  and  determining  whether  that  signal 
level  is  more  indicative  of  a  ‘1’  or  a  ‘O’.  This  is  accomplished  by  using  an  appropriate 
threshold  to  decide  if  a  ‘1’  or  a  ‘O’  was  sent.  Specifically,  if  the  signal  exceeds  the 
threshold,  it  registers  a  ‘1’  and  if  it  does  not,  it  registers  a  ‘O’.  For  a  digital  com¬ 
munication  system,  a  bit  error  occurs  when  a  ‘1’  is  transmitted  and  a  ‘0’  is  detected 
or  when  a  ‘0’  is  transmitted  and  a  ‘1’  is  detected.  Optimal  threshold  determina¬ 
tion  depends  upon  the  conditional  probabilities  of  the  signal  level  at  the  receiver  for 
each  symbol  transmitted  [82],  Since  atmospheric  random  fluctuations  cause  these 
conditional  probabilities  to  vary  over  time,  it  is  beneficial  to  adaptively  change  the 
threshold  given  the  current  signal  level  to  provide  the  maximum-likelihood  optimal 
solution  for  each  time  slice  [14, 18,47].  The  rate  at  which  the  turbulence  causes  fades 
and  bit  errors  at  the  receiver  is  much  slower  than  the  data  rate,  thus  a  system  can 
be  designed  to  adaptively  change  this  threshold  with  time.  This  research  designs  and 
tests  by  simulations  realistic  adaptive  threshold  systems  and  measures  the  bit  error 
rate  (BER)  performance  improvement. 

1-4  Scope  and  assumptions 

The  links  tested  in  this  research  are  primarily:  air-to-air  links  at  distances 
of  no  more  than  100  km  at  10  km  in  altitude  and  air-to-ground  and  ground-to-air 
links  (surface  to  10  km)  at  elevation  angles  between  20  and  90  degrees.  The  sim¬ 
ulations  use  the  Hufnagel- Valley  turbulence  profile  with  W  —  21  m/s  (denoted  by 
HV-21)  [66].  These  links  are  shown  graphically  in  Fig.  1  and  the  HV-21  turbulence 
profile  is  shown  in  Fig.  4  in  Chapter  2.  Most  regimes  tested  in  this  research  are  in 
the  weak-turbulence  regime  where  there  is  little-to-no  saturation  of  the  intensity  vari¬ 
ations.  In  the  moderate-to-strong  turbulence  regime  the  amplitude  effects  (intensity 
variations,  or  scintillation)  of  the  turbulence  saturate,  whereas  the  phase  effects  do 
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Figure  1:  Airborne  laser  communication  link  scenarios  tested  in  this  research.  Air-to- 
air  links  of  100  km  at  10  km  in  elevation.  Surface  to  10  km:  air-to-ground 
and  ground-to-air  links  at  elevation  angles  of:  20-90  degrees  [20]. 
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not.  Moderate-to-strong  turbulence  regimes  are  identified  and  handled  on  a  case-by- 
case  basis. 

In  most  cases,  the  atmospheric  turbulence  is  assumed  to  be  isotropic  (statistical 
properties  independent  of  direction)  and  homogeneous  (statistical  properties  indepen¬ 
dent  of  position).  The  strength  of  turbulence  varies  with  altitude.  At  each  altitude 
the  turbulence  is  assumed  to  be  isotropic  and  homogeneous  [2,66,69]. 

The  turbulent  air  around  an  aircraft  is  neither  isotropic  nor  homogeneous  and 
cannot  be  modeled  in  the  same  way  as  atmospheric  turbulence.  Detailed  analysis  of 
aero-optic  effects  have  not  been  covered  in  this  research.  Losses  due  to  atmospheric 
absorption  and  extinction  are  estimated  and  attributed  as  a  loss  when  calculating  the 
received  power. 

To  compare  transmitter  systems,  one  receiver  size  is  used  for  all  investigations. 
In  practice,  the  collection  optic  should  be  made  as  large  as  possible.  Not  only  do  larger 
telescopes  collect  more  power,  the  angle-of-arrival  and  intensity  variations  decrease  as 
well.  However,  the  larger  the  aperture  the  worse  the  effects  due  to  higher-order  phase 
perturbations.  Aperture  averaging  of  the  intensity  occurs  when  the  receiver  aperture 
is  larger  than  the  correlation  width  pcw  (i.e.  characteristic  size  of  highly  correlated 
intensity  at  the  receiver)  [2]. 

1 . 5  Summary  of  main  results 

The  goal  of  this  research  is  to  develop  techniques  and  systems  to  improve  two- 
way  free-space  laser  communication  between  a  ground  station  and  an  airborne  plat¬ 
form  and  among  airborne  platforms.  The  techniques  developed  in  this  research  are 
designed  to  minimize  the  complexity,  size,  weight,  and  power  requirements,  especially 
for  airborne  platforms. 

Significant  reductions  in  bit  error  rate  are  attained  by  implementing  optimized 
multiple-transmitter  systems  to  average  out  the  deleterious  effects  of  turbulence. 
Through  analysis  and  simulation,  optimal  separations  for  double-transmitter  sys- 
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terns  for  various  geometries  and  tracking  systems  were  determined.  This  research 
derived  for  the  first  time  the  optimal  angular  separations  for  multiple  transmitter 
systems  for  airborne  FSOC  systems.  A  majority  of  previous  research  on  multi¬ 
ple  transmitters  focused  on  satellite  communications  (where  the  turbulence  is  only 
present  over  a  small  part  of  the  propagation  path)  or  constant  turbulence  strength 
paths  [4,5,29,59,60,68].  The  angular  separations  for  all  three  airborne  geometries 
(air-to-air,  air-to-ground,  and  ground-to-air)  are  determined. 

Previous  research  on  isoplanatism  defined  the  maximum  angle  over  which  the 
turbulence  effects  between  two  paths  is  relatively  similar  [26,69].  These  isoplanatic 
angles  have  been  determined  for  the  tilt  variance,  higher-order  phase  variance,  and 
the  scintillation  (intensity  variance).  Here,  that  work  is  extended  to  determine  the 
minimum  angle  at  which  the  paths  are  relatively  different.  The  less  the  effects  between 
the  paths  are  correlated,  the  better  the  averaging  effect  of  multiple  beams.  Since  all 
three  airborne  geometries  are  studied,  a  comparison  of  multiple-beam  airborne  laser 
communication  system  performance  is  presented  for  the  first  time. 

Using  these  optimal  separation  distances,  the  temporal  impacts  of  multiple 
transmitters  are  studied  along  with  an  adaptive  threshold  system.  Wave-optics  sim¬ 
ulations  show  multiple  transmitters,  receiver  and  transmitter  trackers,  and  adaptive 
thresholding  significantly  reduces  BER  (by  a  factor  of  10,000).  Two  transmitters  alone 
provide  a  significant  BER  improvement  over  one  transmitter,  especially  for  the  strong 
turbulence  regime  with  up  to  a  9-decibel  (dB)  improvement  gain.  Adaptive  thresh¬ 
olding  systems  provide  significant  improvement  over  optimal  fixed  thresholds  for  both 
single-  and  double-transmitter  systems,  providing  an  additional  3-4  dB  over  both  sys¬ 
tems.  This  indicates  that  the  improvement  provided  by  these  multiple  techniques  can 
be  combined  to  provide  even  further  improvement.  There  are  further  trickle-down  ef¬ 
fects,  since  two  beams  also  reduce  the  fade  length,  indicating  they  would  likely  provide 
even  greater  improvement  with  interleaving  and  forward-error-correction  coding. 
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1.6  Dissertation  organization 

Background  on  atmospheric  turbulence  effects  and  the  impacts  on  digital  laser 
communication  are  covered  in  Chapter  II.  Next,  Chapter  III  examines  what  others 
have  done  to  address  the  challenges  of  laser  communication.  Chapter  IV  derives  the 
requisite  separation  distances  and  angles  for  averaging  the  effects  by  determining  the 
minimum  separation  for  uncorrelated  phase  and  amplitude  effects.  In  Chapter  V,  this 
separation  is  implemented  in  a  time  series  simulation  to  determine  the  fade  statistics 
as  well  as  the  tracking  system  performance.  In  addition,  an  ideal  and  realistic  adap¬ 
tive  threshold  system  is  implemented  and  compared  with  an  optimal  fixed  threshold 
system.  Finally,  in  Chapter  VI  conclusions  for  the  research  effort  are  summarized  and 
recommendations  for  future  research  are  discussed. 
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II.  Background 


This  chapter  briefly  describes  laser  propagation  through  atmospheric  turbulence 
and  optical  digital  communication.  The  reader  is  referred  to  [2,  66,  69]  for  a 
more  complete  treatment  of  atmospheric  turbulence  effects  and  [1,56]  for  additional 
considerations  of  optical  communications.  In  addition,  subsequent  chapters  of  this 
dissertation  describe  in  more  detail  special  considerations  and  assumptions  used  to 
study  these  effects. 

2.1  Propagation  of  Gaussian-beam  waves 

This  first  section  describes  how  a  laser  beam  propagates  in  a  vacuum.  Laser 
beams  can  propagate  in  different  modes.  In  most  cases,  the  lowest-order  transverse 
electro-magnetic  (TEM)  Gaussian-beam  mode  TEMoo  is  a  good  mode  to  work  with, 
because  the  limiting  cases  include  an  infinite  plane  wave  and  a  spherical  wave.  This 
research  concentrates  its  analysis  on  this  lowest-order  TEM0o  mode. 

A  Gaussian  laser  beam  is  parameterized  by  its  1/e  radius  W0,  its  radius  of 
curvature  F0,  and  the  optical  wavelength  A.  The  initial  optical  field  £/o(r,  0)  of  the 
TEMqo  wave  exiting  the  laser  at  z  —  0,  centered  on  the  optical  axis  (r  =  0)  is  [2] 


(1) 


1  /O 

where  r  =  ( x 2  +  y2)  is  the  radial  distance  from  the  optical  axis  and  a o  is  optical  field 
amplitude  in  (W/m2)1/2.  The  respective  amplitude  and  phase  of  the  Gaussian-beam 


wave  are 
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Solving  the  wave  equation  for  the  Gaussian  beam  propagation  results  in  [2] 
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Figure  2:  Diffraction  effects  on  collimated  laser  beams  of  differ¬ 
ent  beam  widths. 


where  ©o  and  A0  are  the  set  of  input  plane  beam  parameters  defined  by  [2] 
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The  new  spot  size  radius  W(z)  and  the  new  radius  of  curvature  F(z)  at  z  and  can  be 
described  by  the  input-plane  beam  parameters  as  [2] 


W  =  Wo(A20  +  Ql)1'2 


(6) 


and 

f„  (eg  +  Aq)  (go  - 1) 

eg  +  Ag-e,,  '  [n 

Figure  2  shows  the  spot  size  radius  of  three  different  size  collimated  beams  (i.e. 
F{)  — >  oo)  for  a  propagation  distance  out  to  L  —  100  km.  The  plot  shows  the 

very  strong  inverse  relationship  of  input  beam  size  on  the  beamwidth  for  long  laser 
propagations. 
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2.2  Atmospheric  turbulence  effects 

2.2.1  Index  of  refraction  fluctuations  in  the  atmosphere.  The  atmospheric 
index  of  refraction  fluctuations  occur  due  to  temperature  and  velocity  differences 
caused  by  turbulent  air  movement.  Optical  waves  that  propagate  through  the  at¬ 
mosphere  are  distorted  by  these  temporal  and  spatial  refractive  index  fluctuations. 
These  optical  wave  distortions  make  the  stars  twinkle,  limit  the  resolution  of  imaging 
systems,  and  degrade  the  propagation  of  lasers. 

This  atmospheric  turbulence  caused  by  the  heating  and  cooling  of  the  earth 
and  its  atmosphere  lead  to  large-scale  variations  in  the  air  temperature.  Kolmogorov’s 
theory  explains  how  these  temperature  differences  cause  winds  that  break  up  the  large 
scale  inhomogeneities  into  smaller  ones  as  the  laminar  flow  transitions  to  turbulent 
flow.  This  turbulent  flow  causes  random  pockets  of  temperature  differences  that  vary 
randomly  in  time  and  space. 

The  index  of  refraction  n  is  sensitive  to  temperature  differences.  At  optical 
frequencies,  the  index  of  refraction  as  a  function  of  air  pressure  (which  drives  the  air 
density)  and  temperature  is  determined  by  [36] 

77. 6  P 

n  —  tiq  4 - — —  x  1CT6,  A  =  500  nm  (8) 

where  no  —  1,  T  is  the  temperature  in  Kelvins,  and  P  is  the  air  pressure  in  millibars. 
The  temperature  fluctuations  normally  dominate  the  index  variations,  although  some 
turbulence  around  airframes  might  provide  significant  air  density  changes  to  appre¬ 
ciably  affect  the  index  of  refraction. 

Over  the  last  50  years  since  the  invention  of  the  laser  in  1958  [71]  and  the 
first  infrared  sensors,  modeling  optical  propagation  through  atmospheric  turbulence 
has  become  increasingly  important.  Since  the  turbulence  is  random  and  cannot  be 
predicted,  several  models  have  been  developed  to  statistically  characterize  how  the 
turbulence  fluctuates.  A.N.  Kolmogorov’s  theory  is  the  most  widely  accepted,  due  to 
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its  simplicity  and  consistent  agreement  with  observations  [41].  Much  of  the  research 
since  is  based  on  his  theory. 

Air  masses  of  uniform  index  of  refraction  are  commonly  referred  to  as  turbu¬ 
lent  eddies.  Kolmogorov  assumed  that  these  eddies  are  statistically  homogeneous  and 
isotropic.  Assuming  homogeneity  and  isotropy  allow  relations  to  be  developed  in  the 
spatial  and  spectral  domains  to  characterize  the  statistical  nature  of  the  fluctuations 
as  random  processes.  Three  relations  in  particular  are  used  to  characterize  the  ran¬ 
dom  processes  that  describe  turbulence:  the  three-dimensional  covariance,  structure 
function,  and  power  spectral  density  (PSD). 

First,  consider  a  time- varying  complex  random  process  /( r),  where  r  is  a  three- 
dimensional  spatial  vector.  In  many  cases,  the  three-dimensional  (3-D)  auto-covariance 
B(r i,r2)  defined  by 

Bf(r i,r2)  =  ([/(ri)  -/i(ri)][r(r2)  -/i*(r2)]),  (9) 

is  sufficient  to  describe  a  stationary  random  process,  where  //(r *)  is  the  expected  value 
of  /(rj)  and  (•)  refers  to  the  ensemble  average.  A  wide-sense  stationary  (WSS)  random 
process  is  a  system  where  the  variance  and  mean  remain  constant  and  do  not  vary 
with  time.  Any  strictly  stationary  process  which  has  a  mean  and  a  covariance  is  also 
WSS.  If  a  WSS  random  process’s  temporal  variance  and  mean  can  be  approximated 
by  the  ensemble’s  variance  and  mean,  it  is  said  to  be  ergodic  in  the  variance  and 
mean. 

The  relation  that  links  the  spectral  representation  to  the  spatial  one  is  the 
Wiener-Khinchin  theorem.  It  basically  states  the  power  spectral  density  of  a  WSS 
random  process  is  equal  to  the  Fourier  transform  of  the  autocovariance  function.  The 
PSD  <E>^(k)  of  /  is  given  by 

f  Bf(r)  exP  i-JK  '  r)  dr,  (10) 
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where  n  is  the  3-D  angular  spatial  frequency  vector.  If  a  random  process  /( r)  has  a 
slowly  varying  mean,  it  cannot  be  represented  by  the  auto-covariance.  If,  for  small 
separations  r,  the  process  is  homogeneous  in  the  mean,  it  is  said  to  have  stationary 
increments.  Accordingly,  the  3-D  structure  function  Df( r)  adequately  describes  this 
process  and  is  denoted  by 

dAt)  =  (I/Od)  -/(ri  +  r)!J>-  (n) 

Furthermore, if  the  process  is  homogenous  expanding  the  structure  function  yields  the 
following  convenient  relationship  to  the  auto-covariance: 

£>/(r)  =  +  (/2(ri  +  r))  -2(/(ri)/(ri  +  r)) 

=  2  [5,(0)  -  B,{ r)] .  (12) 


If  the  medium  is  isotropic,  meaning  the  relationships  are  the  same  regardless  of 
the  direction,  the  vector  equations  can  be  converted  to  scalar  separation  distances. 
For  the  isotropic  case,  the  structure  function  is  determined  by 


Df(r)  =  87T 


/»oo 

sin  (nr) 
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The  inverse  relationship  is 


(13) 


<*/(k)  = 
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Jo  ^ 

r  —Df(r) 
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nr 


dr. 


(14) 


The  refractive  index  3-D  turbulence  power  spectrum  (n,  z )  is  a  function  of 

both  the  turbulence  along  the  path  and  the  spatial  frequency.  The  variable  z  repre¬ 
sents  the  location  along  the  propagation  path,  and  n  denotes  the  spatial  frequencies 
of  the  turbulence.  The  smaller  frequencies  are  bounded  by  no  =  2n /Lq  where  L0  is 
the  outer  scale,  representing  the  largest  sizes  of  turbulence.  The  larger  frequencies 
are  bounded  by  nm  =  2i r//0,  where  l0  is  the  inner  scale.  As  the  name  suggests,  the 
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inner  scale  refers  to  the  smallest  portions  of  the  turbulence.  If  the  turbulence’s  spatial 
frequency  is  in  the  range  27 r/L0  <  hi  <  2tt/10,  the  turbulence  power  spectrum  can  be 
represented  by  the  Kolmogorov  turbulence  power  spectrum  [36] 


4>*(K,z)=  0.033C„2(z)K-11''3. 


(15) 


where  Cn(z)  is  the  refractive  index  structure  parameter.  The  modified  von  Karman 
spectrum  (hi,  z )  takes  the  inner  and  outer  scales  into  account  and  is  represented 


by  [36] 


(16) 


where  nm  =  27 r/l0  and  k0  —  27 r/L0.  The  original  von  Karman  spectrum  does  not 
include  the  exponential  factor  for  the  inner  scale.  Both  spectrums  are  used  in  this 
research  depending  upon  which  is  appropriate.  If  the  turbulence  is  limited  to  the  range 
27 r/L0  <  n  <  2i r/Z0,  the  Kolmogorov  turbulence  power  spectrum  is  used.  Otherwise, 
the  von  Karman  spectrum  is  used. 

The  inner  scale  plays  a  large  part  into  determining  the  scintillation,  since  these 
small  scale  variations  break  up  the  wavefront  and  cause  amplitude  variations.  The 
inner  scale  length  varies  with  altitude.  Holding  the  energy  dissipation  rate  constant, 
the  inner  scale  varies  nearly  inversely  with  the  air  density  [85] 


(o  cc  p  3/i. 


(IT) 


At  ground  level,  inner  scale  measurements  typically  range  from  2  mm  to  9  mm  [85].  A 
comprehensive  study  accomplished  over  White  Sands,  NM  determined  the  inner  scale 
length  increased  exponentially  from  l0  =  1  cm  at  5  km  to  /0  =  8  cm  at  19  km  [21,85]. 
Applying  the  model  in  Eq.  (17)  and  Z0  =  1  cm  at  5  km  a  model  for  the  inner  scale  vari¬ 
ation  with  altitude  was  determined.  Since  the  air  pressure  at  5  km  is  approximately 
54  kPa,  this  equates  to  an  air  density  p  at  5  km  of  0.843  kg/m3.  The  equation  for 
inner  scale  used  in  this  research  is  determined  by  solving  for  the  constant  multiplier 
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in  Eq.  (17),  leading  to 


l0 

0.01 


C\ 

p3/4 

Cto 

0.8433/4 


where  the  air  density  is  (assuming  an  ideal  gas) 


(18) 


P  = 


PM 

~KT ’ 


(19) 


and  P  is  the  air  pressure  in  pascals,  M  =  0.0289644  kg/mol  is  molecular  weight  of  dry 
air,  R  =  8.31447  is  the  universal  gas  constant  in  J/(mol  K),  and  T  is  the  temperature 
in  Kelvins.  Solving  for  the  constant  yields  Ci0  =  0.0088  m~5/4  kg3//4,  which  is  used  to 
determine  the  inner  scale  for  the  simulations  performed  in  this  research,  as  shown  in 
Fig.  3. 


Large-scale  turbulence  affects  the  phase,  and  following  this  trend,  the  largest 
scale  turbulence  affects  the  wavefront  tilt  (i.e.  largest  scale  phase  effects)  and  causes 
laser  beams  to  randomly  wander.  This  outer  scale  denoted  by  L0  varies  with  altitude 
and  tends  to  be  an  ellipsoidal  function  of  zenith  angle.  Near  the  ground  Lq  ps  OAh, 
where  h  denotes  the  altitude,  as  shown  in  Fig.  3.  At  higher  altitudes,  the  outer  scale 
is  determined  by  the  vertical  outer  scale  and  the  horizontal  outer  scale.  The  vertical 
outer  scale  L0vert  typically  varies  from  10  m  to  70  m  [21].  The  horizontal  outer  scale 
is  very  large.  Aircraft  measurements  have  determined  it  can  be  over  hundreds  of 
kilometers  [85].  For  horizontal  propagation,  an  infinite  outer  scale  is  used  in  this 
research.  However,  for  slant  ranges,  the  outer  scale  is  determined  by  taking  a  slice 
through  the  vertical  outer  scale 

L°  ~  cos(0' 
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Figure  3:  This  plot  shows  the  inner  scale  as  it  varies  with 
altitude  calculated  from  Eq.  (17),  calibrated  by 
Iq  =  1  cm  at  5  km  altitude  from  measurements  taken 
at  White  Sands  Missile  Range,  NM.  The  outer  scale  Lq 
was  determined  from  the  approximation  of  Lq  0.4h 
and  the  measurements  at  White  Sands  [21]. 
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2.2.2  Turbulence  models  and  moments.  The  refractive  index  structure  pa¬ 
rameter  characterizes  the  strength  of  atmospheric  turbulence.  There  are  many  differ¬ 
ent  models  used  to  describe  the  effects  of  atmospheric  turbulence  as  they  vary  with 
altitude.  A  very  common  model  that  is  used  throughout  this  work  to  describe  the 
turbulence  strength  C2(h)  is  the  Hufnagel- Valley  model.  This  turbulence  profile  is 
described  by  [66] 

Cm  =  5.94  x  10- (^)  V  exp  (h^) 

+  2.7X10- exp  (— ^)mexp(—  A),  (21) 


where  h  is  the  height  above  the  ground  in  meters,  A  in  m -2//3  sets  the  turbulence 
strength  near  the  ground,  and  v  represents  the  high-altitude  root-mean  square  (rms) 
wind  speed  in  m/s.  To  change  the  strength  of  turbulence  at  high  altitudes,  the  v  term 
is  typically  varied.  This  research  uses  the  common  value  for  A  of  1.7  x  10  14m  2/3  and 
sets  the  v  variable  to  21  m/s  which  is  often  referred  to  as  either  HV-21  or  HV5/7.  The 
5  and  7  in  HV5/7  refer  to  a  coherence  diameter  of  ro  =  5  cm  and  an  isoplanatic  angle 
of  0o  =  7  /urad  at  a  wavelength  of  0.5  /irn  for  zenith  imaging  through  the  atmosphere. 

The  rms  wind  speed  v  in  Eq.  (21)  is  calculated  from  [2] 


v  = 


r>20x  103 


15  x  lo3  y5xl0  3 


V2{h)  d h 


1/2 


(22) 


and  substituted  back  into  Eq.  (21).  The  propagation  path  moments  fim  frequently 
appear  in  equations  that  describe  specific  atmospheric  turbulence  effects.  The  full 
moments  for  ground-to-space  propagation  are  [31,69] 


roo 

C2(z)zm  dz  =  secm+1  (£)  /  C2(h)hm  dh, 

Jo 


(23) 
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Figure  4:  HV-21  and  the  SLC-Day  strength  of  turbulence  profile 
models. 

where  £  is  the  zenith  angle.  For  horizontal  propagation  or  constant  turbulence 
strength,  the  moments  simplify  to 

=  C,y<”+1>/(m  +  1),  (24) 

where  L  is  the  propagation  path  length.  For  a  slant  propagation  path  length  L 
through  the  atmosphere  to  a  height  77,  where  L  =  sec  (£)  77,  the  turbulence  moments 
can  be  broken  into  these  partial  moments  [69] 

=  /  C’(z)z”  dz  =  sec”+1  (f)  [  Cl(h)hTAh,  (25) 

Jo  Jo 

poo  poo 

/  C2n(z)zm dj  =  Secm+1 «)  /  Cl(h)hTAh.  (26) 

J  L  J  H 

The  equations  for  and  /i+  above  apply  to  a  flat-earth  model.  The  curvature  of 
the  earth  must  be  included  when  propagating  long  distances  at  low  elevation  an¬ 
gles,  especially  when  propagating  horizontally.  Since  light  travels  along  geodesics, 
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the  propagation  path  cuts  through  turbulence  at  a  lower  altitudes,  causing  stronger 
turbulence  in  the  middle  of  the  path.  To  consider  these  differences,  the  values  for  h  in 
the  turbulence  profiles  would  include  h(z)  as  a  function  of  the  propagation  distance 
location  2  [69].  For  the  propagations  and  profiles  modeled  in  this  work,  the  differ¬ 
ence  in  the  strength  of  turbulence  C^(z)  for  flat  earth  and  curved  earth  models  are 
negligible. 


2.2.3  Atmospheric  phase  parameters.  Diffraction  and  turbulence  signifi¬ 
cantly  affect  laser  beam  propagation.  Whether  the  beam  propagates  through  air  or 
through  space,  the  beam  diffracts  and  spreads  as  it  propagates.  Through  the  air,  ther¬ 
mal  mixing  of  the  air  causes  random  index  of  refraction  variations  leading  to  random 
phase  fluctuations  of  optical  beams.  These  random  fluctuations  cause  laser  beams  to 
spread,  randomly  wander,  and  scintillate. 


These  turbulence  effects  also  limit  the  resolution  of  imaging  systems.  Without 
atmospheric  turbulence,  the  resolution  of  an  imaging  system  improves  as  the  receiver 
aperture  increases,  commonly  referred  to  as  the  diffraction-limited  resolution.  In 
turbulence,  the  imaging  system  is  approximately  limited  to  the  diffraction-limited 
resolution  of  a  r0-diameter  aperture.  For  a  plane-wave  source,  this  coherence  diameter 
(a.k.a.  Fried  parameter)  r0tPW  is  [69] 


^0  ,pw 


L  l  -3/5 

2/ 


0.423/c2  /  Cn(z)  dz 
Jo 

2.  — \  —3/5 


(0.423/c2/io ) 


(27) 

(28) 


and  for  a  point  source  (spherical  wave),  r0)SW  is  defined  as 


^ 0,sw 


(29) 

(30) 
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where  k  =  2n/X  is  the  wave  number  and  A  is  the  optical  wavelength.  The  isoplanatic 
angle  #o  defines  the  angular  separation  between  beams  where  the  higher  order  phase 
effects  are  relatively  correlated.  This  angle  is  used  to  determine  the  maximum  angu¬ 
lar  separation  between  a  reference  beam  and  the  transmission  beam  to  measure  the 
higher-order  phase  effects.  It  is  defined  as  [69] 


en  = 


2.91  k2  /  C2n(z)(L  -  zf/3  dz 


-3/5 


(31) 


In  each  case  the  light  propagates  from  the  source  at  z  =  0  to  the  receiver  at  £  =  L  [27]. 


2.2.4  Atmospheric  beam  wander  and  beam  spreading.  In  free-space  laser 
communications  when  the  laser  beam  wanders  off  the  receiver  due  to  turbulence,  a 
fade  occurs.  For  constant  C2  (as  in  horizontal  propagation),  the  Gaussian-beam 
wander  was  characterized  by  Andrews  and  Phillips  by  the  second  moment  of  the 
beam’s  hot-spot  displacement  for  the  infinite  outer  scale  case  as  [2] 

(rc)  =  2A2C2nL3W-1/3  2F1  Q,  1; 4;  1  -  |0O|)  ,  (32) 

where  2F1  is  the  hypergeometric  function,  @0  is  defined  in  Eq.  (5).  The  equation  can 
be  simplified  for  a  collimated  beam 

(rc)  =  2.42  Cp3w~1/3  (33) 

=  0.0657  m2 

and  for  a  focused  beam 

(r2)  =  2.72C2L3Wq~1/3),  (34) 

=  0.0738  nr2, 
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where  Wq  is  the  effective  beam  radius  and  L  is  the  propagation  length.  The  numeric 
answer  given  above  is  for  an  air-to-air  100  km  path  with  C2  =  10-1'  m-2/3  for  an 
infinite  outer  scale  (i.e.  Kolmogorov  turbulence  model)  and  Wo  =  2.5  cm.  These  same 
parameters  are  used  in  simulations  in  subsequent  chapters.  For  the  uplink  channel, 
the  rms  beam  wander  is  [2] 

v^|>=0.73i/SecK)(^)  (^)  .  (35) 


where  £  is  the  zenith  angle. 

In  order  to  calculate  the  centroid  location  probability  density  function  (PDF), 
let  the  received  beam’s  centroid  Cartesian  components  be  xc  and  yc.  If  the  turbu¬ 
lence  is  isotropic,  one  can  assume  that  xc  and  yc  are  independent  and  identically 
distributed  random  variables  with  a  Gaussian  PDF.  The  PDF  of  rc  can  be  calculated 
by  first  determining  the  x  and  y  contributions  to  the  PDF  and  using  those  in  a  jointly 
Gaussian  PDF.  Since  the  x  and  y  centroid  variances,  a2c  and  a2Jc,  are  equal,  they  can 
be  determined  by 

=  yj  (36) 

resulting  in  a2c  =  a2c  =  (r2)  /2.  Using  the  x  and  y  centroid  variances,  the  centroid 
location  PDF  takes  the  form 


fxcyc  (*u 


i 

- - exp 

27TCr3;c,  <Jyc 


1 

27 nr2 


exp 


(  x2  +  y2\ 

\  2<  )■ 


(37) 

(38) 


For  the  100  km  C2  =  10  17  m  2/3  air-to-air  path  cr2^  =  cr,2c  =  (r2)  /2  =  0.033  m2  and 

1  /O 

the  standard  deviation  is  (r2)  '  =  0.256  m. 

1  /Q 

Now  computing  the  PDF  in  terms  of  r  =  ( x 2  +  y2)  7  requires  a  transformation 
of  the  random  variables  x  and  y  to  r  [57].  The  cumulative  density  function  (CDF)  is 
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determined  by 


Fr(r)  =  1  -  exp  i  r  >  0  (39) 

Taking  the  derivative  of  the  CDF  yields  the  PDF 

T 

frir)  =  —  exp 

where  a2c  =  a2Jc  =  (r2)  /2.  The  turbulence-induced  short-term  beam  radius  Wst 
at  the  receiver  combined  with  the  random  centroid  locations  determines  the  average 
power  at  the  receiver.  This  radius  for  a  collimated  beam  is  [2] 


2  a2x 

JjC 


r  >  0 


(40) 


/  A2  \  V6] 

WST  =  W . 

1  +  1.64C'^A;7/6L11/6A5/6 

r°je(i+A§)  | 

(41) 


while  for  a  focused  beam  the  radius  is  [2] 


WST  =  H'yl  + 


(42) 


where  A  =  z/  ( kW 2).  For  the  100  km  air-to-air  path  the  beam-spreading  factor 
Wst/W  does  not  change  much  with  Wq.  Figure  5  shows  the  maximum  beam¬ 
spreading  occurs  when  Wq  =  0.263  m  with  a  beam  spreading  factor  of  Wst/W  =  1.14. 
In  this  case,  for  any  diameter  beam,  most  of  the  beam  spreading  is  caused  by  diffrac¬ 
tion  not  turbulence. 


The  long-term  spot  size  Wlt  is  calculated  using  the  short-term  spot  size  and 
the  random  centroid  locations  to  determine  the  average  power  profile  at  the  receiver 
plane.  The  analytic  long-term  field  follows  a  Gaussian  shape  given  by 


U(x,y) 


WLTx 


1 


2tt  Wltx 

WLTy  = 


wLTy 

Wlt 

2 


exp 


-0.5 


(  x 2 

+ 


y 2  ^ 


(43) 

(44) 
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1.15 


Figure  5:  Turbulence  induced  short-term  beam  spread  for  the  100  km  air-to-air  path 
with  C2  =  10~17  m~2/3. 

where  [2] 


Wlt  =  yjwlr  +  (rl) 

=  wj  1  +  1.33ctJA5/«. 


(45) 

(46) 


A  Monte-Carlo  simulation  with  random  draws  of  centroid  locations  was  performed 
to  verify  Eq.  (46).  As  shown  in  Fig.  6,  the  simulation  matched  very  closely  to  the 
Gaussian  profile  and  the  theoretical  spot  size,  since  Wlt  is  determined  by  the  Wst 
and  (r2).  This  also  confirms  that  the  PDF  used  to  determine  the  centroid  locations 
is  consistent  with  the  definition. 


2.2.5  Beam  wandering  and  spreading  effects  on  laser  communication.  Next, 
the  beam  wandering  is  combined  with  the  beam  spreading  to  determine  the  effect 
on  the  received  power.  For  this  example,  the  receiver  has  a  circular  aperture  with  a 
radius  of  5  cm,  and  all  the  light  collected  by  the  aperture  is  detected.  The  cumulative 
probability  of  collected  power  is  used  to  determine  the  probability  of  fading.  The 
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Figure  6:  Long-term  optical  field  profile  for  a  Wo  =  2.5  cm  and  a  constant  C2  = 
1 0 ~ 1 '  m-2/3  air-to-air  100  km  path.  Plotted  using  short-term  spot  size 
and  random  centroid  locations  (using  rms  beam  wander)  in  a  1000  draw 
Monte-Carlo  simulation  and  compared  to  analytic  long-term  spot  size. 

highest  power  in  this  plot  corresponds  to  a  beam  received  on-axis.  Conversely,  the 
lower  received  powers  occur  as  the  beam  wanders  off  axis,  corresponding  to  a  fade  in 
a  laser  communication  link.  The  probability  is  determined  from  the  hot-spot  location 
PDF  and  the  beam  spreading  effects  mentioned  earlier. 

Taking  beam  spreading  and  wandering  effects  into  account,  it  is  suspected  a 
smaller  Wo  would  improve  signal  fading  for  the  long  propagation  path  of  the  air- 
to-air  100-km  path.  For  the  case  of  Wa  =  10  cm,  the  following  results  were  found 
W st  =  53  cm  and  (r2)  =  521  cm2.  The  beam  wander  for  Wq  =  10  cm  is  less  than  for 
the  Wo  =  0.05  cm  case,  although  the  received  beam  size  played  a  much  larger  factor. 
The  larger  the  transmit  aperture,  the  smaller  the  beam  size  at  the  receiver.  For  these 
smaller  received  beams,  the  beams  walk  off  the  receiver  and  deep  fades  occur.  For 
the  larger  received  beams  the  maximum  power  received  is  less,  but  the  fades  are  not 
as  deep  since  the  beam  does  not  completely  walk  off  the  receiver. 
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Figures  7  and  8  show  how  this  affects  the  collected  power  PDF  and  CDF.  For 
Wo  =  10  cm,  the  probability  that  the  received  power  was  less  than  a  -47  dBm  (150  pho¬ 
tons  at  1  GBit/sec)  threshold  was  over  1012  times  higher  than  the  Wo  =  5  cm  case 
and  1048  times  higher  than  the  Wq  =  2.5  cm  case.  To  achieve  a  fade  probability  of 
less  than  10  10,  the  W0  =  10  cm  transmitter  required  a  threshold  60  times  lower  than 
the  Wo  =  5  cm  case  and  360  times  lower  than  the  Wo  =  2.5  cm  case.  The  lower 
the  required  threshold,  the  less  the  margin  for  error,  leading  to  deeper  fades  at  the 
receiver. 

Another  approach  would  be  to  defocus  the  beam,  spreading  the  beam  out  at 
the  receiver  and  reducing  the  fades.  This  effect  can  be  illustrated  by  comparing  the 
results  for  a  collimated  beam  and  a  focused  beam.  For  this  long-propagation  example, 
the  difference  between  the  beam  size  at  the  receiver  is  still  significant  enough  to  affect 
the  results.  As  expected,  the  collimated  beams  performed  better,  resulting  in  less 
fading  than  the  focused  beam  as  shown  in  Figure  8  for  the  two  Wo  =  5  cm  cases. 

For  long-distance  propagations,  the  beam  width  steadily  increases  as  it  propa¬ 
gates  due  to  diffraction.  The  turbulence  also  spreads  the  beam  width  even  more,  but 
the  primary  cause  for  the  extended  turbulence  scenario  is  diffraction.  Narrow  beam 
widths  tend  to  wander  more  than  wide  beam  widths  since  the  average  tilt  across  these 
narrow  beams  is  higher.  Referring  back  to  Figure  2  on  beam  spread  due  to  diffrac¬ 
tion,  it  is  clear  that  if  the  strength  of  turbulence  varies  over  the  path,  this  beam  width 
affects  the  turbulence-induced  beam  wander  and  beam  spread.  Turbulence-induced 
phase  tilts  at  the  beginning  of  the  propagation  result  in  larger  beam  wander  due  to 
the  long  lever  arm  of  the  turbulence  and  small  beam  size  at  the  transmitter. 

One  phenomenon  not  considered  in  this  analysis  is  the  beam’s  angle  of  arrival 
(AOA)  variance.  Typically  the  receiver  will  have  a  large  collection  optic  that  focuses 
the  beam  onto  the  detector.  As  the  AOA  varies,  the  irradiance  on  the  detector  also 
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Probability  Density 


Received  Power  (mW)  Received  Power  (mW) 


(b)  (c) 

Figure  7:  100  km  Air-to-Air  Scenario.  Collected  power  PDF  in 
milliwatts  for  (a)  Wo  =  10  cm;  (b)  Wo  =  5  cm;  and 

(c)  Wo  =  2.5  cm. 
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Figure  8:  100  km  Air-to-Air  Scenario.  Collected  power  CDF  for 
Wq  =  10,  Wq  =  5,  and  Wq  =  2.5  cm.  A  cumulative 
probability  of  1  is  reached  at  the  maximum  collected 
power. 


varies.  The  Rx  AOA  variance  in  radians  squared  is  [69] 

6.08 


rr\A  _ 

t  ~  Dl/3L2 


z  Cn(z)dz, 


(47) 


where  D  is  the  receiver  diameter.  The  rms  image  jitter  is  just  the  focal  length  times 
the  rms  angle  of  arrival,  f'Tt,  which  determines  an  additional  received  power  variation. 
This  image  jitter  decreases  as  the  receiver  optic  size  increases,  since  the  variances  in 
the  phase  average  over  a  larger  optic.  Small  receiver  optics  may  encounter  large  image 
jitter  due  to  turbulence  and  usually  require  a  receiver  tracker.  The  more  the  incident 
light  is  angularly  off-axis,  the  less  power  is  received  at  the  detector. 


2.3  Amplitude  effects 

There  is  another  critical  effect  that  the  model  in  the  previous  section  of  beam 
spreading  and  wandering  does  not  include  that  must  be  considered  when  determin- 
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ing  signal  fading.  That  effect  is  beam’s  temporal  and  spatial  intensity  variance  (or 
scintillation).  Not  only  does  the  beam  spread  out,  it  also  varies  in  intensity  signif¬ 
icantly  over  the  entire  beam.  These  intensity  variations  increase  with  propagation 
distance  through  the  turbulence  and  are  more  significant  the  further  off  the  beam 
center.  These  variations  are  also  more  significant  for  wider  beams  than  for  narrower 
beams,  leading  to  more  pronounced  fluctuations  as  the  beam  widens  due  to  diffrac¬ 
tion.  To  determine  the  optimum  beam  size,  diffraction  and  the  turbulence  effects  of 
beam  wander,  beam  spread,  and  scintillation  must  be  considered. 

As  the  beam  propagates,  the  phase  aberrations  cause  intensity  variations  at  the 
receiver.  These  effects  are  particularly  important  in  a  laser  communication  system, 
since  these  variations  can  cause  fades  and  bit  errors.  This  turbulence  effect  can  be 
characterized  by  the  log-amplitude  variance,  often  referred  to  as  the  Rytov  number 
1Z,  as  it  is  calculated  using  the  Rytov  approximation.  The  Rytov-theory  expressions 
for  the  scintillation  discussed  here  are  only  valid  for  log-amplitude  variances  less  than 
0.25.  Beyond  this  value,  the  scintillation  begins  to  saturate.  This  saturation  value  is 
typically  between  0.3  and  0.4  [69].  The  plane-wave  log-amplitude  variance  cw  at  the 
receiver  is  [2] 


(48) 


and  the  spherical-wave  log-amplitude  variance  is 


(49) 


Specifically,  for  a  plane-wave  space-to-ground  propagation  [2] 


0.5631fc7/y5/e. 


(50) 


while  for  ground-to-space  propagation 


H 


(51) 
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and  for  horizontal  propagation  with  a  uniform  C \  profile 


0.3071  k7/GC2nL11/6. 


(52) 


The  spherical  wave  Rytov  number  is 


nsph  =  0.5631  k7/6 


(53) 

(54) 


nsph  <  0.25. 


For  weak  turbulence,  the  spherical  wave  Rytov  number  7 Zsph  is  equal  to  the  log- 
amplitude  variance  cU  sw.  Since  the  scintillation  saturates  with  strong  turbulence,  the 
Rytov  number  does  not  equal  the  log-amplitude  variance  for  1Zsph  >  0.25. 

The  detector  at  the  receiver  ultimately  measures  the  irradiance.  Under  weak 
turbulence,  the  normalized  irradiance  variance  (referred  to  as  the  scintillation  index, 
i.e.  the  irradiance  variance  over  the  average  irradiance  1/(1)),  is  [2] 


(C(r))  (/(r  ))*  (/2(r)) 


(55) 

(56) 

(57) 

(58) 


UM>2  UM)2  (/(r)>2 


4cr^(r)  for  <  0.25 

i (r )  for  <  0.25, 


where  r  is  the  radial  distance  from  the  optical  axis  and  afn  r(r)  is  the  log- irradiance 
variance.  Given  this  relationship,  the  plane-wave  scintillation  index  is 


L 


L 


a/nl  =  2.251  k7/6 


Cl(z)(L-zf/edz, 


(59) 


and  for  spherical-wave  propagation,  it  is 


l 


■L 


a/n  j  =  2.251  k7/6 


C2(.-)(2/I)i,S(I-l)S,‘d2. 


(60) 


Andrews  and  Phillips  define  the  plane-wave  scintillation  index  for  a  constant  C2  profile 
as  the  Rytov  variance  4  [2],  From  the  above  relations,  the  Rytov  variance  is  defined 
as 

4  =  1 .23fc7/6C^L11/6.  (61) 

For  example,  the  Rytov  variance  for  the  air-to-air  100  km  path  for  C2  =  10“ 1 '  m_2//3 
is  a\  =  0.924. 

2.3.1  Strong  fluctuation  theory.  The  previous  section’s  relations  for  am¬ 
plitude  and  intensity  variations  were  valid  for  weak  turbulence  only.  They  were 
computed  through  Rytov  theory.  In  strong  turbulence,  the  phase  variance  contin¬ 
ues  to  increase,  while  the  scintillation  or  variance  in  the  amplitude  and  irradiance 
saturates.  There  are  multiple  theories  used  to  predict  the  behavior  as  the  Rytov 
number  increases  [2],  Extended  Rytov  theory  describes  the  scintillation  by  breaking 
it  up  into  large-  and  small-scale  amplitude  variations.  These  irradiance  variances  use 
the  gamma-gamma  distribution  to  describe  the  irradiance  variation  PDF.  The  irra¬ 
diance  variance  due  to  turbulence  for  spherical,  planar,  and  Gaussian  beams  can  be 
described  by  [2] 

ct2i(L)  =  exp  (4  x  +  afnY)  -  1,  (62) 

where  ofn x  and  crfnY  are  the  variances  in  the  small-  and  large-scale  irradiance  fluctu¬ 
ations,  respectively.  The  equations  for  these  variances  differ  depending  upon  the  type 
of  source.  In  the  limiting  case  of  Kolmogorov  turbulence  (i.e  /0  =  0  and  L0  =  oo)  for 
spherical  waves,  the  large  and  small-scale  log  variances  are  represented  by  [2] 
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For  planar  wavefronts  these  equations  become  [2] 
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2.3.2  Temporal  effects.  To  represent  the  temporal  statistics  accurately, 
the  temporal  frequency  of  the  turbulence  must  be  modeled  appropriately.  The  two 
quantities  used  to  describe  the  turbulence  frequency  are  the  Greenwood  frequency  fa 
and  the  Tyler  frequency  fr-  They  are  given  by  [66] 
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respectively.  These  are  the  temporal  power  spectrum  “break  frequencies.”  A  majority 
of  the  tilt  jitter  is  below  the  Tyler  frequency.  Most  of  the  frequencies  for  the  higher- 
and  lower-order  phase  disturbances  are  below  the  Greenwood  frequency. 


The  wind  velocity  V{h)  in  (21)  is  frequently  described  by  the  Button  wind 
model  [2] 


V  ( h )  =  a >sh  +  Vg  +  30  exp 


fh  —  9400  \  2 
v  4800  J  ’ 


(71) 
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where  Vg  is  the  ground  wind  speed  and  us  is  the  slew  rate  of  a  ground  transmitter 
(Tx)  tracking  an  aerial  receiver  (Rx). 


2-4  Layered  atmosphere  model 

Simulations  are  often  used  to  study  atmospheric  turbulence  effects  on  imaging, 
communication,  and  beam-projection  systems.  Analytic  solutions  to  wave  propaga¬ 
tion  oftentimes  require  assumptions  that  limit  the  validity  of  the  results,  like  the 
Rytov  approximation  that  limits  the  turbulence  to  weak  fluctuations.  Not  only  that, 
adaptive  optics  (AO)  systems  and  diversity  techniques  such  as  multiple-transmitter 
performance  cannot  be  calculated  in  closed  form.  The  turbulence  effects  in  these 
simulations  can  be  represented  by  optical  field  screens  placed  along  the  propagation 
path  having  the  appropriate  statistics. 

If  the  layered  model  matches  the  refractive  index  spectrum  and  the  phase  vari¬ 
ance  of  the  propagation  path,  it  can  be  used  in  analytic  calculations,  computer  sim¬ 
ulations,  and  experiments.  The  layers  are  chosen  to  represent  the  continuous  model 
as  closely  as  possible  so  that  several  low-order  moments  of  the  layered  model  match 
the  continuous  model 
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where  N  is  the  number  of  phase  screens  being  used  and  C^.  is  the  effective  turbu¬ 


lence  strength,  A Z{  is  the  turbulence  layer  thickness,  and  is  the  location  of  the  ith 


screen  [66].  The  atmospheric  parameters  for  the  layered  turbulence  model  are 
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which  are  analogous  to  the  continuous  model  versions  in  Eqs.  (28),  (30),  (48),  (49), 
and  (31).  The  full  path  r0  can  be  determined  by  the  phase  screens’  r0.  along  the  path 
by 

N 
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The  screen  r0i  value  is  related  to  C2  by 


t o//3  =  0.423*2Cj(2i)A2j. 


(79) 


Equations  (73)  -  (77)  can  be  written  in  terms  of  r0.,  which  simplifies  the  method  of 
choosing  screen  properties  and  screen  locations. 

To  determine  the  temporal  performance  for  layered  models  in  terms  of  the 
Greenwood  and  Tyler  frequencies  use 
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and 


1/2 


respectively. 
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Many  methods  have  been  developed  to  generate  phase  screens.  The  methods 
can  be  broken  into  two  mathematical  approaches.  The  first  approach  represents  the 
wavefront  phase  with  a  2-dimensional  rectangular  grid  of  points  —  a  sampled-based 
approach.  In  the  second  approach,  the  screen  is  represented  as  a  sum  of  orthogonal 
basis  functions  or  a  modal  basis.  The  most  common  method  is  the  sample-based 
screen  representation  in  which  the  fast  Fourier  transform  (FFT)  is  used  to  compute 
the  screen  realizations.  The  FFT  technique  is  direct  and  computationally  efficient. 
The  problem  with  FFT  based  methods  is  that  the  energy  of  low  spatial  frequencies  in 
the  screen  is  under-represented.  The  low  spatial  frequencies  (e.g.,  tilt)  contain  a  large 
fraction  of  the  power  for  atmospheric-turbulence-induced-wavefront  perturbations.  A 
modal-based  approach  allows  for  much  better  low  frequency  representation  [48,66,83]. 
Other  FFT  techniques  with  low-order  boost  (such  as  sub-harmonic  or  Zernike)  work 
well  for  static  simulations  but  large  screens  must  be  created  (and  stored)  for  long  time 
series  simulations  due  to  discontinuities  at  the  edges. 

The  Fourier-series  (FS)  phase  screen  generation  approach  allows  flexibility  in 
choosing  the  lowest  and  highest  spatial  frequencies  sampled  independent  of  the  phase 
screen  grid  sampling.  The  modal  method  used  here  follows  Welsh’s  development  based 
on  the  FS  [83].  A  FS  expansion  of  the  wavefront  phase  over  a  square  area  is  used 
as  the  basis  for  representing  the  phase  screen.  This  approach  much  more  accurately 
represents  the  low  spatial  frequencies  than  the  sample-based  approach.  The  modal 
phase  screen  is  defined  for  all  space  and  need  only  be  evaluated  at  the  grid  points  of 
interest.  Shifting  the  phase  screen  to  non-integer  multiples  of  samples  is  as  easy  as 
evaluating  the  FS  on  the  shifted  grid,  using  the  same  FS  coefficients.  This  approach 
is  particularly  effective  since  the  fields  of  the  multiple  widely  spaced  beams  must  be 
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calculated  over  long  periods  of  time.  Here,  logarithmically  spaced  frequencies  are 
used,  as  recommended  in  an  MZA  report  by  Magee  [50]. 

The  FS  expansion  can  be  approximated  by  truncating  the  summations  to  a 
finite  number  of  terms.  The  truncated  expansion  for  the  phase  is  [83] 
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where  the  phase  is  represented  in  a  square  of  dimension  Dp  and  c^  n,  is  the  FS  coeffi¬ 
cient  for  the  spatial  frequency  f  =  xnj Dp  +  yn' / Dp.  The  terms  x  and  y  are  the  x— 
and  y— directed  unit  vectors,  and  x  and  y  are  the  components  of  the  spatial  vector  x. 
The  continuous  model  atmospheric  parameters  are  matched  to  the  layered  turbulence 
model  for  the  isoplanatic  angle  6q,  Rytov  number  Tl,  and  coherence  diameter  r0. 


2-4-1  PDF  of  the  scintillation.  The  irradiance  variations  at  the  receiver, 
also  called  scintillation,  include  temporal  variations  as  well  as  spatial  variations,  but 
for  now,  only  the  spatial  variations  are  considered.  In  Section  5.4.2,  the  PDFs  are 
calculated  for  different  scenarios  in  temporal  simulations.  Many  atmospheric  turbu¬ 
lence  experiments  under  weak  turbulence  conditions  have  shown  that  the  natural  log 
of  the  intensity  is  approximated  by  a  Gaussian  distribution,  i.e.  [56] 
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Above,  l  =  In (//(/))  is  the  log  intensity,  /  is  the  optical  intensity  at  a  point,  and 
a2  is  the  variance  of  the  log  intensity.  After  performing  a  transformation  of  random 
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variables,  Eq.  (84)  takes  on  the  familiar  lognormal  PDF 


/(/)  =  — L=^exp[-  [ln(//(/))  +  o-f/2] 2 /2cr|] ,  (85) 

V/2vra^  J 

where  the  log-intensity  mean  is  (f)  =  —  crf/2.  The  log  intensity  variance  follows 
directly  from  the  log-amplitude  variance  which  depends  upon  the  strength  of  the 
optical  turbulence  along  the  path.  Therefore,  for  a  uniform  horizontal  path,  as 
shown  earlier,  the  on-axis  log  intensity  variance  is 

„  _  _  f  1.23  C'^7/6L11/6,  plane-wave 

ali  =  al  =  Aal=\  n  .  (86) 

0.496C^fc7//6Ln/6,  spherical  wave 

The  intensity  variations  across  the  receiver  pupil  were  characterized  in  Andrews  and 
Phillips’  book  by  the  following  general  expression  for  the  scintillation  index  of  a  Gaus¬ 
sian  beam  [2]  : 

crh(r,L)  =  3.62£;7/6L5/6A5/6-^-  [  C^(z)(l  —  zjVf^  dz,  r<W  (87) 

IT  2  J  o 

in  general  and 

2 

r,L)  =  l.ll^A5'6^,  r<W  (88) 

for  a  constant-^  path.  It  is  clear  from  these  expressions  that  scintillation  increases 
with  distance  from  the  beam  center. 

Another  PDF  used  to  model  the  variance  in  the  intensity  due  to  stronger  atmo¬ 
spheric  turbulence  conditions  is  the  gamma-gamma  PDF.  The  gamma-gamma  PDF 
was  developed  under  the  assumption  that  large-  and  small-scale  irradiance  fluctua¬ 
tions  are  governed  by  the  following  gamma  distributions  [2] 
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where  a  and  (5  are  parameters  related  to  the  large-  and  small-scale  scintillations.  They 
are  determined  by 
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where  o2nX  and  o2nY  relate  to  the  intensity  variance  caused  by  the  large-  and  small- 
scale  irradiance  fluctuations  shown  in  Eqs.  (63)  -  (66). 

For  the  varying  turbulence  profiles  of  the  air-to-ground  and  ground-to-air  link, 
the  log-intensity  variance  afnI  can  be  substituted  for  the  Rytov  variance  oR.  Putting 
it  all  together,  the  gamma-gamma  PDF  for  the  irradiance  fluctuations  is  [2] 
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where  the  on-axis  mean  irradiance  (J(0,L))  7^  1  and  Kv(x )  is  the  modified  Bessel 
function  of  the  first  kind. 


2. 5  The  PDF  of  intensity  for  multiple  transmitters 

What  happens  to  the  irradiance  distribution  if  multiple  beams  are  received 
through  independent  turbulence  realizations?  It  is  a  reasonable  hypothesis  that  there 
should  be  a  significant  variance  reduction  and  ultimately  a  much  lower  bit  error  rate. 
There  is  no  closed-form  solution  for  the  distribution  of  the  sum  of  multiple  lognormal 
distributions  or  for  even  identically  distributed  ones.  If  these  multiple  laser  sources 


are  sufficiently  separated  they  could  be  considered  independent  identically  distributed 
(iid)  random  variables  (RV)  in  the  limiting  case.  A  number  of  folks  have  attempted 
to  estimate  the  distribution  of  the  sum  of  lognormal  variables  and  in  most  cases  have 
used  a  lognormal  PDF  in  the  estimate.  Fillio  et  al.  presented  a  very  accurate  estimate 
by  matching  the  first  two  moments  of  the  inverse  of  the  exact  sum  with  those  of  the 
lognormal  approximation  inverse,  i.e.  [22] 


E  [A-] 

=  E  [S_1]  and 

(94) 

E  [A"2] 

=  E[S“2], 

(95) 

where  S  is  the  sum  of  Al  lognormal  variables  and  X  is  the  proposed  approximation  to 
S  and  E[-]  is  the  expectation  operation.  Using  the  fact  that  the  lognormal  distribution 
kth  moment  is  given  by 

E  =  e^+fc2ff2/2,  (96) 

they  solve  for  /i  and  cr,  resulting  in 

Us  =  0.5  In  E  [S'-2]  —  21nE  [S'-1]  (97) 

<x2  =  lnE  [S~2]  -21nE  [S'"1]  .  (98) 

This  estimate  is  accurate  for  moments  much  higher  than  the  second  moment.  The 
problem  with  this  approach  is  there  is  no  closed-form  solution  to  the  moments  E  [S^1] 
and  E  [S~2] . 

A  straightforward  approach  to  determine  the  distribution  of  the  sum  of  iid 
signals  up  to  the  second  moment  can  be  determined  since  the  sum  of  lognornrals 
follows  an  approximately  lognormal-type  distribution.  For  the  sum,  s  =  YliLi  xi 
of  independent  signals  xt,  the  variance  of  s  is  the  sum  of  the  individual  variances 
VAR  [ay] ,  and  the  mean  is  the  sum  of  the  individual  means  E[xj],  where  VAR  [■]  is 
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the  variance  operation.  That  is  to  say  [57] 
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and  when  they  are  independent  and  identically  distributed  (i.e.  iicl)  they  become 


E  [s]  =  M-E[Xi] 

VAR  [s]  =  M  ■  VAR  [xt]  . 

Lognormally  distributed  independent  variables  have  the  PDF  [57] 
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E  [Xi]  =  e(Mi+0'*?/2) 


VAR[xi]  =  (ea?  -  l)  e2w+< 

Therefore,  the  mean  and  variance  of  the  sum  of  M  iid  variables  are 
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and 

VAR  [s]  =  M  -  l)  e2fli+<Ti .  (107) 

Assuming  the  approximation  of  the  sum  of  lognormally  distributed  independent  vari¬ 
ables  is  also  lognormally  distributed,  the  parameters  /is  and  as  can  be  represented  in 
terms  of  the  E  [s]  and  VAR  [s]  .  Using  the  general  form  of  a  lognormal  distribution, 
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the  parameters  /is  and  as  are  represented  by 
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Therefore,  in  general  for  s  =  'Yh^Lixii  where  X;  are  lognormally  distributed  inde¬ 
pendent  variables  with  parameters  /Xj  and  of,  the  PDF  of  s  can  be  approximated 
as 
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For  iid  signals,  the  parameters  /i.s  and  a2  can  be  further  simplified  to 
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As  M  gets  large,  the  distribution  approaches  a  Gaussian  distribution,  as  it  should 
according  to  the  central  limit  theorem  (CLT)1.  The  CDF  of  a  lognormal  RV  is 
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Similarly,  the  CDF  of  the  sum  of  lognormal  variables  is  given  by 
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where  /is  and  as  are  defined  in  Eq.  (Ill)  and  Eq.  (112). 

The  problem  with  this  approach  is  that  it  is  only  accurate  for  the  first  two 
moments  of  the  lognormal  PDF  and  breaks  down  for  higher-order  moments.  These 
higher-order  moments  tend  to  affect  the  tails  of  the  distribution.  As  seen  earlier 
for  the  beam  wander  and  beam  spreading  example  in  Section  2.2.5,  the  tails  of  the 
distribution  significantly  affect  the  cumulative  distribution  used  to  calculate  the  fade 
probability. 


lrThe  central  limit  theorem  states  that  the  distribution  of  the  sum  of  a  sufficient  number  of 
independent  RVs  with  finite  mean  and  variance  is  approximately  normally  distributed.  The  CLT 
explains  why  the  Gaussian  RV  adequately  describes  so  many  different  natural  random  processes  [57, 
74], 
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Going  back  to  the  fundamental  approach,  the  PDF  for  the  sum  of  M  indepen¬ 
dent  RVs  can  be  determined  by  the  convolution  of  the  PDFs  of  the  individual  PDFs, 
noted  by  [57] 

fa  =  fl  ®  fl  ©  '  '  ’  ©  fn,  (121) 

where  the  convolution  operation  symbol  is  ©.  The  convolution  can  be  calculated 
by  Fourier  transforming  the  PDFs,  multiplying  the  resultant  functions  and  inverse 
transforming  the  result  as  follows 

ps  =  pi  ®  P2  ®  ■  ■  ■  ®  Pm  (122) 

F{Pi}  J,  (123) 

where  J-  is  the  Fourier  transform  operator  and  T  1  is  the  inverse  Fourier  transform 
operator.  Since  there  is  no  closed-form  solution  to  the  Fourier  transform  of  the  log¬ 
normal  PDF,  one  can  use  the  discrete  Fourier  transform  (DFT)  to  approximate  the 
result,  where  Tkn  is  the  DFT  operator  and  T  1  is  the  inverse  DFT  operator.  Using 
the  DFT,  the  distribution  can  be  approximated  by 

/  1  \M_1  f  M  1 

Ps  (xk)  =  (  A;  J  T~kn  \PiM]>  ,  (124) 

where  Pi{;xk )  are  the  discretized  independent  PDFs  (sampled  N  times)  of  the  multiple 
transmitters  and 
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The  DFT  operation,  Thnl  and  inverse  DFT  operation,  Tk  1,  are  dehned  by 
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and 


1  7V”1 

KS-}  =  <  <127> 

n= 0 

respectively.  The  finer  the  sampling  of  the  individual  PDFs,  the  more  accurate  the 
PDF  estimate  of  the  sum.  This  relation  applies  to  independent  random  variables, 
assuming  the  Txs  are  adequately  spaced  to  be  independent.  The  intensity  distribution 
of  each  independent  laser  path  varies  depending  upon  how  far  off  axis  they  are  located. 
As  M  gets  large,  the  distribution  approaches  a  Gaussian  distribution,  as  it  should, 
according  to  the  CLT. 

Due  to  the  periodic  nature  of  the  DFT,  any  aliasing  must  be  mitigated.  In  this 
case,  one  needs  to  ensure  there  is  finite  support  to  the  PDF  and  that  the  PDF  is 
“zero-padded”  so  that  the  PDF  extends  to  at  least  2  times  the  support,  to  ensure 
minimal  wrap-around. 

To  test  the  accuracy  of  this  approach,  the  PDF  is  calculated  for  the  sum  of  inde¬ 
pendent  lognormal  variables  in  a  Monte-Carlo  simulation  and  compared  to  the  DFT 
approach  mentioned  here.  This  was  accomplished  by  generating  5  independent  log¬ 
normal  random  variables  (RVs),  adding  them  together,  and  calculating  the  resulting 
PDF.  In  Fig.  9  the  average  intensity  of  the  sum  of  the  lognormal  variables  is  compared 
to  a  single  lognormal  variable.  The  means  of  the  lognormal  variables  used  in  the  sum 
were  1/5  the  mean  of  the  individual  lognormal  variable.  The  plot  shows  that  for  the 
sum  of  the  RVs  the  mode  is  greater  and  the  variance  has  been  reduced  as  well.  As 
M  gets  very  large  the  PDF  of  the  sum  looks  Gaussian  with  a  mean  centered  around 
1/(1)  =  e.  In  most  practical  cases  the  turbulence  effects  of  two  or  more  beams  are 
partially  correlated.  The  less  correlated  the  effects,  the  more  these  multiple  beams 
will  reduce  signal  variability. 

This  DFT  approach  also  works  with  the  gamma-gamma  PDF.  Here,  the  prob¬ 
ability  of  a  fade  for  a  100  km  FSOC  system  operating  a  horizontal  link  with  a 
C'i  =  10-'7  m  2//3  is  calculated.  The  transmitter  uses  a  W0  =  2  cm  collimated 
beam  operating  at  a  wavelength  A  =  1.55  /ini  and  a  receiver  aperture  Dq  =  10  cm  for 
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Figure  9:  The  average  sum  of  lognormal  RVs  is  compared  with 
a  single  lognormal  variable,  of  =  1,  fii  =  1/2,  M  =  5 
RVs,  and  with  the  PDF  sampled  N  =  10000  times. 

the  Kolmogorov  turbulence  spectrum.  The  total  transmitted  power  of  the  single-Tx 
and  five-Tx  systems  were  identical.  A  fade  in  this  example  is  defined  as  when  the 
signal  drops  6  dB  below  the  mean.  The  probability  of  a  fade  for  a  single-Tx  system  is 
0.019.  Assuming  independent  transmitters,  the  probability  of  a  fade  for  the  five-Tx 
case  is  approximately  6.6  x  10  7,  reducing  the  fade  probability  by  a  factor  of  over 
28,000.  These  results  are  shown  graphically  in  Fig.  10  where,  distinctly,  the  single 
transmitter  has  much  higher  probability  of  being  to  the  left  of  the  threshold  than  the 
five-Tx  case.  Clearly,  multiple-Tx  systems  can  provide  significant  improvement. 

2.6  Digital  communication  and  detection  theory 

2.6.1  Modulation.  With  any  communication  system,  the  signal  at  the  trans¬ 
mitter  (Tx)  must  be  modulated  (i.e.  varied  in  some  way)  to  efficiently  represent  the 
information  and  propagate  it  effectively  through  the  channel  to  the  receiver  (Rx).  At 
the  Rx,  the  signal  is  received,  demodulated,  and  detected  to  extract  the  information 
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Figure  10:  PDF  of  a  single  gamma-gamma  RV  (solid  line)  and 
a  PDF  of  a  sum  of  5  gamma-gamma  RVs  (dashed 
line).  The  gamma-gamma  RV  models  the  received 
intensity  for  a  100  km  C2  =  10-1'  m-2/3  path  for  a 
Wq  =  2  cm-Tx  and  a  D  =  10  cm-diameter  Rx.  The 
threshold  Ft  is  a  solid  vertical  line  6  dB  below  the 
mean.  Multiple  Txs  are  scaled  to  one  fifth  the  power 
of  the  single  Tx. 
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from  the  signal.  The  information  can  be  represented  by  binary  (2  symbols)  or  M-ary 
(M  symbols)  modulation.  The  information  can  be  encoded  on  the  frequency,  polarity, 
signal  level,  phase,  the  pulse  width  of  the  signal,  etc.,  or  with  any  combination  of  these 
used  to  represent  2  to  M  symbols.  This  research  effort  uses  the  binary  modulation 
scheme  of  on/off  keying  (OOK),  which  is  accomplished  by  turning  on  the  laser  to 
transmit  a  ‘1’  and  turning  off  the  laser  to  transmit  a  ‘O’. 

2.6.2  Noise  sources  and  measured  signal  probability  density  functions.  The 
process  of  measuring  the  signal  at  the  receiver  is  inherently  noisy.  That  is,  there 
are  multiple  factors  that  contribute  to  the  measured  signal  uncertainty  at  the  Rx 
for  each  symbol.  Each  symbol  transmitted  has  a  different  PDF  associated  with  the 
measured  receiver  signal.  The  Gaussian  distribution  is  often  used  to  model  system 
noise  because  of  the  CLT,  defined  in  footnote  1.  Even  if  the  individual  noise  sources 
are  not  Gaussian,  often  the  total  noise  can  be  approximated  as  Gaussian  [74], 

There  are  many  potential  noise  sources  in  the  optical  measurement  process, 
such  as  electronic  thermal  (or  Johnson)  noise,  photon  (or  shot)  noise,  generation 
recombination  noise,  1/f  noise,  background  noise,  dark  current,  amplifier  noise  etc.  [18, 
19].  More  detail  is  given  in  Section  4.3.  For  optical  receivers  there  are  typically 
three  sources  of  noise  that  dominate:  thermal  noise,  shot  noise,  and  amplifier  noise. 
Thermal  noise  is  caused  by  thermal  motion  of  the  electrons  in  the  wires  and  resistor 
in  a  system.  Shot  noise  is  caused  by  the  random  arrival  of  photons  at  the  receiver  and 
is  characterized  by  a  Poisson  distribution,  which  can  be  approximated  by  a  Gaussian 
distribution  if  the  signal  levels  are  large  enough.  Inherently  there  are  noise  sources  in 
amplifiers,  and  optical  amplifiers  are  no  exception.  In  the  case  of  erbium  doped  fiber 
amplifiers  (EDFAs),  they  exhibit  amplified  spontaneous  emission  (ASE)  noise  [1]. 

2.6.3  Detection  theory.  For  the  binary  OOK  modulation,  the  laser  turns  on 
to  transmit  a  ‘1’  and  turns  off  to  transmit  a  ‘O’.  The  transmission  of  a  ‘1’  or  ‘0’  is 
equally  likely  and  denoted  by  the  events  Hi  and  H0j  respectively.  The  likelihood  ratio 
test  (LRT)  determines  the  optimal  decision  threshold  based  upon  the  PDFs  of  the 


45 


measured  signal  level  irn  for  the  transmission  of  a  ‘1’  p{im\Hi)  and  the  transmission 
of  a  ‘O’  p(im\H0).  In  general,  the  likelihood  ratio  test  is  [82] 


A  (im) 


pick  Hi 

pjimlHr)  >  P(H0) 

p(tm\H0 )  <  P{HiY 

pick  H0 


(128) 


For  the  equally  likely  case,  since  the  a  priori  probabilities  are  P(Hq)  =  P(Hi)  = 
0.5,  the  above  inequality  simplifies  to 


if  p(im\Hi)  >  p(im\H0)  pick  1/, 


(129) 


or 

if  p(im\H0)  >  piimlHx)  pick  H0.  (130) 

The  optimum  detection  criteria  can  best  be  described  graphically.  Given  equally 
likely  signaling,  the  optimum  detection  criterion  is  simply  the  point  at  which  the 
two  probability  densities  intersect  as  shown  in  Fig.  11.  This  maximum-likelihood 
solution  minimizes  the  total  probability  of  an  error.  The  PDF  of  the  receipt  of  a  ‘0’  is 
relatively  constant,  whereas  the  PDF  for  the  receipt  of  a  ‘1’  depends  upon  the  channel 
conditions.  This  channel  could  be  affected  by  an  obstruction,  the  weather  (clouds, 
rain,  fog,  etc.),  or  changes  in  atmospheric  turbulence  conditions.  These  turbulence 
conditions  vary  significantly  over  time  and  thus  could  benefit  from  a  threshold  that 
varies  with  the  optical  signal  level  [13,18].  The  optimal  fixed  threshold  is  determined 
in  Section  5.2.2. 1  and  implemented  in  Section  5.4.  The  optimal  adaptive  thresholds 
are  determined  in  Section  5. 2. 2. 2  and  implemented  in  Section  5.4. 

This  chapter  described  how  FSOC  systems  are  hampered  by  atmospheric  tur¬ 
bulence.  The  next  chapter  summarizes  the  state-of-the-art  in  FSOC  research  and  lays 
the  groundwork  for  the  research  described  in  Chapters  IV  through  VI. 
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Photocurrent 


Figure  11:  Decision  threshold  for  an  00K  digital  receiver.  If 
a  signal  is  detected  above  the  threshold  it  is  called 
a  ‘1’  and  below  the  threshold  it  is  called  a  ‘O’.  The 
shaded  regions  indicate  the  probability  of  an  error. 
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III.  Previous  Work 


FSOC  research  is  broad  and  extensive,  but  it  rarely  focuses  on  airborne  scenarios. 

Small-size  and  low-weight  constraints  are  imperative  for  making  airborne  laser 
communication  feasible.  The  key  to  this  research  is  to  use  an  integrated  approach 
to  engineer  systems  with  as  minimal  hardware  complexity  while  achieving  acceptable 
performance.  At  each  step  along  the  way,  the  technique  is  optimized  to  address  the 
appropriate  system  effect.  This  chapter  surveys  current  free-space  laser  communica¬ 
tion  research  and  identifies  opportunities  for  the  greatest  improvement  for  the  least 
cost  in  size,  weight,  and  complexity. 

Sections  3.1  -  3.3  describe  Tx  and  Rx  designs  that  increase  Rx  power  and  de¬ 
crease  power  fluctuations  in  atmospheric  turbulence.  In  Section  3.4,  modulation  tech¬ 
niques  are  identified  with  the  goal  of  finding  one  that  is  less  sensitive  to  low  power  and 
large  fluctuations.  Finally,  Section  3.5  delineates  different  signal  processing  techniques 
to  determine  how  to  decipher  the  bits  more  accurately. 

3.1  Optical  transmitter  and  receiver  design 

Transmitter  and  receivers  need  to  be  designed  for  a  particular  scenario  and  range 
of  atmospheric  conditions.  Much  of  the  research  has  been  focused  on  the  ground- 
to-space  and  the  ground-to-ground  propagation  paths,  rarely  covering  the  air-to-air 
scenario. 

Traiismitter  size.  The  transmitter  aperture  size  and  divergence  drives  the  size 
of  the  beam  at  the  receiver  due  to  diffraction  and  plays  a  large  part  in  the  type 
and  degree  of  intensity  variations  caused  by  atmospheric  turbulence.  Atmospheric 
turbulence  causes  the  beam  to  wander,  spread  out,  and  break  up,  resulting  in  intensity 
variations  at  the  receiver.  The  larger  the  beam,  the  less  the  beam  wanders,  but  the 
more  the  beam  breaks  up  and  the  lower  the  average  received  power.  Churnside 
determined  a  relation  for  beam  wander  in  weak  turbulence  [16].  Yenice  and  Evans 
proposed  adaptive  adjustment  of  the  beam  size  to  reduce  the  intensity  fluctuations 
on  a  ground-to-satellite  laser  uplink  [88].  They  mentioned  a  beam  size  adjustment  of 
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no  more  than  a  factor  of  two  conld  be  feasible  and  effective  [88,89].  This  technique 
requires  real-time  knowledge  of  the  atmospheric  conditions,  but  one  lesson  from  their 
research  is  to  properly  design  and  optimize  the  transmit  beam  size  for  the  type  of 
turbulence  most  likely  encountered  in  the  link.  A  transmitter  could  be  defocused  to 
increase  the  beam  size  at  the  receiver  thereby  reducing  signal  fades  due  to  beam  walk- 
off,  but  as  the  beam  size  increases  the  intensity  fluctuations  of  the  beam  also  increase. 
As  always,  a  trade-off  must  be  made  between  beam  walk-off,  average  received  power, 
and  intensity  fluctuations  [75].  Using  the  analysis  accomplished  in  Section  2.2.5  and 
Yenice’s  work  on  adaptive  Tx  beam  size  [88,89],  a  Wo  =  2.5  cm  collimated  source 
optimized  for  the  100  km  air-to-air  path  is  used  for  the  individual  Tx  size,  which 
meets  requirements  for  smaller  and  lighter  weight  components. 

Partially  coherent  beams.  The  intensity  fluctuations  of  an  optical  field  after 
propagating  through  turbulence  depend  upon  the  turbulence  along  the  path  and  the 
spatial  coherence  of  the  source.  Therefore,  one  approach  to  reduce  these  fluctuations 
is  to  decrease  the  spatial  coherence  of  the  beam  by  placing  a  phase  diffuser  in  front  of 
the  transmitting  beam  [64,65].  Reducing  the  coherence  of  the  beam  not  only  spreads 
out  the  beam,  but  also  reduces  the  intensity  fluctuations  caused  by  the  interference 
(often  referred  to  as  laser  speckle)  of  the  coherent  beam.  This  technique  is  very  effec¬ 
tive  for  short  propagation  paths,  but  for  long  propagation  paths  the  beam  spreading 
and  power  reduction  can  be  quite  severe.  Therefore,  it  is  not  appropriate  for  the 
air-to-air  scenario. 

Optical  amplification.  An  Erbium-doped  fiber  amplifier  (EDFA)  could  reduce 
the  power  required  at  the  detector  [29].  As  with  any  amplifier,  it  must  be  analyzed 
and  optimized  so  as  not  to  amplify  the  noise  level  so  much  that  it  becomes  counter¬ 
productive.  Razavi  and  Shapiro  studied  the  link  margin  improvement  of  an  amplified 
spontaneous  emission  (ASE)  optical  preamplifier,  varied  the  receiver  size  and  num¬ 
ber  of  receivers,  and  tested  OOK  and  binary  pulse  position  modulation  schemes  with 
adaptive  and  constant  thresholds  [63] .  They  showed  optical  pre-amplification  of  30  dB 
was  sufficient  to  overcome  thermal  noise,  increasing  the  receiver  sensitivity  by  over 
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20  dB.  The  sensitivity  was  still  below  the  shot-noise-limited  case  by  about  a  5  to  8  dB 
due  to  the  noise  hgure  of  the  amplifier  (primarily  driven  by  the  ASE  noise).  Often¬ 
times,  an  amplifier  is  used  at  the  transmitter  end  to  reduce  unwanted  amplification  of 
noise  sources  like  ASE  noise.  An  EDFA  was  used  in  this  research  due  to  its  low-noise 
performance  and  bandwidth  capability  in  the  Terahertz  [28]. 

Receiver  aperture  size  ( aperture  averaging).  As  the  receiver  optic  size  increases, 
more  power  is  collected  at  the  receiver.  In  addition,  the  total  received  power  fluctua¬ 
tion  also  decreases,  especially  when  the  size  of  the  receiver  increases  beyond  the  cor¬ 
relation  width.  The  correlation  width  is  pcw  =1  to  3[LA/(27t)]1/2  for  weak  turbulence. 
For  apertures  larger  than  pcw,  different  parts  of  the  beam  experience  uncorrelated 
variations  in  intensity  over  the  aperture.  As  a  result,  the  aperture  essentially  aver¬ 
ages  these  variations  by  focusing  the  intensity  pattern  onto  the  detector,  commonly 
referred  to  as  aperture  averaging.  Many  have  studied  this  phenomenon  in  regard  to 
FSOC  performance  [2,15,24,63,91-94].  This  research  uses  the  principle  of  reciprocity 
to  determine  the  optimal  separation  distance  for  multiple  transmitters  using  pcw.  This 
approach  is  first  discussed  in  Section  4.1.5. 

To  collect  as  much  power  as  possible  and  reduce  intensity  variations,  the  receiver 
optic  should  always  be  made  as  large  as  practical  for  any  given  system.  In  this  research 
a  single  nominal  receiver  size  is  used  to  evaluate  multiple-transmitter  system  designs. 

3.2  Wavefront  control 

Wavefront-control  systems  measure  and  correct  for  fades  in  real  time.  Typically, 
wavefront  sensors  measure  the  real-time  aberrations  of  the  propagation  path,  and  a 
closed-loop  AO  system  applies  a  correction  to  pre-compensate  the  beam  in  real  time. 
Most  simple  wavefront-control  systems  require  the  following:  a  beacon  on  the  receiver 
station,  a  sensor,  and  a  feedback  loop  to  adjust /control  the  transmitted  signal  in  real 
time  to  improve  the  link.  The  simplest  of  these  systems  involves  only  a  tracking 
system  with  a  fast-steering  mirror  to  keep  the  beam  centered  on  the  receiver.  Beam 
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walk-off  results  in  the  deepest  fades,  and  therefore  a  system  to  keep  the  beam  centered 
on  the  receiver  can  significantly  improve  performance. 

Tracking  and  Pointing.  Tracking  refers  to  the  real-time  estimation  of  the  trans¬ 
mitter  direction  needed  to  center  the  beam  on  the  receiver.  Pointing  refers  to  the 
actual  transmitter  direction  which  includes  the  jitter  of  the  transmitter  platform  and 
the  tilt  compensation  applied  to  a  fast-steering  mirror  to  account  for  the  atmospheric 
turbulence.  For  long  transmission  lengths,  seemingly  minor  platform  vibrations  and 
errors  in  the  control  system  at  the  transmitter  can  have  a  significant  effect  on  the 
direction  of  the  transmitted  beam.  This  overall  error  in  the  tracking-control  system 
is  referred  to  as  jitter.  For  a  100  km  path  an  uncontrolled  vibration  on  the  order  of 
10  /rrad  at  the  transmitter  translates  into  a  1  meter  beam  movement  at  the  receiver. 
Jet  Propulsion  Laboratory  and  the  California  Institute  of  Technology  presented  their 
approach  for  a  tracking  and  pointing  system  for  a  FSOC  system  for  the  International 
Space  Station  (ISS)  in  2000  [44],  Their  goal  was  a  downlink  capability  of  rates  up 
to  2.5Gbit/s  using  adaptive  optics,  tracking,  and  pointing  systems  to  compensate  the 
beam  for  the  turbulence  effects  [8,32,33,44],  Many  others  have  studied  and  designed 
various  tracking  systems  for  FSOC  systems  [35,62,79]. 

Platform  jitter  control.  In  Arnolds  research  on  receiver  jitter,  he  assumed  the 
beam-tracking  sensor  and  the  data  receiver  used  the  same  detector  [6].  The  satellite 
vibrations  at  the  Tx  caused  the  beam  to  non-uniformly  spread  across  the  tracking/- 
data  detector  array.  Since  the  data  receiver  used  the  same  detector  as  the  tracker, 
he  was  able  to  vary  the  gain  of  the  four  detectors  in  his  quad  cell  sensor  depend¬ 
ing  upon  the  received  tracker  signal,  improving  communication  system  performance. 
As  mentioned  earlier,  jitter  at  the  transmitter  can  be  very  severe.  For  instance,  the 
primary  cause  of  beam  pointing  error  for  the  proposed  ISS  downlink  is  platform  jit¬ 
ter,  not  the  atmosphere  [44].  This  is  due  to  the  long  propagation  length  and  the 
turbulence  layer  near  the  receiver,  not  the  transmitter.  Toyoshima  determined  the 
optimal  ratio  of  the  beam  divergence  to  the  angular  platform  jitter  to  reach  the  de¬ 
sired  BER  [75].  Sophisticated  tracking  systems  have  been  used  in  previous  research, 
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but  for  the  multiple-transmitter  case,  a  simple  centroid  tracker  is  used  as  a  baseline 
from  which  to  compare  future  work.  Many  more  sophisticated  tracking  systems  are 
used  to  help  counter  scintillation  effects.  Here,  the  multiple  Txs  help  to  alleviate 
those  effects.  In  addition,  this  research  implements  a  Rx  tracker  which  is  especially 
important,  since  the  Rx  aperture  is  smaller  than  coherence  diameter  r0  (defined  in 
Eq.  (30))  and  encounters  significant  AO  A  variations. 

AO  systems.  Higher-order  phase  aberrations  cause  the  beam  to  break  up  and 
spread  out  and  can  be  measured  with  a  wavefront  sensor  and  corrected  by  a  deformable 
mirror  [66,69,80].  Some  researchers  have  proposed  complicated  wavefront  control 
systems  to  sense  and  correct  the  transmitted  signals  in  real  time  to  reduce  fading 
at  the  receiver  [7,33,77,78].  Tyson  et  al.  reduced  the  BER  by  a  factor  of  41  in  a 
hardware-in-the-loop  experiment  using  AO  wavefront  control  [81].  These  bulky  AO 
systems  were  used  in  high  phase  perturbation  cases  of  D/tq  =  2  to  9  often  found  in 
near  field  turbulence  regimes.  These  phase  correction  systems  would  not  be  effective 
for  the  air-to-air  case  studied  here,  where  D/tq  =  1/2. 

Receiver  adaptive  optic  systems.  Wavefront-sensorless  systems  reduce  AO  sys¬ 
tem  complexity  by  removing  the  requirement  for  a  wavefront  sensor.  These  AO  sys¬ 
tems  attempt  to  maximize  the  received  power  by  optimally  adjusting  an  adaptive  el¬ 
ement.  Booth  proposed  an  efficient  method  to  measure  the  wavefront  Zernike  modes 
and  optimally  control  the  receiver  AO  [9].  Others  have  proposed  other  methods  and 
techniques  [53,54,87].  With  these  systems  there  are  no  complicated  adaptive  elements 
on  the  transmitter,  significantly  reducing  transmitter  size  and  weight.  However,  both 
types  of  higher  order  AO  systems  only  correct  for  the  phase  cannot  correct  the  strong 
scintillation  in  the  air-to-air  scenario.  The  multiple-Tx  approach  used  here  averages 
out  the  strong  scintillation  effects. 
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3.3  Diversity  techniques 

To  avoid  these  bulky  control  systems,  diversity  techniques  (i.e.  multiple  Txs  and 
Rxs,  and  interleaving)  have  been  used  in  the  literature  to  average  out  the  spatial  and 
temporal  variations  at  the  receiver.  Each  path  should  be  separated  from  the  others  to 
maximize  the  differential  variance  between  them.  Much  of  the  research  has  focussed 
on  multiple  Rxs.  Some  research  has  been  accomplished  on  multiple-Tx  systems,  but 
surprisingly  little  research  can  be  found  on  the  theoretical  angular  separation  required 
to  take  full  advantage  of  these  spatial  diversity  techniques,  especially  for  the  airborne 
application  [34,58,59,68].  Rather,  most  literature  on  differential  statistics  for  angular 
separations  is  concerned  with  an  isoplanatic  angle  (relatively  similar  paths)  for  sensing 
and  wavefront  control. 

Diversity  techniques  take  advantage  of  the  atmosphere’s  randomness  by  propa¬ 
gating  through  different  atmospheric  conditions  and  averaging  the  result.  Multiple- 
transmitter  systems  take  advantage  of  uncorrelated  spatial  atmospheric  effects  to 
average  out  the  phase  and  amplitude  variations  at  the  receiver.  This  diversity  can  be 
accomplished  in  either  time  or  space.  These  techniques  include  multiple  transmitters, 
multiple  receivers,  and  time  interleaving.  In  each  case,  the  goal  is  to  separate  the 
signals  far  enough  apart  that  the  fades  are  reduced  by  an  averaging  effect. 

Multiple-input  and  multiple-output  (MIMO)  FSOC  systems  have  been  studied 
for  both  coherent  and  direct-detection  (incoherent)  systems.  For  example,  spatial 
diversity  could  be  accomplished  at  the  receiver  by  a  single  large  aperture  (larger  than 
pc  or  multiple  small  apertures  [3] .  Some  researchers  have  studied  multiple-transmitter 
approaches  [34,58,59,68]  for  the  space-to-ground  and  the  ground-to-ground  scenarios, 
but  few  have  studied  the  air-to-air  scenario.  Much  of  the  spatial  diversity  research  has 
focused  on  MIMO  techniques  or  just  multiple  receivers  to  reduce  signal  fades,  with¬ 
out  defining  the  requisite  separation  distance.  In  this  research,  the  spatial  diversity 
techniques  are  optimized  for  the  particular  geometry  and  turbulence  profile. 
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There  are  some  advantages  to  multiple  small  apertures  on  the  airborne  platform 
over  a  single  optic  since  they  could  be  incorporated  into  conformal  optics.  Conformal 
optics  lend  themselves  well  to  airborne  platforms  for  size,  platform  aerodynamics,  as 
well  as  reductions  in  aero-optic  effects  around  the  optic  [52], 

Fortunately,  as  this  document  shows  later,  multiple-transmitter  techniques  on 
airborne  platforms  are  feasible  due  to  the  relatively  small  separation  distance  required 
for  uncorrelated  scintillation  effects  between  the  beams  (i.e.  31  cm  for  the  100  km 
path).  Multiple  transmitters  average  out  the  effects  of  scintillation  rather  than  using 
bulky  or  complicated  AO  systems  that  only  correct  for  the  phase.  The  scintillation 
effects  of  the  extended  turbulence  for  the  air-to-air  path  are  very  strong  (i.e.  scintil¬ 
lation  effects  weighted  heavily  in  the  center  of  the  path  as  shown  in  (48)  and  (49)), 
making  this  scenario  particularly  suited  for  a  multiple-Tx  system. 

3-4  Modulation  techniques 

Many  different  modulation  schemes  can  be  used  to  encode  the  information  on 
an  optical  signal.  There  are  coherent  techniques,  direct-detection  techniques  (incoher¬ 
ent),  polarization-modulation  schemes,  and  quantum  crypto-keying  techniques.  Some 
benefits  and  drawbacks  of  these  techniques  are  discussed  below. 

Coherent  techniques  require  extremely  precise  timing  since  optical  signals  are 
at  extremely  high  frequencies  (for  380-1550  nm  wavelengths,  the  frequencies  are  200- 
800  Terahertz).  I11  2006,  Lange  et  al.  used  coherent  modulation,  namely  homodyne 
binary  phase-shift  keying  (BPSK),  in  a  142  km  link  between  two  islands  with  a  bea¬ 
con  and  tracking  system.  They  transmitted  5.625  Gbit/s  with  a  BER  that  varied 
between  10~4  to  10”6,  showing  the  potential  robustness  of  coherent  techniques  [42], 
They  cited  its  immunity  to  sunlight  interference  as  a  decisive  factor  in  its  potential 
for  commercial  use.  There  are  coherent  techniques  that  can  take  advantage  of  the 
diversity  of  multiple-transmitter  techniques.  For  example,  Haas  showed  a  marked  im- 
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provement  for  a  coherent  heterodyne  modulation  scheme  by  using  multiple  space-time 
coded  transmitters  [29]. 

Direct-detection  systems  modulate  the  transmitted  power  incoherently  and  do 
not  require  the  precise  timing  of  coherent  techniques.  One  of  the  simplest  of  these 
incoherent  techniques  is  on/off  keying  (OOK)  [74],  Many  have  used  this  binary  digital 
communication  technique  in  FSOC  systems  in  which  a  ‘1’  is  represented  by  a  pulse 
and  the  ‘0’  is  represented  by  the  absence  of  a  pulse  [1,39,95].  Another  direct-detection 
technique  modulates  the  polarization  state  of  the  transmitted  optical  signal.  One  such 
technique  is  polarization  shift  keying  (PolSK)  [76].  Since  turbulence  is  isotropic  and 
does  not  appreciably  affect  the  polarization  state,  it  seems  to  be  an  ideal  modulation 
scheme  at  first.  In  this  case  when  a  fade  occurs,  it  affects  all  symbols  equally.  Unfor¬ 
tunately,  the  key  problem  is  during  a  signal  fade  the  signal  might  not  be  high  enough 
to  adequately  determine  which  polarization  was  sent.  Often,  the  signal  level  required 
to  determine  the  polarization  is  much  higher  than  the  level  needed  to  determine  if 
the  laser  is  on  or  off  as  in  OOK. 

Another  laser  communication  technique  is  quantum  cryptography  with  entan¬ 
gled  photons,  where  its  promise  lies  in  its  security  aspects.  The  polarization  state  of 
polarization-entangled-photon  pairs  is  modulated  and  received.  The  Heisenberg  un¬ 
certainty  principle  guarantees  that  an  eavesdropper  disturbs  the  signal  and  therefore 
is  detected  [38,61,73].  Current  systems  can  only  propagate  at  very  low  bit  rates. 

For  simplicity,  ease  of  use  with  multiple-transmitter  systems,  sensitivity  at  very 
low  power  levels,  and  to  focus  on  mitigation  of  turbulence  effects,  OOK  is  used  ex¬ 
clusively  in  this  research.  After  implementing  these  techniques,  there  are  still  signal 
fades  that  need  to  be  mitigated  in  some  way.  Signal  processing  techniques  mentioned 
in  the  next  section  can  allay  the  effects  of  these  fades  with  little  to  no  added  size  or 
weight. 
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3.5  Signal  processing  techniques 

Once  the  signal  has  been  received,  signal-processing  techniques  like  adaptive 
thresholding,  automatic  gain  control,  interleaving,  and  forward  error-correction  coding 
can  be  used  to  decipher  symbols  more  accurately,  resulting  in  improved  BER.  Most 
of  these  techniques  have  been  applied  to  a  single-transmitter  system,  while  others 
have  been  applied  to  a  coherent  phase  multiple-transmitter  system.  Still  others  have 
been  proposed  for  the  short,  constant  turbulence  strength  of  the  “last-mile”  horizontal 
propagation.  These  techniques  are  driven  by  the  temporal  statistics  and  conditional 
probabilities  of  the  FSOC  system. 

The  BER  of  a  digital  communication  system  depends  upon  the  conditional  prob¬ 
ability  distributions  of  the  received  signals  [82],  The  communication  system  can  be 
improved  by  determining  the  optimum  decision  threshold.  Since  these  probabilities 
vary  widely  (but  slowly  in  time  compared  with  the  signal  frequency)  due  to  the  at¬ 
mospheric  conditions,  adaptive  detection  thresholds  can  reduce  the  BER  of  FSOC 
systems  [12-14,18].  Researchers  have  successfully  used  various  estimation  techniques 
to  adaptively  determine  the  optimum  threshold  to  include  least-mean-square  predic¬ 
tors  [13],  Kalman  filters  [12-14],  and  maximum  likelihood  estimators.  The  temporal 
statistics  of  the  received  signal  dictate  what  type  of  estimator  is  adequate.  Shown 
later  a  low  sampling  rate  linear  maximum  likelihood  estimator  was  sufficient,  in  this 
research,  since  the  signal  variation  was  reduced  by  the  multiple-transmitter  system. 
A  lower  sampling  rate  leads  to  a  higher  SNR,  allowing  for  less  signal  diverted  to  the 
estimator. 

Many  investigators  have  proposed  forward  error  correction  codes  (FECs)  and 
interleavers  for  FSOC  systems.  Interleavers  rearrange  the  symbols  in  time  so  that 
the  errors  caused  by  a  long,  deep  fade  are  spread  out  so  the  decoder  can  handle  the 
errors  as  if  they  were  random  errors  [74].  Zhu  and  Kahn  determined  an  upper  bound 
for  the  code  gain  of  a  weak-turbulence  optical  channel.  They  considered  block,  con¬ 
volutional,  and  Turbo  codes  in  conjunction  with  varying  the  interleaver  length  and 
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determined  that  BER  performance  continues  to  improve  with  increasing  interleaver 
length  [96] .  However,  practical  interleaver  length  is  limited  due  to  system  complexity 
and  delays  in  coding  and  decoding.  Turbo  codes  were  shown  to  out-perform  convolu¬ 
tional  codes  for  low  signal  to  noise  FSOC  systems  [55].  Yu  et  al.  demonstrated  that 
Reed-Solomon  codes  could  improve  performance  for  log-normal  atmospheric  intensity 
variation  statistics  [90].  Since  interleaver  length  is  determined  by  the  duration  or 
span  of  the  channel  memory,  not  the  particular  statistical  nature  of  the  effects,  fade 
statistics  can  be  used  to  determine  interleaver  length.  In  Section  5.4.1  the  multiple-Tx 
approach  is  shown  to  significantly  reduce  fade  lengths,  thereby  reducing  the  required 
interleaver  length.  Interleavers  and  FEC  systems  are  beyond  the  scope  of  this  re¬ 
search,  but  they  could  be  implemented  in  conjunction  with  the  techniques  presented 
in  this  research. 

Automated  gain  control  (AGC)  systems  are  used  to  optimize  particular  detector 
and  digital  receiver  performance.  AGC  systems  can  be  either  optical  or  electrical. 
Power  levels  at  the  receiver  vary  greatly,  especially  for  the  long  extended  turbulence 
case  of  air-to-air  laser  communication.  An  optical  AGC  can  reduce  this  fluctuation 
so  that  the  optical  detector  can  be  tailored  to  a  particular  optical  signal  range  [10, 
37].  Since  this  research  does  not  prescribe  a  particular  detector,  an  AGC  was  not 
implemented.  In  addition,  the  multiple  Tx  system  reduces  the  variability  and  the 
adaptive  threshold  system  adjusts  the  threshold  with  time.  Therefore,  likely  reducing 
the  effectiveness  and  need  for  an  AGC  system.  The  EDFA  model  used  in  this  research 
could  be  refined  to  include  the  optical  gain  saturation  inherent  in  EDFAs,  providing 
some  AGC  to  further  improve  the  bottom-line  performance. 

3.6  Assessment  of  best  areas  for  further  research 

This  research  strives  to  reduce  the  hardware  complexity,  while  maintaining  the 
performance  required  for  an  airborne  laser  communication  system.  This  requires 
optimization  of  each  subsystem.  The  goal  is  to  achieve  the  highest  data  rate  at  the 
lowest  BER,  while  maintaining  a  level  of  simplicity  (i.e.  size,  weight,  and  power) 
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suitable  for  an  airborne  platform.  Specifically,  the  goal  of  this  research  is  to  show 
how  a  multiple-beam  system  coupled  with  signal  processing  techniques  can  attain 
reliable  communication  at  high-bit  rates  without  requiring  complicated  higher-order 
control  systems. 

The  techniques  used  in  this  research  work  well  together,  since  the  physical  layer 
improvements  reduce  the  variability  of  the  received  power  and  the  signal  processing 
approaches  decipher  the  bits  more  effectively.  The  Rx  and  Tx  tracking  systems  correct 
for  the  low-order  phase  (AOA  variance  and  beam  wander),  the  multiple  Txs  correct 
for  the  high-order  phase  and  amplitude  effects,  and  the  adaptive  threshold  counters 
the  temporal  variations  of  the  signal  level  regardless  of  the  turbulence  severity.  As 
is  shown  in  Section  5.4.2,  each  of  these  techniques  contributes  to  improve  the  per¬ 
formance,  pushing  the  BER  for  the  high-scintillation  case  below  much  more  benign 
turbulence  cases.  All  of  these  approaches  improve  the  performance  at  the  physical 
layer,  so  that  interleavers  and  FEC  systems  could  be  added  to  reduce  the  BER  even 
further. 


58 


IV.  Anisoplanatic  Turbulence  Effects 


The  same  atmospheric  turbulence  effects  that  limit  the  resolution  of  optical  sys¬ 
tems  and  make  the  stars  twinkle  can  severely  reduce  the  amount  of  power 
received  in  an  FSOC  system.  The  atmospheric  turbulence  in  the  propagation  path 
causes  the  laser  beam  to  wander,  spread,  and  break  up.  These  effects  can  cause  the 
received  signal  power  to  drop  below  the  receiver’s  threshold  for  milliseconds  at  a  time. 
For  an  FSOC  system,  a  millisecond  fade  means  millions  of  bit  errors.  Since  these  op¬ 
tical  power  fades  are  often  very  deep,  simply  turning  up  the  power  in  this  case  would 
not  be  very  effective. 

Two  different  ways  to  improve  this  condition  are  to  increase  the  signal  diversity 
to  average  out  the  effects  or  compensate  for  the  turbulence  conditions  in  real  time. 
In  the  first  approach,  the  temporal  and  spatial  statistics  of  the  turbulence  for  the 
propagation  are  estimated.  Then,  techniques  are  devised  to  overcome  these  effects 
by  applying  multiple  uncorrelated  realizations.  In  the  second  approach,  typically, 
wavefront  sensors  measure  the  real-time  phase  aberrations,  and  a  closed-loop  AO 
system  applies  a  correction  to  pre-compensate  the  transmitted  beam’s  phase  in  real 
time. 

Multiple-transmitter  systems  increase  signal  diversity  and  average  out  the  dele¬ 
terious  effects  of  turbulence  without  bulky,  complicated  AO  systems,  making  it  an 
appropriate  choice  for  airborne  laser  communications.  In  addition,  AO  systems  do  not 
correct  for  scintillation.  Through  analysis  and  simulation,  the  optimal  configurations 
for  multiple-Tx  airborne  FSOC  systems  are  determined  for  various  geometries  and 
tracking  systems. 

This  research  derives  the  requisite  angular  and  parallel  separations  for  multiple- 
Tx  systems  for  airborne  and  ground-to-ground  laser  communication.  A  majority  of 
the  previous  research  on  multiple  transmitters  has  focused  on  satellite  communications 
(in  which  the  turbulence  is  only  present  over  a  short  part  of  the  propagation  path) 
or  ground-based  constant-turbulence-strength  paths  [4,5,30,40,59,60,68].  Here,  this 
work  presents  these  angular  separations  for  three  airborne  geometries  (air-to-air,  air- 
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to-ground,  and  ground-to-air)  through  extended  turbulence  and  determines  practical 
configurations.  Chapter  III  provided  more  details  of  previous  research  and  other 
approaches. 

Previous  research  on  isoplanatism  has  defined  the  maximum  angle  over  which 
the  turbulence  effects  between  two  paths  is  relatively  similar  [26,69,70].  These  isopla- 
natic  angles  have  been  determined  for  the  tilt  variance,  higher-order  phase  variance, 
and  scintillation  (intensity  variance).  This  research  goes  further  to  determine  the 
minimum  angle  at  which  the  paths  are  relatively  different  (i.e.  anisoplanatic).  The 
less  correlated  the  amplitude  and  phase  perturbations  between  the  paths,  the  better 
the  averaging  effect  for  multiple  beams. 

With  these  separations  computed,  wave-optics  simulations  were  conducted  to 
explore  how  separation  distances  affect  the  BER  for  multiple-Tx  FSOC  systems.  The 
simulations  are  performed  for  multiple  scenarios  and  tracking  systems  to  determine 
how  effective  these  multiple-Tx  techniques  might  be  for  airborne  platforms. 

4-1  Uncorrelated  paths 

This  section  determines  the  separation  required  to  attain  uncorrelated  turbu¬ 
lence  effects  between  two  laser  beam  paths.  To  investigate  this,  it  is  instructive  to  first 
determine  when  the  paths  are  relatively  similar.  If  a  system  effect  is  space-invariant, 
it  is  called  isoplanatic  [27].  Therefore,  if  two  laser  beam  paths  are  considered  iso- 
planatic  in  terms  of  any  particular  turbulence  effect,  the  effects  of  the  two  paths  are 
highly  correlated. 

Most  AO  systems  have  a  beacon  path  to  measure  the  turbulence.  Sensors  at  the 
imaging  system  or  laser  transmitter  measure  how  the  turbulence  affects  the  beacon. 
If  the  differences  between  the  phase  effects  (wavefront  variations)  of  the  propagation 
path  and  beacon  path  are  negligible,  the  phase  correction  can  potentially  be  imple¬ 
mented  effectively.  That  is  to  say,  the  phase  effects  of  the  paths  are  isoplanatic.  The 
phase  isoplanatic  angle  9q  is  the  largest  angle  between  two  paths  for  which  the  wave- 
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front  variations  in  the  two  paths  are  relatively  similar  [66].  The  relation  for  this  angle 
is  presented  in  Eq.  (31)  in  Chapter  II.  If  the  paths’  effects  are  significantly  different, 
the  paths  are  anisoplanatic. 

There  are  three  different  types  of  isoplanatism  of  interest  in  this  research:  the 
tilt  6tai  phase  9q,  and  scintillation  isoplanatic  9X0  angles.  Tilt  refers  to  the  direction 
of  propagation  and  deals  with  tracking  a  wandering  beam  or  a  jittering  image.  Phase 
incorporates  both  the  tilt  and  the  higher-order  phase  aberrations.  Scintillation  cor¬ 
responds  to  the  variations  in  intensity  over  the  pupil.  Typically,  the  tilt  isoplanatic 
angle  is  larger  than  the  phase  isoplanatic  angle,  which  is  larger  than  the  scintillation 
isoplanatic  angle.  Using  these  isoplanatic  conditions  as  a  starting  point,  the  aniso¬ 
planatic  conditions  are  determined  for  the  phase  and  amplitude  effects.  The  phase 
anisoplanatic  angle  9^ind  was  first  derived  in  support  of  this  research  [46,47,49]. 

The  refractive-index  fluctuations  drive  the  phase  and  amplitude  turbulence  ef¬ 
fects.  For  the  derivations  in  this  section,  the  von  Karman  refractive  index  fluctuation 
PSD  models  these  fluctuations  [2, 66] 

0.033Cfc) 

J  (ft2  +  K§)n/6’ 

where  k  is  the  3-D  radial  spatial  frequency  and  Ko  =  2tt/Lq.  This  PSD  is  equivalent 
to  Eq.  (16),  with  the  inner  scale  l0  set  to  zero.  This  PSD  is  the  most  appropriate 
since  it  includes  the  outer  scale  L0  which  limits  the  size  of  the  large-scale  phase  effects 
(i.e.  turbulent  eddies)  which  drive  the  phase  anisoplanatic  conditions.  The  refractive 
index  PSD  used  to  derive  the  isoplanatic  conditions  consists  of  only  the  numerator  in 
Eq.  (131),  since  the  outer-scale  does  not  affect  the  isoplanatic  conditions. 
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Figure  12:  Isoplanatic  angle  geometry. 


4-1.1  Phase  isoplanatism.  Fried  derived  the  phase  perturbation  structure 
function  in  1966  as  [23,66] 


A/,  (Ax) 


(hHx)  -  V’(x  +  Ax)]2) 

2B^ (0)  —  2B^ (Ax),  for  stationary  random  processes, 

r-L 

2.91k2  (Ax)5/3  /  C2(z)dz, 


(132) 

(133) 

(134) 


where  (•)  is  the  expectation  operator  and  B^  is  the  auto-correlation.  The  ^  term 
denotes  the  pupil  phase  perturbation  and  Cl(z)  is  the  strength  of  turbulence  along 
the  path.  It  can  be  shown  from  Eq.  (134)  that  the  phase  structure  function  at  the 
receiver  for  two  point  sources  separated  by  angle  6  as  viewed  by  the  receiver  is 


D46,L) 


2.91  k2  [sin  (0)]5/3 


(L  -  z)5/3  C2n  (z)  dz. 


(135) 


This  angular  separation  is  shown  in  Fig.  12.  For  this  geometry,  most  of  the  literature 
has  defined  the  isoplanatic  angle  to  be  the  angle  at  which  the  structure  function  is 
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less  than  or  equal  to  unity  [2,26,66].  Applying  this  condition, 


L  j  =  1  rad2, 


(136) 


and  solving  for  the  angle  results  in  the  familiar  isoplanatic  angle  relation  defined 


by  [26] 


L 


0O=  2.91k2  /  C2(z)(L  -  z)s/3  dz  .  (137) 


(137) 


It  is  important  to  note  which  isoplanatic-angle  definition  is  used.  This  definition 
assumes  the  two  point  sources  are  in  the  z  =  0  plane.  Many  definitions  in  the 
literature  define  the  z  —  0  point  as  the  receiver  location  [2,66,69].  To  adjust,  let 
z'  =  L  —  z  in  Eq.  (137). 

4-1.2  Angular  phase  independence  of  two  beams.  Now,  these  concepts  are 
applied  to  the  statistical  independence  to  determine  the  phase  independence  angle. 
The  phase  structure  function  in  Eq.  (135)  increases  with  separation  angle,  approaching 
a  maximum  value  at  two  times  the  mean  square  phase  or  2cr^,  as  the  two  paths  are 
placed  far  apart.  This  time,  setting  the  structure  function  equal  to  the  maximum 
value  for  the  phase  structure  function  yields  [46] 


Dll>{QlpindlL)  — 


(138) 


Combining  Eqs.  (136)  and  (138)  leads  to  a  way  to  solve  for  the  phase  independence 
angle  yielding  [46] 
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Using  a  geometrical  optics  plane-wave  propagation  approximation,  the  phase  variance 
for  a  point  receiver  can  be  written  as  [2] 

pL  poo 

a^  pl  =  4ir2k2  /  /  z)<\k,  dz  (140) 

Jo  Jo 

rL 

=  0.78/c2Kq  5/3  /  C2(z)  dz.  (141) 

Jo 

For  horizontal  propagation  (i.e.  constant  C2),  the  independence  angle  simplifies  to 
[46,47,49] 

(142) 

As  first  derived  in  support  of  this  research  [46,47,49],  this  relation  for  9^ind  defines  the 
angle  over  which  the  phase  effects  between  the  propagation  paths  of  two  point  sources 
are  nearly  uncorrelated.  It  follows  that  the  phase- independence  separation  distance 
can  be  defined  as  d</,ind  =  L9^irid.  As  expected,  the  phase-independence  angle  increases 
with  outer  scale.  At  this  angular  separation,  the  beams  should  wander  independently, 
and  the  higher-order  phase  perturbations  should  be  uncorrelated  as  well.  At  this 
separation  a  fixed  multiple-transmitter  LCS  (e.g.  last-milc-type  communications) 
could  be  designed  so  that  at  least  one  beam  with  sufficient  power  remains  on  the 
receiver  at  all  times  without  the  need  for  tracking.  This  independence  angle  is  highly 
dependent  on  the  outer  scale,  which  varies  near  the  ground  as  L0  cs  OAh  [85].  For 
example,  two  A  =  1.55  /irn  transmitters  would  need  to  be  separated  by  dyind  =  43  cm 
{9^ind  =  213  /rrad)  for  a  2  km  path  located  1  m  above  the  ground  with  a  turbulence 
strength  of  C2  =  1.71  x  10”14  m-2//3.  For  a  4  km  path,  that  separation  would  need 
to  approach  65  cm.  Figure  13  shows  the  requisite  transmitter  separations  d^ind(/r)  for 
four  different  path  lengths  given  a  particular  height  above  the  ground  for  each  path. 

At  high  altitudes,  the  effective  outer  scale  is  determined  by  the  vertical  outer 
scale  and  the  horizontal  outer  scale.  The  vertical  outer  scale  typically  varies  from  10  to 
70  nr  [21],  while  the  horizontal  outer  scale  can  be  much  larger.  Aircraft  measurements 
have  determined  the  horizontal  outer  scale  can  be  over  hundreds  of  kilometers  [85]. 
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Figure  13:  The  phase  independence  distance  d^ind  is  plotted  versus  the  height  above 
the  ground.  Horizontal  path  with  turbulence  strength  for  each  path  corre¬ 
sponding  to  the  HV57  profile,  L  —  1,  2,  4,  6  km,  and  A  =  1.55  /ini  [46]. 

For  horizontal  propagation  simulations  in  this  work  at  altitude,  an  infinite  outer  scale 
is  used  because  L0  D.  When  a  finite  L0  is  needed  with  slant  ranges,  the  effective 
outer  scale  is  determined  by  taking  a  slice  through  the  vertical  outer  scale 

L0  =  (143) 

cos  4 

where  Lovert  is  the  outer  scale  for  vertical  propagations  and  4  is  the  zenith  angle.  Both, 
the  outer  scale  and  inner  scale  Iq  vary  with  altitude.  In  this  research,  these  bounds 
on  the  turbulence  are  consistent  with  atmospheric  data  presented  by  Wheelon  [85]. 
See  Section  2.2.3  in  Chapter  2  for  more  details. 


4-1.3  Parallel  path  isoplanatism.  Using  Eq.  (134)  again,  this  work  deter¬ 
mines  the  parallel  path  isoplanatic  distance.  Referring  to  the  phase  structure  function 
at  the  receiver  for  two  parallel  path  point  sources  separated  by  Ax  yields  [46, 66] 


A/,(  Ax,  L)  =  2.91  k2  (Axf3  /  C2n(z)  dz. 


(144) 
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As  for  the  angular  isoplanatic  angle  condition,  one  can  determine  the  separation  Axq 
at  which  the  structure  function  is  unity.  The  parallel  isoplanatic  distance  is  [46] 


Axn  = 


2.91  k2  [  C2(z)  dz 

Jo 


-3/5 


=  0.6611p0  =  0.3148r0, 


(145) 


where  p0  is  the  spatial  coherence  radius.  [2, 66]  For  a  constant- C'^  path,  [46] 


Ax0  =  0.5268  k-d/5C-e/5L-3/b  =  (3/8  )3/5L60. 


(146) 


Interestingly  enough,  this  separation  is  simply  (3/8)3//5  times  the  separation  for  an¬ 
gularly  separated  paths.  Similar  to  Eq.  (138),  the  plane-wave  independent  phase 
separation  distance  A Xind  is  determined  by  setting  the  structure  function  equal  to  the 
maximum  value  and  solving  for  the  separation  [46] 


Axind  =  2alplAx0  =  0.4109ac0  5/3 


L  -l  2/5 

Cl{z)  dz 


(147) 


For  a  constant-C'Z  prohle,  it  simplifies  to  [46] 


Axind  =  0Amk4/5C*/5L2/5K~5/3  =  (3/8  )3/5M 


'tftind  * 


(148) 


4- 1.4  Tilt  isoplanatism.  Sasiela  developed  relationships  for  the  differen¬ 
tial  tilt  variance  which  can  also  be  referred  to  as  the  structure  function  a'^(x)  = 
E[T(xi)  —  T(x  1  +a;)]2  of  the  Zcrnike  tilt  T.  This  work  uses  Sasiela’s  notation  to 
allow  the  reader  to  follow  this  work  and  refer  back  to  Sasiela’s  [69].  From  those  rela¬ 
tions,  he  determined  a  relation  for  the  tilt  isoplanatic  angle  for  an  astronomical-seeing 
geometry.  To  simplify  the  relations  for  tilt,  the  Kolmogorov  refractive  index  PSD  is 
used  since  it  does  not  include  the  outer  scale.  Sasiela  investigated  the  effect  of  outer 
scale  on  the  tilt  isoplanatic  angle.  The  outer  scale  greatly  affects  tilt  variance,  but 
does  not  appreciably  affect  tilt  isoplanatism  (especially  when  the  outer  scale  is  much 
larger  than  the  receiver  aperture)  [70].  The  differential  tilt  (i.e.  the  difference  between 
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the  Z-tilts)  consisted  of  two  contributions:  one  for  Arc  <  D  (beams  overlap)  called  the 
lower  contribution  and  the  other  for  Ax  >  D  (beams  do  not  overlap)  called  the  upper 
contribution,  where  D  is  the  receiver  diameter.  The  differential  tilt  variance  differs 
with  each  axis:  the  beam  displacement  axis  is  denoted  by  the  parallel  symbol  (||),  and 
the  perpendicular  symbol  (_L)  denotes  the  other  axis.  Even  though  this  difference  in 
tilt  variance  can  be  quite  significant  [49,51,69],  often  these  two  orientations  are  added 
to  determine  the  total  differential  tilt,  since  the  tilt  effects  in  the  x  and  y  directions 
are  independent  [69].  The  total  differential  tilt  is  the  sum  of  the  lower  and  upper 
contributions  for  each  axis  [69] 
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Sasiela  determines  the  lower  and  upper  contributions  by  calculating  the  differ¬ 
ential  tilt  variance  for  the  lower  and  upper  condition.  The  lower  contribution  for 
Ax  <  D  (when  the  beams  overlap)  is  [69] 
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where  Lc  =  D/d  and  Hc  =  Dcos(£/)/9.  The  function  pFj(a);  (6);  z\  is  the  generalized 
hypergeometric  function  where  a  and  b  are  p—  and  q— dimensional  arrays,  respectively. 
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The  upper  contribution  for  Ax  >  D  (beams  do  not  overlap)  is  [69] 
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where  L  is  propagation  length  and  H  is  the  transmitter  height.  Above  about  30  km 
the  turbulence  is  small  and  practically  negligible.  If  the  turbulence  strength  is  in 
terms  of  the  propagation  distance  z,  select  Lc  =  D /9,  and  if  the  relations  are  in  terms 
of  the  height  h,  use  Hc  =  Dcos(£,)/9. 

Sasiela  goes  on  to  derive  a  relation  for  the  tilt  isoplanatic  angle  QTa  using  the 
lower  contribution  for  an  astronomical  seeing  geometry.  For  small  displacement  angles 
where  Hc  is  higher  than  the  uppermost  turbulence,  the  tilt  isoplanatic  angle  is  defined 
in  terms  of  Eq.  (150)  using  the  following  differential  tilt  jitter  [69]: 


2.67 !A(Hc)  (eV 
D1/3  \D)  ' 

3.68rt  fe\\^  ,  2.35/i-/3  fey,,3f20\ 

— Dvr(cJ  (6)  +  -d^[d)  { j)  + 


(152) 


This  angle  is  characterized  as  the  angle  in  which  the  differential  tilt  jitter  is  equal  to 
one-half  the  diffraction-limited  beam  width.  For  the  definitions  of  the  moments  of 
the  propagation  path  /xTO,  /i“,  and  refer  to  Eqs.  (23)  -  (26). 

Using  the  first  factor  in  Eq.  (152)  for  6  <  .D/40,  000  and  setting  a^(9)  equal  to 
one  half  of  the  diffraction-width  angle  yields 


0.61A\ 2  _  2.67/12  (0ta\2,,, 
D  )  DV 3  V  D  J  U 


(153) 


Solving  for  9ta  as  Sasiela  did,  the  tilt  isoplanatic  angle  is  defined  as  the  angle  between 
two  sources  for  which  the  tilt  jitter  is  equal  to  one-half  the  diffraction  limited  spot 
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size.  Therefore,  Ota  is  [69] 


Ota 


0.184A.D1/6 

bb1/2 


(154) 


Using  results  published  by  Sasicla  for  the  differential  tilt  between  two  beams, 
i.e.  Eqs.  (150)  and  (151),  an  approximation  is  developed  here  for  the  differential  tilt 
for  a  constant-^  profile  or  horizontal  propagation  that  takes  into  account  both  the 
lower  (Ax  <  D)  and  upper  (Ax  >  D)  contributions.  Starting  with  Ax  <  D,  one 
must  only  consider  the  lower  contribution  for  6  <  D/L ,  where  0  =  A x/z.  After 
performing  the  integration  over  2  from  0  to  the  propagation  length  L,  the  differential 
tilt  becomes  [46] 
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where  Ota  <  D/L.  As  0  approaches  zero,  the  differential  variance  approaches  zero, 
as  expected.  To  define  the  tilt  isoplanatic  angle  as  Sasiela  did,  the  infinite  sum’s  first 
term  is  set  equal  to  one  half  of  the  diffraction-width  angle, 
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Solving  for  Ota  yields  the  tilt  isoplanatic  angle  for  horizontal  propagation  as  approx¬ 
imately  [49] 

nsunnve 

(157) 


eTA  =  °-^l,  0  <  D/L. 


CM'2  ’ 


This  straightforward  equation  for  constant  turbulence  strength  can  only  be  used  when 
the  beams  overlap  and  Ax  <  D.  For  the  astronomical  viewing  condition  the  simplifi¬ 
cation  to  only  use  the  lower  contribution  is  valid  since  turbulence  above  about  30  km 
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is  negligible.  Now  for  the  horizontal  propagation,  considering  9  >  D / L,  and  following 
Eq.  (149),  the  most  significant  terms  must  be  added  to  the  lower  and  upper  portions 
to  determine  the  overall  differential  tilt  variance 
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The  parallel  and  perpendicular  tilt  variances  are  added  to  determine  the  overall  tilt 
variance: 
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(159) 

Unlike  Eqs.  (153)  and  (156)  the  terms  in  Eq.  (159)  include  the  upper  and  lower 
significant  contributions  and  therefore  includes  additional  terms.  As  before,  one  could 
solve  for  9ta,  this  time  numerically,  to  determine  the  tilt  isoplanatic  angle  for  9  > 
D/L. 


A  1.5  Scintillation  anisoplanatism.  Stars  twinkle,  but  the  moon  and  even 
the  planets  do  not  twinkle  in  the  night  sky  because  their  angular  extents  are  much 
larger  than  the  scintillation  independence  angle.  In  weak  turbulence,  the  angle  at 
which  two  point  sources  scintillate  independently  was  postulated  by  Fried  to  be  9Xind  = 
0.8 (Lfc)-1/2,  [25]  corresponding  to  a  separation  distance  of  dXind  =  0.8  (L/k)1^2.  This 
relation  is  very  similar  to  the  correlation  width  pcw  defined  as  the  1/e2  point  of  the 
normalized  irradiance  covariance  function.  [2]  Since  pcw ,  for  weak  turbulence  varies 
between  1  to  3  Fresnel  zones  (L/k)1^2  depending  on  beam  size,  pc  is  referred  to  in  this 
work  as  simply  [2,46] 


pc  =  yJTJk. 


(160) 
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For  strong  turbulence  (lZSph  >  0.25)  the  scintillation  saturates  and  the  correlation 
width  of  irradiance  fluctuations  pcw  is  driven  by  the  spatial  coherence  radius  po  and 
the  scattering  disk  L/{kpf). 

The  correlation  width  pcw  is  often  used  to  describe  the  receiver  size  at  which 
aperture  averaging  occurs  as  the  receiver  size  increases.  Here,  the  principle  of  reci¬ 
procity  determines  the  angular  transmitter  separation,  referring  to  this  relation  as 
the  scintillation  correlation  angle  9Xc  =  (Lk)~1/'2.  The  analogous  angular  relation  to 
the  correlation  width  pcw  is  the  scintillation  independence  angle  9xind  defined  as  the 
angle  at  which  the  scintillation  between  the  two  paths  are  relatively  uncorrelated. 
The  values  of  9Xc  for  propagation  lengths  of  100  km  and  29  km  at  A  =  1.55  pm 
are  1.57  prad  and  2.91  prad,  respectively.  The  values  of  9xin(i  are  determined  by 
simulation  in  relation  to  scalar  multiples  of  9Xc  and  pc  in  Section  4.5. 


4-1-6  Considerations  of  isoplanatic  and  anisoplanatic  effects.  As  mentioned 
in  Section  4.1.2,  the  anisoplanatic  condition  can  be  determined  by  analyzing  the 
structure  functions  of  the  effects.  In  previous  work,  the  analytic  log- amplitude  Dx(d ) 
and  phase  structure  functions  D^(d)  were  determined  for  a  horizontal  path  as  [45,46] 
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where  Jo  is  the  zeroth-order  Bessel  function  of  the  first  kind  and  the  von  Karman 
PSD  from  Eq.  (131)  was  used  to  model  the  turbulence.  The  structure  functions 
are  plotted  with  the  corresponding  isoplanatic  and  anisoplanatic  distances  in  Fig.  14. 
Starting  with  the  phase  effects  plotted  against  the  right  scale  with  a  dash  dot  line,  the 
isoplanatic  angle  occurs  when  the  phase  structure  function  is  unity.  As  the  separation 
widens,  the  tilt  effects  are  isoplanatic  until  the  tilt  isoplanatic  angle  is  reached.  The 
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only  significant  difference  in  these  phases  is  due  to  the  higher-order  phase.  Finally 
at  separations  on  the  order  of  about  pcw  =  2 L0,  all  of  phase  effects  including  tilt  are 
anisoplanatic  between  the  two  paths.  The  amplitude  effects  are  plotted  against  the 
left  scale  with  a  solid  line.  At  about  2 pc,  the  structure  function  reaches  a  maximum 
and  settles  into  a  value  of  two  times  the  mean  square  log-amplitude  variance,  as 
the  amplitude  effects  become  uncorrelated.  These  separations  are  determined  by  the 
Fresnel  zone  ( L/k )1/2  and  are  consistent  with  results  for  weak  turbulence,  [46]  i.e. 
Rytov  number  lZsph  <  0.25  mentioned  in  Section  4.1.5  [2], 

Anguita  et  al.  simulated  much  stronger  turbulence  in  a  ground-to-ground  prop¬ 
agation  scenario  where  TZSph  =  1-6.  Their  uncorrelated  separation  distances  were 
greater  (approximately  6  pc )  due  to  the  long  correlation  tail  of  the  strong  turbulence 
characterized  by  the  scattering  disk  L  ( kp0 )  [2,5].  Polynkin  and  Peleg  also  simulated 
strong  turbulence  in  a  ground-to-ground  configuration  to  study  the  scintillation  re¬ 
duction  of  multiple  transmitters,  but  their  research  also  did  not  cover  the  air-to-air 
turbulence  effect  regime  [59,60]. 

The  relations  for  isoplanatic  and  anisoplanatic  effects  are  compared  for  different 
scenarios  in  Fig.  15.  Horizontal  propagation  near  the  ground  is  shown  in  Fig.  15b. 
If  the  transmitters  are  separated  by  the  phase-independence  angle  for  the  ground-to- 
ground  scenario,  then  tracking  might  not  be  required.  Provided  there  is  a  sufficient 
number  of  transmitters,  the  beams  would  wander  independently  with  at  least  one 
beam  on  the  receiver  at  any  given  moment.  The  fixed  pointing  angle  could  be  de¬ 
termined  by  maximizing  the  long-term  irradiance  for  each  beam.  The  isoplanatic 
angle  and  the  scintillation  correlation  angle  cross  at  about  2.5  km.  For  propagations 
beyond  the  cross-over  point,  scintillation  is  more  correlated  than  phase  effects.  In 
Fig  15c,  these  terms  cross,  too,  this  time  after  propagating  about  100  km.  This  also 
corresponds  very  well  with  Fig.  14  where  the  isoplanatic  angle  and  the  scintillation 
correlation  angle  nearly  coincide  with  90  slightly  smaller  than  6Xc  for  the  100  km 
air-to-air  scenario  [46]. 
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(p)  a 


Figure  14:  The  phase  (right  scale)  and  log-amplitude  (left  scale)  structure  functions 
are  plotted  for  a  100  km  horizontal  propagation  at  10  km  altitude,  with 
angularly  separated  beams.  The  strength  of  turbulence  is  L0/r0  =  286  and 
L^/ (XL)  =  23225  [46].  Isoplanatic  and  anisoplanatic  separation  distances 
are  shown  for  the  amplitude  and  phase  effects. 
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Table  1:  Required  beam  separation  for  averaging. 


(6  a  >  6^irid  ~  6Tind  >  6TA  >  6b  >  60  >  6C  >  6Xind  >  9Xo) 


Angular 

Separation 

Technique 

Applicable  Sce¬ 
nario 

Effects  Averaged 

6  A  ~ 

Non-tracked 
Multiple  Txs 

G-to-G 

AOA  and  beam 
wander,  phase, 
and  amplitude 

6  b  >  60 

Tracked, 

Multiple  Txs 

A-to-A,  A-to-G, 
G-to-A 

phase, 

amplitude 

6c  ~  6Xind 

Tracked, 

Multiple  Txs 

A-to-A,  A-to-G, 
G-to-A 

amplitude 

N/A 

Tracked,  Parallel 
Multiple  Txs 

A-to-A,  A-to-G, 
G-to-A 

amplitude 

For  a  mobile  Tx  and/or  Rx,  the  beams  must  be  tracked.  For  these  tracked-beam 
cases  in  Fig.  15c  air-to-air  path,  Fig.  15(d)  air-to-ground  path,  and  Fig.  15e  ground-to- 
air  path,  separations  beyond  the  isoplanatic  angle  up  to  approximately  the  tilt  isopla- 
natic  angle  should  average  out  the  higher-order  phase  effects.  Separations  larger  than 
the  tilt  isoplanatic  angle  require  separate  trackers.  This  occurs  for  longer  propaga¬ 
tions  and  near-transmitter  turbulence,  since  the  phase  tilt  effects  are  large  due  to  the 
long  lever  arm  of  the  turbulence.  Small  isoplanatic  and  tilt-isoplanatic  angles  occur 
for  propagations  longer  than  100  km  in  Fig.  15c  and  for  the  ground-to-air  propagation 
shown  in  Fig.  15e.  Fig.  15d  for  the  air-to-ground  scenario  shows  that  the  correlation 
angles  get  smaller  as  the  propagations  get  longer,  but  as  the  transmitter  altitude  gets 
above  the  turbulence  at  about  12  km,  the  angles  remain  relatively  similar  [46]. 


4-2  Simulation  set-up  and  validation 

Simulations  of  different  scenarios  and  separation  distances  show  how  much  mul¬ 
tiple  transmitters  improve  BER  performance.  The  turbulence  effects  explored  in 
simulated  scenarios  were  represented  by  random  optical  field  screens  with  the  correct 
statistics  placed  along  the  path.  The  layers  for  this  research  were  chosen  to  simu- 
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Ground  Range  (km) 
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Propagation  length  (km) 
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Figure  15:  (a)  Phase  isoplanatic  angle  (0O),  scintillation  correlation  angle  (0Xc),  tilt 
isoplanatic  angle  (Ota),  and  phase  uncorrelated  angle  (0^ind)  are  shown 
for  a  receiver  diameter  of  Dr  =  20  cm  for  HV-57  profile,  (b)  Horizontal 
propagation:  altitude  h  —  1  m,  L0  =  40  cm,  (72  =  10~14  m-2/3, 

and  L  =  0  to  10  km.  (c)  Horizontal  propagation:  altitude  h  —  10  km, 
L0  =  100  km,  C2  =  10-17  m-2/3,  and  L  =  0  to  300  km.  (d)  Air-to-ground 
path:  Transmitter  height  Htx  =  4  to  20  km,  £  =  70°,  and  receiver  height 
Hrx  =  0  km.  (e)  Ground-to-air  path:  Htx  =  0  km,  £  =  70°,  and  Hrx  =  4 
to  20  km  [46]. 
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Table  2:  Atmospheric  parameters  for  the  scenarios  used  in  the  BER  calculations. 


Scenario 

rovl  (cm) 

rosvh  (cm) 

^ sph 

60  (prad) 

0TA  (//rad) 

Ao 

G-to-G 

2.5 

4.5 

1.08 

0.437 

3.5 

11.7 

3.2 

G-to-A 

10 

85 

0.911 

0.0461 

1.2 

11 

A-to-A 

23 

41 

0.384 

0.155 

1.3 

3.0 
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late  the  continuous  model  so  that  several  low-order  moments  of  the  layered  model 
accurately  match  the  continuous  model  [46]. 

4-2.1  Simulation  set-up.  In  this  research,  10  random  phase  screens  were 
used  to  model  the  turbulence  along  varying-turbulcnce-strength  paths  and  five  screens 
along  the  constant-turbulence-strength  paths.  The  layered  analytic  planar  and  spher¬ 
ical  coherence  diameter  vq,  planar  and  spherical  Rytov  numbers  7 Z,  and  isoplanatatic 
angle  6q  matched  within  1%  of  the  full  path  continuous  atmospheric  turbulence  pa¬ 
rameters.  Table  2  summarizes  the  atmospheric  parameters  for  the  simulations  used  to 
calculate  the  BER.  In  the  simulations  that  follow,  a  Gaussian  beam  with  a  1/e  field 
radius  Wo  =  2.5  cm  and  A0  =  2L/(kWo)  propagates  to  the  receiver  aperture.  An¬ 
drews  and  Phillips  call  beams  with  A0  >  100  approximately  spherical  and  A0  <  0.1 
approximately  planar.  Therefore,  the  equations  in  previous  sections  where  a  point 
source  or  spherical  wave  are  used  are  a  reasonable  approximations,  especially  for  the 
air-to-air  propagation.  Earlier,  the  von  Karman  turbulence  power  spectrum  was  used 
to  model  the  phase  effects  [46]. 

Andrews  and  Phillips’  modified  turbulence  power  spectrum  is  used  in  the  sim¬ 
ulations  performed  in  this  research  because  it  includes  L0  and  Iq,  and  gives  the  best 
agreement  with  collected  atmospheric  data  for  phase  and  amplitude  effects  [2],  The 
Hufnagel- Valley  turbulence  prohle  was  used  in  this  research  with  the  parameters 
set  to  the  HV-57  moderate  turbulence  strength  (i.e.  turbulence  at  the  ground  is 
A  =  1.7  x  10-14m-2/3  and  the  effective  wind  at  altitude  is  W  =  21  m/s)  [46]. 

FS-based  phase  screens  were  used,  since  it  allows  for  better  low  spatial  frequency 
representation  than  other  techniques  [45,48,66,67,83].  This  modal  phase  screen 
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is  defined  in  terms  of  the  spatial  frequencies  of  the  refractive-index  PSD  and  need 
only  be  evaluated  at  the  grid  points  of  interest.  See  Section  2.4  for  details.  This 
approach  is  particularly  effective  when  calculating  the  fields  of  widely  spaced  beams 
over  long  periods  of  time.  Although,  in  this  work  independent  random  realizations 
of  turbulence  are  used  in  the  simulations  by  assuming  the  turbulence  is  an  ergodic 
random  process  [46]. 

Split-step  Fresnel  propagations  were  performed  for  the  Wo  =  2.5  cm  collimated 
Gaussian  beam.  Great  care  was  taken  to  adequately  sample  the  Fresnel  propagation 
between  the  screens  as  well  as  the  turbulence  effects  as  the  beam  propagates.  First, 
the  geometric  constraints  must  be  satisfied  to  avoid  aliasing  in  the  region  of  interest. 
Then  one  must  adequately  sample  the  quadratic  phase  term  in  the  Fresnel  propagation 
so  that  the  phase  in  that  same  of  region  of  interest  is  not  corrupted  by  aliasing.  This 
required  satisfying  these  four  sampling  constraints: 


s  .  AL  —  DRx5\ 

5n< 

(163) 

n>DTx  DRx  XL 

~  25!  25n  2 

(164) 

L  A  z 

(165) 

^<[1  +  ^1-—, 
n  l>Tx 

XL 

(166) 

N  —  T~T“  ’ 

0l0n 

where  8i  is  the  sampling  size  in  the  transmitter  plane  and  8n  is  the  sampling  size  in  the 
receiver  plane  [17].  Additionally,  DRx  is  the  receiver  diameter,  Dtx  is  the  transmitter 
diameter,  L  is  the  propagation  length,  R  is  the  beam’s  radius  of  curvature,  and  A  is  the 
laser  wavelength.  The  fourth  and  most  restrictive  constraint  listed  above  Eq.  (166) 
was  satisfied  by  performing  multiple  partial  propagations  (sometimes  referred  to  as 
split  step  propagations)  to  propagate  the  full  distance.  Finally,  the  turbulence  must 
be  adequately  sampled  and  the  wandering  and  beam  spreading  must  be  taken  into 
account  by  choosing  large  enough  screens. 
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Figure  16:  The  plot  shows  the  irradiance  and  the  phase  of  a  Gaussian  beam  with 
Wo  =  2.5  cm  after  propagating  100  km  with  no  turbulence  using  the 
propagation  code,  compared  against  the  analytic  solution. 

4-2.2  Validation  of  simulation  results.  Validating  simulations  is  important 
to  ensure  the  calculations  and  simulations  accurately  predict  the  diffraction  and  atmo¬ 
spheric  turbulence  effects.  Figure  16  shows  the  irradiance  and  phase  of  a  Wq  =  2.5  cm 
e~l  field  radius  Gaussian  beam  after  propagating  100  km  through  vacuum  using  the 
propagation  code  compared  with  the  analytic  solution.  The  two  discontinuities  in  the 
derivative  of  the  phase  show  the  boundary  of  accurate  phase  representation.  Beyond 
that  width  there  is  aliasing  in  the  quadratic  phase  factor,  but  it  is  beyond  the  region 
of  interest  and  does  not  affect  the  signal  received  by  the  20  cm  aperture.  Now  that 
correct  propagation  is  demonstrated,  the  atmospheric  effects  are  verified.  The  struc¬ 
ture  function  of  each  of  the  phase  screens  was  consistent  with  the  theoretical  values. 
Figure  14  shows  the  structure  function  averaged  over  40  screen  realizations  compared 
to  the  theoretical  value  for  different  turbulence  strengths. 

The  long-term  spot  size  Wlt  is  measured  by  averaging  the  spot  size  of  a  long 
period  of  time  or  over  many  iid  realizations.  This  spot  size  was  consistent  with 
theory  for  all  scenarios.  For  the  ground-to-air  propagation  the  long  term  spot  size 
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Figure  17:  Structure  functions  of  three  representative  screens  averaged  for  40  random 
screen  realizations.  The  analytic  structure  functions  are  plotted  with  solid 
lines. 

was  Wlt  =  67  cm  using  the  propagation  code,  compared  to  a  theoretical  spot  size 
using  Eq.  (46)  of  62  cm. 

The  measured  scintillation  is  consistent  with  the  theoretical  values  for  the  spher¬ 
ical  and  planar  Rytov  numbers.  The  theoretical  log-amplitude  variance  is  0.046  for  a 
spherical  wave  and  0.91  for  a  planar  wave  for  the  0  to  10  km  in  altitude  propagation 
at  a  70°  zenith  angle.  The  centroid-tracked  scintillation  index  for  the  center  pixel 
was  0.073.  For  the  untracked  case  the  scintillation  index  was  0.14  which  more  closely 
matches  the  theoretical  value  of  4  x  0.046  =  0.18  [46]. 

For  the  0  to  10  km  70°-zenith-angle  propagation  path  a  simulated  point  source 
was  propagated  to  the  receiver  and  the  statistics  were  calculated.  The  scintillation 
index  for  the  untracked  case  was  0.12.  The  scintillation  index  is  a  little  lower  than  the 
theoretical  value  of  4x0.046  =  0.18.  This  might  be  due  to  an  aperture  averaging  effect 
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caused  by  the  resolution  of  the  optical  field  at  the  receiver.  These  results  confirm  the 
simulation  operates  as  expected  and  should  adequately  model  the  turbulence  [46]. 

For  each  propagation  length,  the  irradiance  and  phase  of  the  Gaussian  beams 
after  propagating  through  a  vacuum  matched  the  analytical  solution.  The  structure 
function  of  each  of  the  phase  screens  was  also  consistent  with  the  theoretical  values. 
The  scintillation  index  at  the  receiver  for  a  simulated  point  source  was  consistent 
with  the  Rytov  approximation  for  scintillation.  For  each  turbulence  simulation,  the 
measured  long-term  spot  size  was  consistent  with  the  analytic  spot  size.  These  re¬ 
sults  confirmed  the  simulation  operated  as  expected  and  should  adequately  model  the 
turbulence  [46]. 

4-3  Modeling  optical  receiver  signals  and  noise  sources 

The  received  signal  is  converted  from  optical  power  to  receiver  current  in  the 
detector.  The  receiver  current,  which  is  called  the  measured  current  im,  is  the  sum  of 
the  signal  current  and  noise  current  and  is  given  by 

im  =  is  +  iN-  (167) 

Since  simulations  are  used  to  “measure”  the  performance,  it  is  imperative  to  model 
these  realistically.  First,  the  signal  current  can  be  represented  by  [18, 19] 

rjqs 
%s  =  ~h^ 

where  r)  is  quantum  efficiency  (electrons/photon),  q  =  1.602  x  10~19  C  is  the  ele¬ 
mentary  charge,  s  is  optical  power  at  the  detector  (Watts),  h  =  6.626  x  10-34  J-s  is 
Planck’s  constant,  and  v  is  optical  frequency  (Hz).  The  product  hv  gives  the  energy 
in  joules  of  a  single  photon.  The  noise  current  i^  is  modeled  as  a  zero-mean  Gaussian 
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random  variable  with  a  variance  o\  equal  to 


E  K]  =  (T1  =  °2elec  +  °2shot  +  °ASE>  (169) 

where  aflec  is  the  thermal  (Johnson)  noise,  cr2shot  is  the  shot  noise,  and  cr\SE  is  the  ASE 
noise  associated  with  an  EDFA.  Since  each  of  these  noise  current  sources  are  inde¬ 
pendent  the  overall  noise  current  variance  erf  is  the  sum  of  each  noise  source  variance. 
Modeling  receiver  noise  sources  is  essential  to  accurately  represent  a  communication 
system.  Two  fundamental  noise  sources  for  optical  receivers  are  the  signal-level- 
dependent  shot  noise  and  the  temperature-dependent  Johnson  noise.  Shot  noise  cr2shot 
is  fundamentally  a  Poisson  random  process  as  the  photo-electrons  are  generated  for 
random  arrivals  of  photons.  The  number  of  photons  per  bit  is  well  over  100  in  all  cases 
studied  here.  Furthermore,  as  the  captured  power  approaches  zero,  the  shot  noise  ap¬ 
proaches  zero  much  faster  than  the  signal.  Therefore,  shot  noise  current  statistics  are 
well  approximated  as  a  zero- mean  Gaussian  random  process  [1,3].  The  mean-square 
current  due  to  shot  noise  is  given  by  [19] 

E  [ilhot]  =  a2hot  =  2qisB=  2mh*  B ,  (170) 

where  B  is  the  electrical  bandwidth.  The  signal  current  is  assumed  to  be  constant 
during  any  given  integration  period  corresponding  to  a  single  bit.  Comparing  the 
frequency  of  atmospheric  change  (~  kHz)  with  bandwidths  studied  (3>  MHz),  this  is 
a  reasonable  assumption. 

Johnson  (thermal)  noise  current  is  typically  modeled  as  a  zero-mean  Gaussian 
random  variable,  with  mean-square  current  determined  by  [19] 

r-9  i  9  4  KTB 

E  [*eZec]  =  ®  elec  =  >  (l?l) 
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where  K  =  1.381  x  10~23  J/K  is  Boltzmann’s  constant,  T  is  temperature  of  the 
electronics  (K),  and  R  is  the  effective  input  resistance  (ohms).  These  are  the  primary 
noise  sources  in  the  receiver. 

There  is  also  noise  due  to  the  type  of  amplifier  or  gain  mechanism.  In  this 
research,  the  received  optical  signal  is  coupled  into  an  optical  fiber  to  be  amplified 
by  an  EDFA.  The  EDFA’s  advantage  is  the  capability  to  achieve  high  gain  at  very 
high  bandwidths.  In  addition,  the  gain  in  an  EDFA  saturates,  affording  some  gain 
control  to  reduce  optical  signal  variation.  This  reduces  the  required  dynamic  range 
of  the  detector.  The  EDFA’s  noise  is  due  to  amplified  stimulated  emission  and  can 
be  modeled  as  a  signal-dependent  zero-mean  Gaussian  noise  source  given  by  [1] 

p 

E  Vase\  =  Vase  =  M2nspriinr]2outG(G  -  1  )-^B 

=  4nspr]outq(G  -  1  )iaB,  (172) 

where  G  is  the  gain,  r]in  and  r]out  are  the  input  and  output  losses,  and  nsp  is  the 
spontaneous  emission  factor.  Other  gain  mechanisms  like  avalanche  photo  diodes 
(APD)  are  limited  to  about  100-200  GHz  gain-bandwidth  product,  whereas  an  EDFA 
with  a  PiN  photodiode  would  be  at  least  an  order  of  magnitude  higher  [1],  For  lower 
bandwidths  an  APD  is  advantageous  due  to  a  much  higher  coupling  efficiency. 

4-4  Simulation  approach 

The  receiver  consists  of  a  20-cm-diameter  lens  with  a  1  m  focal  length  focused 
onto  a  single-mode  optical  fiber.  The  fiber  core’s  diameter  was  3  pm  and  the  numerical 
aperture  (NA)  was  0.20,  consistent  with  a  commonly  available  EDFA.  For  the  100  km 
propagation,  the  collimated  Gaussian  beam  at  the  Rx  is  much  larger  than  the  aperture 
and  since  r$>  D  the  spot  size  at  the  focal  plane  of  the  lens  is  determined  by  2.44A//D, 
limited  only  by  diffraction.  For  the  ground-to-ground  and  ground-to-air  cases  the 
turbulence- induced  spot  size  is  approximately  2.44/A/ro,  since  r0  <  D  for  those 
cases.  See  Table  2  and  Fig.  18. 
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Figure  18:  When  D  <  ro,  the  spot  size  is  determined  primarily  by  diffraction 
2.44A//.D.  Where  as  when  D  >  r0,  the  spot  size  spreads  primarily  due  to 
the  turbulence,  leading  to  a  much  wider  spot  determined  by  2.44A  f/r0  [46]. 


For  the  angularly-separated-beam  scenario,  two  Gaussian  beams  are  displaced 
half  the  separation  distance  in  opposite  y-directions  and  a  linear  phase  is  applied  to 
“aim”  both  beams  at  the  center  of  the  receiver  aperture.  For  the  parallel-beam  case, 
the  two  Gaussian  beams  are  displaced,  and  each  beam  remains  off-axis  by  half  the 
separation  distance. 

A  coarse  tracking  system  was  simulated  for  the  transmitter  and  receiver  for 
the  ground-to-air  and  the  air-to-air  paths  by  implementing  an  ideal  centroid  tracker 
and  adding  random  tracking  system  errors.  The  errors  in  the  Tx  tracker  cr|  are 
driven  by  tilt  isoplanatism  aPA  due  to  the  point-ahead  angle,  temporal  errors  aPT  in 
the  controller,  platform  jitter  <jpj,  and  measurement  error  crPM.  Therefore  the  total 
transmitter  tracker  error  is  [31,46] 


aj  —  <JPA  +  CTPT  +  <7pj 


+  < 7 ' 


TM  i 


(173) 


measured  in  rad2  of  the  tilt  angle.  The  tilt  anisoplanatic  error  is  driven  by  the 
separation  angle  size  compared  to  the  tilt  isoplanatic  angle  and  given  by  [31] 
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in  rad2  of  the  tilt  angle  where  9  is  the  angle  between  the  propagation  path  and 
the  beacon  path.  This  tilt  anisoplanatic  error  requires  a  different  tilt  isoplanatic 
angle  OpAtirad)  definition  than  the  one  previously  mentioned  in  Eq.  (154).  This  new 
definition  of  9TA{irad)  is  the  angle  between  two  beams  at  which  the  wavefront  phase 
error  due  to  the  tilt  anisoplanatism  is  1  rad,  resulting  in 

0.1947A 

9tA(1  rad)  -  - Yj2 - 


The  required  point-ahead  angle  determines  the  inherent  angular  separation  of 
the  beacon  and  the  transmit  path.  The  point-ahead  angle  for  the  ground-to-air  path 
is  [31] 

Opa  —  —  sin  £,  (176) 

c 

where  Vp  is  the  platform  velocity  and  c  is  the  speed  of  light.  This  angle  can  be 
substituted  for  6  in  Eq.  (174).  The  temporal  error  is  [31,66] 


2  _  9 
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and  the  tilt  measurement  error  is  [31] 


2  _  9 
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(178) 


The  temporal  error  is  dependent  on  the  tracking  system  bandwidth,  and  the  tilt 
measurement  error  depends  upon  the  signal-to-noise  ratio  (SNR).  The  platform  jitter 
is  driven  by  the  transmitter  and  receiver’s  residual  vibrations.  The  receiver  tracker 
error  includes  the  last  three  terms  of  Eq.  (173). 


4-5  Simulation  results 

The  simulations  used  independent  random  realizations  of  the  appropriate  tur¬ 
bulence  statistics  to  determine  the  optimal  separations  for  a  two-transmitter  system 
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in  three  scenarios:  a  ground-to-ground  fixed  transmitter  and  receiver  4  km  link  with 
a  Rx  and  Tx  height  of  1  m,  a  ground-to-air  (h  =  1  m  to  10  km)  30.2  km  path  with 
a  zenith  angle  £  =  70°,  and  an  air-to-air  100  km  path  at  10  km  in  altitude  [46].  This 
approach  assumes  that  the  turbulence  can  be  approximated  by  an  ergodic  random 
process  where  the  temporal  statistics  can  be  approximated  by  the  ensemble  statistics. 

4-5.1  Differential  irradiance  and  tilt  variance.  A  number  of  performance 
measurements  were  calculated,  but  the  differential  irradiance  variance  between  the 
two  beams  o\irr  =  E  [(A  —  If)2]  —  E  [A  «  A]2  best  indicated  when  adequate  aver¬ 
aging  would  occur.  [5,60,68]  The  larger  the  differential  irradiance  variance,  the  less 
correlated  the  irradiance  fluctuations  become.  For  uncorrelated  beams  this  variance 
should  approach  two  times  the  irradiance  variance  of  a  single  beam.  If  the  angular 
separation  is  much  beyond  this  point,  the  power  received  at  the  detector  or  fiber  is 
reduced  due  to  the  difference  in  the  AOA  of  the  beams.  As  parallel  beams  move 
farther  off-axis,  power  reduces  and  the  variability  of  the  constituent  beams  increases. 
This  effect  increases  the  double-Tx  system  BER,  driving  the  performance  below  the 
single-Tx  case  [46]. 

A  4  km  ground-to-ground  propagation  was  performed  with  the  same  beam  pa¬ 
rameters  as  the  29.2  km  propagation  path.  This  beam  wandered  off  the  receiver 
frequently.  The  beam  walk-off  standard  deviation  was  0.095  m  while  the  short  term 
beam  half-width  was  0.083  m.  Again,  the  optimal  separation  for  this  scenario  was 
about  3 pc.  The  irradiance  variance  was  much  higher  than  the  other  scenarios  since  the 
beam  wandered  off  the  receiver  often.  In  Fig.  19a  the  differential  irradiance  variance 
approached  two  times  the  single  irradiance  variance  at  about  3 pc.  Figure  19b  shows 
that  in  the  y-direction  (parallel  to  the  separation)  the  two  beams  start  to  wander  in¬ 
dependently  for  narrow  separations  whereas  in  the  x-direction  (perpendicular  to  the 
beam  separation)  the  beams  wander  independently  at  approximately  d^ind/ 2.  The 
differential  tilt  variance  approached  two  times  the  single  tilt  variance  very  quickly 
in  the  direction  parallel  to  the  separation,  but  for  the  perpendicular  direction  the 
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Figure  19:  Ground-to- Ground  untracked  case.  Propagation  length  of  4  km,  at  1  m 
above  the  ground,  collimated  beam  with  Wo  =  0.025  m.  (a) This  plot 
shows  how  the  differential  irradiance  variance  approaches  2  times  the  single 
beam  irradiance  variance,  (b)  This  plot  shows  that  in  the  y-direction 
(parallel  to  the  separation)  the  two  beam  start  to  wander  independently  for 
narrow  separations  whereas  in  the  x-clirection  (perpendicular  to  the  beam 
separation)  the  beams  wander  independently  at  approximately  dpind/ 2  [46]. 

required  separation  distance  for  uncorrelated  phase  was  10  to  15  pc  or  about  d^ind/2. 
This  is  consistent  with  the  phase  independence  angle  which  includes  both  separation 
directions.  This  is  also  consistent  with  Sasicla’s  relations  in  Eqs.  (150)  and  (151), 
where  the  tilt  component  parallel  to  the  displacement  is  1.73  times  the  perpendicular 
component  for  small  displacements  and  approaches  equality  for  very  large  displace¬ 
ments  [46,69]. 

Next,  Fig.  20  shows  the  differential  irradiance  variance  and  the  differential  cen¬ 
troid  variance  for  the  air-to-air  scenario  shown  in  Table  2.  In  Fig.  20a  at  2  and  3  pc, 
the  differential  variance  reaches  a  maximum.  Fig.  20b  shows  the  differential  centroid 
variance  continues  as  the  separation  between  the  transmitters  increases  because  the 
outer  scale  L0  and  corresponding  dphiind  for  this  propagation  is  very  large  [46]. 


As  shown  in  Fig.  21,  beams  approached  uncorrelated  irradiance  variance  at 
about  2-3  pc  for  angularly  separated  beams  for  both  the  air-to-air  tracked  system  and 
the  ground-to-air  tracked  system.  In  addition,  for  the  air-to-air  scenario  the  amplitude 
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Figure  20:  100km  air-to-air  propagation  path  at  a  10  km  altitude  for  a  collimated 
Wo  =  2.5  cm  beam  (a)  differential  beam  irradiance  variance,  and  (b) 
differential  centroid  location  variance  [46]. 


structure  function  (shown  in  Fig.  14)  has  a  peak  before  settling  into  the  asymptotic 
value  of  two  times  the  the  single-beam  irradiance  variance  [46]. 


4-5.2  Bit  error  rates  for  independent  realizations.  The  BER  for  the  different 
scenarios,  tracking  systems,  and  separation  distances  are  shown  in  Figs.  22  -  25.  The 
BERs  were  calculated  using  the  PDF  of  the  measure  signal  current  (determined  by 
the  histogram  of  the  received  signal),  accounting  for  the  shot,  thermal,  and  ASE  noise. 
Since  the  shot  and  ASE  noise  are  signal-dependent,  their  variances  changed  for  each 
independent  realization,  while  the  thermal  noise  variance  was  fixed.  After  solving  for 
the  optimal  fixed  threshold,  the  probabilities  of  missed  detection  and  false  alarm  were 
determined.  The  total  power  in  the  single-Tx  system  was  1  Watt,  and  the  total  power 
in  the  double-Tx  system  was  also  1  Watt  (0.5  Watts  per  transmitter).  For  the  BER 
plots,  the  EDFA  gain  remained  constant  at  30  dB.  The  signal  level  differences  shown 
in  the  plots  vary  due  to  differences  in  propagation  attenuation,  coupling  efficiencies, 
transmitter  levels,  etc.,  but  do  not  depend  on  the  EDFA  gain  [46]. 

First,  the  BER  for  the  ground-to-ground  path  for  single  and  double-Tx  systems 
were  calculated.  No  tracking  system  was  used  in  this  case,  and  the  beams  walked  off 
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Figure  21:  Differential  irradiance  variance  for  two  angularly  separated  beams.  Irradi- 
ance  is  taken  from  the  center  of  the  untracked  beams,  separately  tracked 
beams,  and  a  single  combined  beam  tracker.  The  solid  blue  line  is  two 
times  the  variance  of  on-axis  irradiance  of  a  single  beam.  The  differen¬ 
tial  variance  approaches  two  times  this  value  as  the  separation  increases, 
(a)  Air-to-air  path  angular  separation  (b)  Ground-to-air  path  angular  sep¬ 
aration  [46]. 


of  the  receiver  often.  To  quantify  this,  note  that  the  beam  walk-off  standard  devi¬ 
ation  was  9.5  cm  and  the  short-term  beam  half-width  was  8.3  cm.  The  differential 
tilt  variance  approached  two  times  the  single  tilt  variance  in  the  direction  parallel  to 
the  separation  for  very  small  separations  (i.e.  3  cm).  However,  for  the  perpendic¬ 
ular  direction  the  required  separation  distance  for  uncorrelated  phase  was  about  10 
pc  or  about  dpinj2.  This  phenomenon  is  consistent  with  the  differential  tilt  relations 
presented  by  Sasiela  [69].  Figure  22  shows  slight  improvement  for  the  double  trans¬ 
mitter  case  over  the  single  beam  when  the  transmitters  propagated  in  parallel,  but 
no  improvement  for  angular  separations.  Due  to  the  inherent  non-zero  mean  AOA 
for  angularly  separated  beams,  the  focal  spots  of  both  beams  on  average  miss  the 
fiber,  reducing  the  coupling  efficiencies.  If  the  fiber  core  is  small,  as  in  this  case,  this 
could  severely  limit  the  coupling  of  both  beams.  Whereas  with  the  parallel  beams, 
both  focal  spots  were  on  average  on  the  center  of  the  fiber,  allowing  for  much  better 
coupling  [46]. 


Power/Loss  Factor  (dB) 

(a) 


(b) 


Figure  22:  These  plots  show  the  BER  for  a  ground-to-ground  link.  In  plot  (a)  the 
beams  were  angularly  separated  and  in  plot  (b)  the  beams  were  separated, 
but  traveled  in  parallel  [46]. 


Next,  the  BER  plots  for  the  tracked  ground-to-air  scenario  in  Fig.  23  show  that 
there  is  an  improvement  afforded  by  using  two  transmitters  of  about  3  dB  for  the 
ideal  tracker.  Interestingly,  two  transmitters  also  improved  performance  for  the  non¬ 
ideal  case  for  a  tracker  error  of  c jj  =  A/(4 D).  Parallel  beams  were  used  with  a  center 
tracker  system.  For  this  case,  Fig.  24  shows  an  improvement  again  for  a  tracking 
system  with  a3  =  A/(4 D).  This  improvement  reduces  in  both  cases  as  the  tracking 
system  performance  degrades  [46]. 

Finally,  the  BER  plots  in  Fig.  25  for  the  air-to-air  100  km  path  shows  the 
best  improvement  for  a  separation  distance  of  2  to  3  pc.  This  is  consistent  with 
the  differential  scintillation  measurements.  The  largest  improvement  (approximately 
4  dB)  occurs  for  the  finest  tracking  system.  As  the  tracking  degrades,  the  improvement 
due  to  the  signal  diversity  decreases. 

For  all  cases  studied  here,  a  separation  of  a  small  multiple  of  pc  was  sufficient 
to  provide  adequate  averaging.  These  separations  were  consistent  with  the  coherence 
width  described  by  Andrews  and  Phillips  [2],  The  amplitude  independence  separa¬ 
tion  distance  A xXind  and  angle  9Xind  for  these  scenarios  are  approximately  2 pc  and 
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Figure  23:  Bit  error  rate  for  a  groimd-to-air  link  with  angularly  separated  beams 
using  various  tracking  systems  (a)  ideal  centroid  tracker,  (b)  aj  =  A/(8 D), 
(c)  aj  =  A/(4D),  and  (d)  aj  =  3 A/ (8 D)  [46]. 
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Figure  24:  Bit  error  rate  for  a  ground-to-air  link  with  parallel  separated  beams  us¬ 
ing  various  tracking  systems  (a)  ideal  centroid  tracker,  (b)  Uj  =  A/(8 D), 
(c)  <7j  =  A/(4 D),  and  (d)  aj  =  3 A/ (8 D)  [46]. 
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Figure  25:  BER  for  Air-to-air  100km  path  at  10km  altitude  using  various  tracking 
systems  (a)  ideal  centroid  tracker,  (b)  aj  =  A/(8 D),  (c)  aj  =  A/(4Z)), 
and  (d)  <jj  =  3A/(8 D)  [46]. 
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Tabic  3:  Optimal  separation  distances  for  multiple  Tx  systems. 


Scenario 

Tx  Alt 

Rx  Alt  (km) 

Prop  Dist 

Best  Tx 

Separation 

G-to-G  Angular 

1  nr 

1  nr 

4  knr 

not 

conclusive 

G-to-G  Parallel 

1  nr 

1  nr 

4  knr 

2.5  pc  = 

7.9  cm 

G-to-A  Angular 

1  nr 

10  knr 

29  knr 

Pc 

8.6  enr 

G-to-A  Parallel 

1  nr 

10  knr 

29  knr 

Pc/2 

4.3  cm 

A-to-A  Angular 

10  km 

10  knr 

100  knr 

2pc  = 

31  cm 

2 (Lk)  1/2,  respectively.  Table  3  shows  the  optimal  separation  distances  for  the  three 
scenarios  tested  in  this  section. 

The  optimal  separation  for  the  air-to-air  100  km  path  scenario  at  10  km  in 
altitude  was  only  31  cm.  The  scintillation  for  this  air-to-air  scenario  is  relatively 
high  at  Tl.Sph  =  0.155  whereas  the  phase  effects  characterized  by  D/r0sph  =  0.49  are 
relatively  low.  This  indicates  that  a  technique  to  reduce  scintillation  like  the  multiple- 
transmitter  system  tested  here  in  this  chapter  has  advantages  over  AO  systems  for 
the  air-to-air  scenario.  The  next  chapter  studies  this  air-to-air  scenario  further  by 
implementing  a  time-series  simulation  to  study  the  temporal  effects  of  multiple  trans¬ 
mitters.  Chapter  V  also  studies  the  trickle-down  effects  of  additional  techniques  to 
include  centroid  tracking  and  adaptive-threshold  systems. 
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V.  Temporal  Analysis  and  Signal  Processing 


This  research  shows  how  multiple  transmitters  reduce  received  signal  variability 
and  the  length  of  a  fade  caused  by  long  propagation  paths  through  extended 
turbulence  in  FSOC  systems.  Many  researchers  have  studied  how  adaptive  optics 
systems  can  improve  FSOC  system  performance  especially  in  satellite  communica¬ 
tions  [7,77,81].  These  conventional  AO  systems  correct  for  the  phase  only  and  cannot 
correct  for  strong  scintillation,  but  here  multiple  transmitters  “average  out”  strong 
scintillation  effects  by  incoherently  summing  up  multiple  beams  at  the  receiver.  Oth¬ 
ers  have  done  research  on  MIMO  and  multiple  transmitter  FSOC  systems  [30,43,60], 
but  have  not  delved  into  the  anisoplanatic  separations  required  and  the  temporal 
considerations  of  tracking  systems  and  fade  statistics  for  the  airborne  regime. 

Here,  two  mutually  incoherent  laser  beams  angularly  separated  by  the  irra- 
diance  independence  (anisoplanatic)  angle  [46]  average  the  scintillation  effects  and 
reduce  the  fade  length,  depth,  and  number  per  second.  Wave-optics  simulations  show 
that  a  combination  of  transmitter  diversity,  adaptive  thresholding,  and  tracking  sys¬ 
tems  significantly  reduce  the  BER  in  air-to-air  FSOC.  The  reduction  in  fade  length 
indicates  multiple  beams  might  provide  even  greater  improvement  when  coupled  with 
interleaving  and  FEC  coding. 

5.1  Spatial  statistics 

The  spatial  statistics  of  the  turbulence  effects  are  required  to  determine  the 
simulation  parameters.  This  research  simulates  a  100  km  air-to-air  path  and  isolates 
the  phase  (spherical  coherence  diameter  r0sp)  and  amplitude  (spherical  Rytov  num¬ 
ber  lZSp)  effects  by  adjusting  the  turbulence  profile.  The  !Zsp  is  equal  to  the  spherical 
log-amplitude  variance  cr'^  sp  for  weak  turbulence  [69].  The  parameters  were  chosen 
to  emulate  this  air-to-air  horizontal  scenario  with  aircraft  velocities  between  56  and 
280  m/s  and  altitudes  between  4  and  15  km.  These  parameter  ranges  were  chosen  be¬ 
cause  of  the  results  in  Fig.  25  of  Section  4.5  that  show  a  relatively  small  separation  of 
about  31  cm  (2 pc)  is  required  to  average  scintillation  effects  for  the  air-to-air  scenario. 
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Table  4:  Atmospheric  parameters  for  the  scenarios  used  in  the  simulations.  The  op¬ 


tical  wavelength  was  A  =  1.55  //rri. 


Scenario 

Alt.  (km) 

Velocity  (m/s) 

At  (/xs) 

Tlsp 

fa  (Hz) 

ro3P  (cm) 

1 

15 

113 

133 

0.0596  (L) 

117  (L) 

74  (H) 

2 

4 

56 

137 

0.1827  (H) 

113  (L) 

38  (L) 

3 

N/A 

280 

30 

0.0854  (L) 

518  (H) 

38  (L) 

4 

N/A 

280 

30 

0.1979  (H) 

518  (H) 

55  (H) 

5 

4 

225 

30 

0.1827  (H) 

518  (H) 

38  (L) 

The  experimental  design  fully  investigates  the  different  scenarios  and  conditions  of 
an  air-to-air  scenario.  Since  r0ap,  lZsp,  and  fo  adequately  describe  the  spatial  and 
temporal  turbulence  effects,  the  simulated  conditions  consist  of  a  one  half  fractional 
factorial  of  these  three  factors.  Designing  the  test  in  this  way  enables  the  determina¬ 
tion  of  the  primary  driving  factors  for  fades  and  bit  errors.  Table  4  summarizes  the 
atmospheric  parameters  for  the  simulations.  There  are  five  different  scenarios  with 
different  altitudes  (km),  air  velocities  (m/s),  sampling  times  (/xs),  spherical  Rytov 
numbers  lZsp,  Greenwood  frequencies  fc  (Hz),  and  spherical  coherence  diameters  r0ap 
(cm). 

The  spatial  statistics  of  the  turbulence  effects  also  determine  how  far  apart  the 
transmitters  must  be  to  get  good  averaging.  The  farther  the  two  transmitters  are 
separated,  the  more  uncorrclated  the  effects  become  (i.e.  anisoplanatic).  This  is 
important  when  multiple  transmitters  are  used  to  average  out  the  turbulence  effects. 
Next,  the  anisoplanatic  phase  and  amplitude  effects  are  considered. 

5.1.1  Anisoplanatic  effects.  For  a  multiple-Tx  system,  the  phase  and  am¬ 
plitude  fluctuations  decorrelate  as  the  beam  separation  increases  in  distance  or  angle. 
Starting  at  small  separations,  high-order  and  low-order  phase  effects  are  highly  cor¬ 
related  for  Tx  separations  less  than  or  equal  to  the  isoplanatic  angle  6$.  Most  of 
the  literature  defines  this  angle  as  when  the  phase  perturbation  structure  function  is 
less  than  or  equal  to  unity,  leading  to  Eq.  (137)  in  Chapter  IV  [26].  Applying  the 
maximum  value  of  the  phase  perturbation  structure  function  determines  the  phase 
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independence  angle.  At  this  separation,  the  high-order  and  low-order  phase  effects  of 
multiple  beams  are  relatively  uncorrelated,  spreading  and  wandering  independently. 
This  angle  was  defined  in  Eq.  (139)  and  shown  here  again  as  [46,49] 

0*^  =  (m) 

where  cr2^  is  the  phase  variance  for  a  plane-wave  source  and  a  point  receiver,  given 
in  Eq.  (140). 

As  for  the  amplitude  effects,  the  correlation  width  pcw  is  often  used  to  determine 
how  large  receivers  need  to  be  to  provide  some  degree  of  aperture  averaging  of  the 
scintillation  effects.  The  correlation  width  is  defined  as  the  1/e2  point  of  the  normal¬ 
ized  irradiance  covariance  function  [2],  Since  pcw  for  weak  turbulence,  varies  between 
1  to  3  Fresnel  zones  ( L/k )1//2  depending  on  beam  size,  [2]  in  this  work  it  is  referred 
to  as  simply  pc  =  (L/k)1/2 .  In  recent  work,  the  principle  of  reciprocity  was  used  to 
illustrate  how  pc  could  determine  how  far  apart  transmitters  need  to  be  to  provide 
adequate  averaging  in  the  receiver  [46,49].  Due  to  AOA  considerations,  the  increase 
in  off-axis  irradiance  variance,  and  negatively  correlated  amplitude  effects  near  pcw, 
very  wide  separations  are  not  necessarily  the  optimal  configuration  [2,5,46,49,60]. 
Therefore,  this  research  separates  the  beams  angularly  by  2 6XC  =  2(Lk)~ 1/2,  or  a 
separation  of  2 pc  at  the  transmitters,  as  determined  in  Chapter  IV. 

5. 2  Temporal  considerations 

Thus  far,  spatial  statistics  were  used  to  describe  the  effects  of  atmospheric  tur¬ 
bulence.  In  this  section,  the  temporal  statistics  are  considered  to  determine  BER 
improvement  afforded  by  tracking  systems  and  adaptive  thresholding.  Taylor’s  frozen 
flow  hypothesis  states  that  the  turbulence  structure  is  essentially  frozen  as  it  moves 
across  the  propagation  path  for  small  time  intervals  [66].  The  random  screens  “scroll” 
across  the  laser  beam  at  different  points  along  the  path  to  generate  a  time-series  of 
the  turbulence  in  the  simulations  described  in  Sections  5.3  and  5.4. 


96 


5.2.1  Frequency  of  the  turbulence  effects.  The  first  thing  to  consider  when 
building  a  temporal  simulation  is  the  sampling  frequency.  One  approach  is  to  de¬ 
termine  the  frequency  of  the  irradiance,  given  by  firr  =  1/7W,  passing  over  the  Rx 
aperture  where  Tjrr  is  the  characteristic  time  of  the  irradiance.  Over  the  time  Tjrr 
the  turbulence  evolves  so  that  the  effects  are  only  slightly  different  than  the  previous 
time  slice  (isoplanatic  time  difference).  The  average  time  it  takes  one  anisoplanatic 
irradiance  patch  to  scroll  across  the  aperture  is 

T*“  “  wm  “  (180) 

where  |U(z)|  is  the  transverse  velocity  of  the  turbulence  and  fa  is  the  Greenwood  fre¬ 
quency  (characteristic  atmospheric  frequency  defined  in  terms  of  the  3dB  bandwidth 
for  a  high-order  phase  controller)  [66,69].  The  irradiance  at  the  edges  of  this  time 
period  rXind  are  anisoplanatic.  That  is  to  say,  after  this  time  period  the  irradiance  is 
relatively  uncorrelated  with  the  previous  time  slice. 

These  simulations  use  a  conservative  estimate  of  Tirr  to  ensure  they  include  all 
potential  signal  variations.  Using  fa  as  a  reference  and  varying  the  temporal  sampling 
frequency  of  the  simulations  enables  the  determination  of  an  adequate  sampling  rate. 
The  received  signal  power  spectral  density  (PSD)  for  increasingly  finer  resolutions  is 
calculated  until  the  PSD  is  relatively  similar  over  20  dB  down  from  the  maximum 
value.  This  determination  is  shown  graphically  in  Fig.  26.  These  PSDs  were  consistent 
for  different  random  realizations.  The  resulting  sampling  frequency  is  fSirr  =  64 fa. 
At  first,  it  seems  like  this  frequency  is  very  high,  but  the  Greenwood  frequency  is  a 
measure  of  the  higher-order  phase  temporal  variations.  Amplitude  variations  occur  at 
much  smaller  scale  than  phase  effects  and  therefore  require  finer  spatial  and  temporal 
resolutions  [2,46].  These  results  were  determined  for  a  100  km  path  with  a  turbulence 
velocity  of  225  m/s  (503  miles/hr)  and  fa  =  518  Hz.  Temporal  plots  of  the  signal 
variation  appear  smooth  to  the  eye,  indicating  again  that  the  sampling  is  adequate 
(see  Fig.  29). 
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How  does  this  temporal  separation  Tjrr  relate  to  the  spatial  correlation  of  the 
irradiance  in  terms  of  pcl  To  adequately  sample  the  irradiance  requires  sampling  at 
the  Nyquist  frequency,  which  requires 

fsirr  =  2  x  firr.  (181) 


Since  fSirr  =  64/c  is  known  to  be  adequate  from  the  PSD  analysis  in  Fig.  26,  it  is 
instructive  to  determine  what  fraction  of  pc  is  adequate  to  satisfy 

Pc/X 

Tirr~\V{z)[ 

Solving  for  x  in  the  above  equation,  results  in 


x  = 


PcJ  Sirr 

TO I 


12. 


Finally,  this  yields  the  following  irradiance  time  constant 


(183) 


pj  12 

"  in*)  r 


(184) 


This  relation  can  be  used  to  determine  the  sampling  rate  fs  in  simulations  with  sig¬ 
nificant  irradiance  fluctuations  without  having  to  calculate  the  PSD  of  the  receiver 
irradiance.  The  separation  pc/ 12  is  shown  in  Fig.  27  using  the  structure  functions 
previously  shown  in  Fig.  14.  The  separation  of  pj 12  could  be  considered  the  scintilla¬ 
tion  isoplanatic  separation  distance  for  angular  paths,  corresponding  to  a  scintillation 
isoplanatic  angle  of  9Xo  =  (pc/12)  /L. 


5.2.2  Threshold  determination.  For  a  binary  symbol  system  like  the  one 
used  here,  once  the  signal  is  received  a  decision  must  be  made  based  on  a  threshold 
whether  a  ‘1’  or  a  ‘O’  was  sent.  The  transmitter  modulates  with  OOK,  where  the  laser 
turns  on  to  transmit  a  ‘1’  and  turns  off  to  transmit  a  ‘O’.  The  transmission  of  a  ‘1’ 
or  ‘0’  is  equally  likely,  denoted  by  the  events  Hi  and  H0,  respectively.  The  likelihood 


Covariance  Power  Spectral  Density  Estimate 


Figure  26:  Received  signal  power  spectral  density  (PSD)  for  an  air-to-air  100  km  path 
at  4km  in  altitude.  The  PSD  of  the  signal  is  sampled  at  fs  —  64  x  fa  shown 
with  a  solid  line,  at  fs  —  32  x  fc  with  a  dotted  line,  and  at  fs  =  5.6  x  fa 
with  a  dashed  line.  The  horizontal  dot  dash  line  is  20  dB  down  from  the 
peak  value.  The  vertical  dashed  line  is  the  Greenwood  frequency  fa- 


99 


(p)  a 


Figure  27:  Phase  (right  scale)  and  log-amplitude  (left  scale)  structure  functions  are 
plotted  for  a  100  km  horizontal  propagation  at  10  km  altitude,  with  an¬ 
gularly  separated  beams.  The  strength  of  turbulence  is  L0/r0  =  286  and 
Lq/(\L)  =  23225.  [46]  Isoplanatic  and  anisoplanatic  separation  distances 
are  shown  for  the  amplitude  and  phase  effects.  The  amplitude  isoplanatic 
separation  distance,  pj  12  is  circled  with  a  dashed  line. 
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ratio  test  (LRT)  from  Eq.  (128)  determines  the  optimal  decision  threshold  based 
upon  the  discretized  PDF  of  the  measured  current  level  im  of  the  transmission  of  a  ‘1’ 
p(im\Hi)  and  transmission  of  a  ‘0’  p(im\Ho).  Using  the  LRT  and  the  assumption  that 
P(Hq)  =  P(Hi)  (equally  likely  signaling)  leads  to  the  following  two  relations;  [82]  if 

p(im\H{)  >  p(im\H0)  (185) 

then  the  algorithm  picks  H\  and  if 

p(im\Ht)  <p(im\H0)  (186) 


then  the  algorithm  picks  H0.  The  optimum  detection  criteria  can  best  be  described 
graphically  as  shown  in  Fig.  11  and  reported  in  Section  2.6.3.  These  turbulence 
conditions  vary  significantly  over  time  and  thus  could  benefit  from  a  threshold  that 
varies  with  the  optical  signal  level  [13,18]. 


5.2.2. 1  Fixed  Threshold.  The  fixed  threshold  calculation  takes  into 
account  the  PDF  of  the  signal  p(s)  due  to  variations  caused  by  channel  conditions. 
In  this  case,  the  channel  conditions  are  dictated  by  the  atmospheric  turbulence  con¬ 
ditions.  To  simplify  the  notation,  let  pi(im)  =  p{im\Hi)  and  po(im)  =  p(im\Ho)-  The 
LRT  for  this  scenario  is  [18, 72] 


P(H0)p0(im) 

1  (  -i\\ 

- eXP  o  2 

17  elec  K^elecJ 


P(Hi)  /  pi(im\s)p(s)ds 


r  p(*)  _ 

-  (ir  -  im(s))2 

/  /  \  cxp 

Jo  <u(s) 

2al(s) 

ds. 


(187) 

(188) 


The  threshold  %t  current  can  be  solved  for  numerically  whether  the  PDF  of  the  tur¬ 
bulence  induced  power  fluctuations  p(s)  is  analytic  or  calculated  from  the  discretized 
PDF  of  the  simulated  received  power  before  the  measurement  noise  is  applied.  Since 
this  is  a  simulation,  the  PDF  p(s)  is  calculated  from  the  raw  received  signal  before  the 
noise  is  applied.  The  noise  associated  with  detecting  a  ‘0’  is  primarily  due  to  thermal 
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noise  as  shown  in  Eq.  (171).  The  noise  associated  with  detecting  a  ‘1’  can  be  broken 
into  the  sum  of  the  thermal,  shot,  and  amplifier  noise,  defined  by 

°1  =  allec  +  alhot  +  aASEi  (189) 


where  cr2hot  is  the  shot  noise  due  to  the  random  arrival  of  photons  [see  Eq.  (170)]  and 
aASE  [See  Eq.  (172)]  is  the  amplified  stimulated  emission  noise  associated  with  an 
Erbium-doped  fiber  amplifier  [1,46].  The  probability  of  an  error  Pe  is  the  probability 
of  a  missed  detection  Pmd  plus  the  probability  of  a  false  alarm  Pf„; 

Pe  =  P{H\)Pmd  +  P(H0)Pfa  =  ^  (190) 

where 

P^  =  \L  erfc  (kSrf ) p(s)ds’  (191) 

Pfa  =  2erfc  ( '  (192) 
The  erfc(x)  is  the  complementary  error  function  in  Eqs.  (191)  and  (192). 


5. 2. 2. 2  Adaptive  threshold.  For  this  temporally  varying  turbulence,  an 
ideal  optimal  adaptive  threshold  shows  the  best  possible  BER  improvement.  Since  the 
threshold  is  determined  for  each  current  level,  the  PDF  of  the  received  signal  level  p(s) 
is  not  required  for  this  calculation.  Only  the  estimates  of  the  mean  and  variance  of 
the  two  conditions  are  required  to  set  the  threshold.  Solving  for  the  optimal  adaptive 
threshold  current  yields  [14, 18] 


ir 


ho^i  ~  RiQ'o 


+ 


cr0ai 

aj  -  a l 


ho)2  +  2  (erf 


(193) 


The  work  here  assumes  /i0  =  0  and  a0  =  a ezec,  since  <Jshot  =  &ase  =  0  when  a  ‘O’  is 
sent.  This  ideal  adaptive  threshold  system  calculates  the  optimal  adaptive  threshold 
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for  each  time  slice  with  the  corresponding  raw  received  signal  level  s  in  the  simulation 
and  implements  that  threshold  to  determine  whether  it  is  a  ‘1’  or  a  ‘O’. 

For  the  adaptive  threshold  case  the  probability  of  a  missed  detection  and  the 
probability  of  false  alarm  now  have  a  threshold  that  varies  with  the  signal  level  along 
with  all  of  the  other  signal  dependent  terms.  The  Pm(i  becomes 


(194) 


where  the  threshold  now  becomes  a  function  of  the  received  power  s.  The  Pfa  also 
becomes  a  function  of  s  given  by 


(195) 


In  a  real  system,  an  estimator  needs  to  be  implemented  to  determine  what 
threshold  It  is  used  for  the  next  particular  time  slice.  The  performance  of  this  es¬ 
timator  is  driven  by  the  measurement  noise  and  the  estimator’s  sampling  frequency. 
Since  the  mean  and  variance  of  the  transmission  of  a  ‘0’  are  relatively  constant  (for 
a  fixed  temperature),  Eq.  (193)  becomes  a  function  of  the  mean  and  variance  of  the 
signal  level  of  a  ‘1’.  Because  the  variance  o\  is  signal- dependent  and  the  signal  vari¬ 
ation  is  slow  compared  to  the  data  rate,  the  variation  in  the  adaptive  threshold  is 
only  a  function  of  the  signal  level  for  the  transmission  of  a  ‘1’.  In  addition,  since  the 
transmission  of  a  ‘1’  or  ‘0’  is  equally  likely,  the  mean  signal  level  for  the  transmission 
of  a  ‘1’  can  be  determined  by  ji\  ~  2 /J,rCvd  —  Po,  if  averaged  over  a  short  period  of 
time  with  respect  to  the  turbulence  (i.e.  r  <  Tjrr).  Therefore,  the  estimated  optimal 
adaptive  threshold  can  be  deduced  in  a  simulation  by  using  the  estimate  of  the  current 
signal  level  is  and  the  estimated  measurement  noises  allec  and  af. 

In  this  work,  the  received  power  is  split  into  two  branches  with  99  %  of  the 
power  used  in  the  digital  Rx  and  1  %  used  in  the  estimator.  The  estimator  measures 
the  current  in  the  previous  time  slice  (iEm-)  and  the  differential  signal  in  the  previous 
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Figure  28:  Adaptive  threshold  estimator. 


two  measurements  (i eui _ )  to  determine  the  estimated  signal  level.  The  differential 

of  the  measured  signal  Am=  _ )  determines  the  trend  in  the  previous 

two  estimator  measurements  and  _ ,  respectively,  to  further  refine  the  es¬ 

timate.  Figure  28  illustrates  how  it  operates.  In  these  simulations  the  temperature 
and  bandwidth  are  constant,  so  cb)ec  remains  constant.  The  estimated  current  is  is 
determined  by 


is  =  100[iB_  +  n(iE~)  +  Am]  (196) 

=  100[i^_  +  Aro],  (197) 

where  the  noise  in  the  measurement  n{iE-)  is  a  zero- mean  Gaussian  random  variable 
with  a  variance  equal  to 

°~1  —  allec  +  alhot  +  aASE-  (198) 

The  estimator  bandwidth  and  the  signal  level  drive  these  noise  sources  in  Eq.  (198), 
but  since  the  estimator  bandwidth  need  only  be  in  the  kHz  range  to  keep  up  with  the 
turbulence,  the  noise  power  is  relatively  low.  Reducing  bandwidth  of  the  estimator 
further  increases  the  latency  of  the  estimator  and  degrades  the  performance  of  the 
estimator.  If  /q  in  Eq.  (193)  is  set  to  equal  the  estimated  signal  is  and  /x0  ~  0,  the 
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equation  becomes 


<2  +  2  (t?  -  <*\ 


.2 

e/ec 


:)  In 


(199) 


where  Sq  is  the  estimate  of  ny  using  is. 

5.3  Simulation  set-up 

5.3.1  Modeling  the  turbulence.  The  turbulence  effects  explored  subsequently 
in  simulated  scenarios  were  generated  using  10  Fourier-series-based  random  phase 
screens  with  the  correct  statistics  placed  along  the  path  [48,83].  The  layered  ana¬ 
lytic  spherical  coherence  diameter  r0sp,  spherical  Rytov  number  !Zsp,  and  isoplanatatic 
angle  9q  matched  within  1%  of  the  full  path  continuous  atmospheric  turbulence  pa¬ 
rameters. 

The  simulations  propagate  a  collimated  Gaussian  beam  with  a  1/e  field  radius 
of  Wo  =  2.5  cm  using  a  split-step  Fresnel  propagation  to  a  20  cm  diameter  receiver 
aperture.  Great  care  was  taken  to  adequately  sample  the  Fresnel  propagation  between 
the  screens  to  avoid  aliasing  in  the  beam  as  well  as  the  quadratic  phase  term  [17].  See 
Section  4.2.1  for  the  constraints. 

5.3.2  Temporally  modeling  tracker  jitter.  A  coarse  tracking  system  was 
simulated  for  the  transmitter  and  receiver  for  the  five  air-to-air  paths  detailed  in 
Table  4  by  implementing  an  ideal  centroid  tracker  and  adding  random  tracking  system 
errors.  The  errors  in  the  Tx  tracker  are  driven  by  tilt  isoplanatism  crPA  due  to  the 
point-ahead  angle,  temporal  errors  a'PT  in  the  controller,  measurement  error  and 
platform  jitter  cr2PJ.  Therefore,  the  total  transmitter  tracker  error  is  [31] 


(200) 
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The  tilt  anisoplanatic  error  is  driven  by  how  close  the  beams’  separation  angle  is  to 
the  tilt  isoplanatic  angle.  The  receiver  tracker  error  includes  the  last  three  terms  of 
Eq.  (200). 

The  simulated  proportional-integral  (PI)  tracker  controller  attains  temporally 
correlated  jitter  rather  than  using  random  draws  to  simulate  tracker  jitter.  The  tilt 
corrector  simulates  the  tracker  by  centering  the  centroid  of  the  irradiance  at  the  pupil 
plane.  In  order  to  emulate  the  Tx  and  Rx  trackers,  temporal  errors  are  incorporated  in 
the  tilt  correctors  by  adjusting  the  loop  gain  of  the  feedback  controller.  The  simulated 
different  residual  beam  wander  cr6u)  and  beam  angle-of- arrival  (AO A)  jitter  ctaoa  are 
multiples  of  the  entire  beam  wander  a\-  for  the  Tx  and  AOA  jitter  ot  for  the  Rx, 
respectively.  These  runs  included  an  ideal  Tx  tracker  to  control  the  beam  wander  (i.e. 
each  centroid  is  shifted  to  the  center  of  the  aperture)  such  that  abw  =  0  and  non-ideal 
controllers  with  residual  beam  wanders  of  (7})W  =  0.25cx_y,  =  0.75 ax,  and  a^)W  =  ax. 
Likewise,  the  Rx  tracker  included  an  ideal  tracker  to  track  the  tilt  at  the  aperture  and 
the  resulting  focal  spot  jitter,  as  well  as  non-ideal  Rx  trackers  with  residual  tilt  jitters 
of  crAOA  =  0.25oy,  =  0.75crr,  =  aT.  For  the  worst  cases  of  a ^  =  ax  and  a AOA  =  Ot , 
the  Tx  beam  was  pointed  on-axis  to  the  Rx  and  the  beams  were  allowed  to  wander 
due  to  turbulence  without  a  tracking  system.  The  one-axis  beam  wander  variance  for 
a  collimated  beam  is  [69] 

°X  =  |)f7k(L2/  Cl(Z)dz~2Ljo  zCn(Z)d  Z  +  jo  (201) 

in  meters  squared.  The  Rx  AOA  variance  Tf2  in  radians  squared  is  given  by  Eq.  (47) 
in  Section  2  [69]. 

5.3.3  Modeling  fiber  coupling.  At  the  Rx,  the  light  is  coupled  into  a  single- 
mode  fiber  to  be  amplified  by  an  EDFA.  There  are  two  primary  considerations  when 
coupling  into  a  fiber;  the  angle  at  which  the  light  enters  the  fiber  core  and  the 
size  of  the  source  at  the  fiber  end  face.  The  coupling  efficiency  due  to  ray  angle 
is  unity  for  light  that  enters  a  fiber  at  less  than  the  fiber’s  numerical  aperture  angle 
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Ona  —  arcsin(A54).  For  single-mode  fibers,  the  efficiency  is  also  limited  by  an  ap¬ 
proximately  Gaussian  field  mode  guided  through  the  fiber.  This  held  can  be  described 
by  [11] 


Elp01 
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£oJo  (  —  ) 
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for  r  <  a 
for  r  >  a, 


(202) 

(203) 


where  u  =  1.55,  w  =  1.1428V  —  0.996,  and  V  =  2.405  is  the  normalized  frequency 
parameter  of  a  single  mode  fiber.  This  parameter  is  known  as  the  V  number  and  is 
given  by 

V  =  2naNA/X.  (204) 

In  Eqs.  (202)  and  (203),  the  Bessel  function  of  the  first  kind  is  denoted  by  J0,  and  K0 
is  the  modified  Bessel  function  of  the  second  kind.  Using  Eq.  (203)  determines  the 
irradiance  profile  in  the  fiber  and  enables  the  program  to  sum  up  the  optical  power 
coupled  in  the  fiber.  To  do  so,  the  held  at  the  image  plane  is  written  as  a  linear 
combination  of  basis  functions 

N 

E(x,  y)  =  Y^  a*Mx,  v),  (205) 

i= 1 


including  the  held  that  couples  into  the  hber 


fi(x,y)  =  ELP01(x,y) 


(206) 


and  the  held  that  fails  to  couple  into  the  hber 


f2(x,y)  =  1  -  ELP01(x,y). 


(207) 


Specifically,  the  coupled  held  can  be  determined  by  evaluating 


ECoupled  Q>lfl  {.Xj  y) , 


(208) 
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where 


IS  E(x,y)h(x,y) dz  dy 
./'/  f\  (x,  y  )  dx  d y 


(209) 


Entering  the  fiber  at  the  numerical  aperture  angle  9ma  yields  the  optimal  coupling 
efficiency  since  it  provides  the  smallest  spot  and  enters  the  fiber  end  face  within  9na- 
Therefore,  the  optimal  focal  length  /  of  a  lens  of  diameter  D  to  couple  a  uniform 
planar  field  into  a  fiber  is 


/  = 


D/2 


(210) 


tan  (6^) ' 

As  the  AOA  variance  increases,  the  optimal  focal  length  shifts  to  slightly  longer  focal 
lengths.  A  trade-off  could  be  made,  since  increasing  the  focal  length  also  increases 
the  spot  size.  Although  in  this  research,  Eq.  (210)  determines  the  focal  length  of  the 
Rx  telescope. 


5-4  Temporal  results 

Computer  simulations  of  airborne  single-Tx  and  double-Tx  FSOC  systems  were 
performed  for  the  scenarios  described  in  Table  4.  The  separation  distance  for  all  five 
scenarios  for  the  double-Tx  system  was  2 pc  =  31  cm.  Section  5.4.1  describes  the 
resulting  temporal  fade  statistics  of  the  detected  signal.  Section  5.4.2  plots  the  BER 
for  all  five  scenarios  and  for  different  techniques  used  to  improve  their  performance. 

5-4-1  Fade  statistics.  First,  a  fade  definition  is  required  to  determine  the 
fade  statistics.  In  Fig.  29,  the  plot  shows  the  optimal  fixed  threshold,  measured  signal 
(which  includes  the  variation  due  to  the  turbulence),  and  the  noise  <j\  associated  with 
the  measurement  of  a  ‘1’.  Notice  how  the  measured  signal  near  the  threshold  crosses 
the  threshold  often  due  to  the  measurement  noise  ay.  Instead  of  counting  each  one  of 
those  crossings  as  a  fade,  a  fade  occurs  only  if  the  variation  due  to  the  turbulence  in 
the  channel  causes  it  to  cross  the  threshold.  In  this  plot,  according  to  the  definition 
above,  there  are  two  distinct  fades  centered  at  about  3.9  and  4.05  ms  as  the  dashed 
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Figure  29:  Scenario  5  (HHL).  Measured  current  (solid  dark  line),  received  current 
(dashed  line),  optimal  fixed  threshold  (solid  light  line),  and  received  cur¬ 
rent  reduced  by  the  standard  deviation  of  the  measurement  noise  on  =  04 
(dotted  line). 

line  is  drops  below  the  wide  solid  line  It-  It  is  also  clear  to  see  in  this  plot  how  as  the 
received  signal  increases,  the  signal- dependent  measurement  noise  <j\  increases. 

Scenarios  1  and  3  did  not  experience  any  fades  for  the  double-Tx  cases  according 
to  the  definition  given  above  so  the  fades  are  not  plotted.  However,  these  scenarios 
did  experience  bit  errors  and  their  performance  are  shown  in  Section  5.4.2.  For  the 
other  three  scenarios,  Figs.  30,  31,  and  32  show  the  double-Tx  cases,  denoted  with 
dashed  line  have  shorter  and  less  frequent  fades  than  the  single-Tx  cases  denoted  by 
solid  line.  This  is  due  to  the  fact  that  the  fade  depths  are  reduced  when  the  multiple 
Txs  “smooth  out”  the  variation  in  the  received  power.  As  for  the  tracking  systems, 
the  better  the  tracker  the  shorter  the  fades,  but  the  greater  the  number  of  fades  per 
second.  The  tracking  systems  have  essentially  the  same  effect  as  just  turning  up  the 
power.  Each  of  the  scenarios  included  an  ideal  centroid  tracker  (o3  =  0)  denoted 
by  a  plain  line,  a  tracker  with  a  jitter  of  Oj  =  err/ 4  denoted  by  a  diamond  symbol, 
and  a  tracker  with  a  jitter  of  o3  =  3err/4  denoted  by  an  asterisk.  Figures  30  and 
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Figure  30:  Scenario  2  (HLL).  The  mean  fade  length  (a)  and  number  of  fades  per 
second  (b)  are  plotted  vs.  SNR.  The  solid  lines  are  for  a  single-Tx  and 
the  dashed  lines  are  for  a  double-Tx.  A  fade  is  defined  here  as  when  the 
raw  received  signal  s  drops  below  the  optimal  fixed  threshold  %t ■  The  data 
rate  was  Rb  =  B  =  8  GHz. 


32  show  results  for  the  most  severe  turbulence  cases  and  consequently  they  have  the 
longest  and  highest  number  of  fades.  The  only  difference  between  these  two  cases  is 
the  Greenwood  frequency  of  the  turbulence.  For  the  low  frequency  case  in  Fig.  30, 
the  fade  length  is  about  4.4  times  as  long,  but  there  are  4.0  times  fewer  fades  than  in 
the  high-/(5  case  in  Fig.  32. 


5-4-2  Bit  error  rate.  First,  the  discretized  PDFs  of  the  raw  received  signal 
p(is),  i.e.  the  variations  caused  by  atmospheric  turbulence,  for  each  of  the  scenarios 
mentioned  in  Table  4  for  the  single  transmitter  cases  are  compared  with  their  double 
transmitter  cases  for  all  five  scenarios.  If  the  PDF  is  heavily  weighted  to  the  left,  the 
chances  of  a  missed  detection  are  greater,  as  it  might  not  reach  above  the  threshold. 
Figure  33  shows  that  the  PDFs  of  the  received  signal  for  all  of  the  scenarios  shifted 
to  the  right  when  two  transmitters  were  used.  Even  for  the  \ow-R,sp  cases,  the  PDFs 
markedly  shifted  to  the  right  to  improve  the  performance.  This  shift  to  the  right 
reduces  the  probability  of  error  since  it  reduces  the  probability  of  a  missed  detection. 
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Figure  31:  Scenario  4  (HHH).  The  mean  fade  length  (a)  and  number  of  fades  per 
second  (b)  are  plotted  vs.  SNR.  The  solid  lines  are  for  a  single-Tx  and 
the  dashed  lines  are  for  a  double-Tx.  A  fade  is  defined  here  as  when  the 
raw  received  signal  s  drops  below  the  optimal  fixed  threshold  ix-  The  data 
rate  was  Rb  =  B  =  8  GHz. 
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Figure  32:  Scenario  5  (HHL).  The  mean  fade  length  (a)  and  number  of  fades  per 
second  (b)  are  plotted  vs.  SNR.  The  solid  lines  are  for  a  single-Tx  and 
the  dashed  lines  are  for  a  double-Tx.  A  fade  is  defined  here  as  when  the 
raw  received  signal  s  drops  below  the  optimal  fixed  threshold  %t ■  The  data 
rate  was  Rb  =  B  =  8  GHz. 
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Figure  33:  This  plot  shows  the  discretized  PDFs  of  the  received  signals  due  to  tur¬ 
bulence  p(is)  for  each  of  the  scenarios  for  the  single-Tx  (solid  line)  and 
double-Tx  (dashed  line)  cases.  Subplots  (a)-(d)  were  calculated  for  sce¬ 
narios  1,  2  &  5,  3,  and  4,  respectively. 


It  also  shifts  the  optimal  threshold  to  the  right,  reducing  the  probability  of  a  missed 
detection.  These  performance  improvements  are  quantified  with  the  BER  calculations 
in  the  next  section. 

Next,  the  BERs  of  ideal  and  realistic  adaptive  thresholds  are  compared  to  the 
optimal  fixed  case.  Recall  this  optimal  fixed  case  takes  into  account  the  PDF  of  the 
received  signal  p(is)  over  the  ensemble  of  the  runs  calculated  (see  Eq.  (188)).  There 
were  10  independent  realizations  with  1000  time  slices  for  each  realization.  The  time 
slices  were  determined  by  rs  =  1/(64/g)  for  each  of  the  scenarios.  Therefore,  each 
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independent  realization  covered  a  time  frame  of  over  15  Greenwood  time  constants, 
resulting  in  well  over  150  relatively  independent  realizations  per  scenario. 

It  is  clear  from  Fig.  34  that  the  BER  significantly  decreases  when  two  transmit¬ 
ters  are  used  (3-10  dB  depending  upon  the  scenario).  The  ideal  adaptive  threshold 
systems  improved  performance  by  up  to  5  dB  for  the  high-Rytov  cases  in  plots  (c) 
and  (d).  This  is  substantial  since  this  was  compared  to  the  optimal  fixed  threshold 
case  for  the  particular  scenario  which  used  the  actual  PDF  of  the  received  signal  to 
determine  the  optimal  threshold.  In  most  cases,  the  fixed  threshold  is  not  chosen 
in  such  a  precise  manner.  The  double-Tx  systems  outperformed  all  other  techniques 
even  though  improvements  due  to  the  adaptive-threshold  technique  were  up  to  5  dB. 
As  expected,  the  system  with  a  Rx  and  Tx  tracker,  an  ideal  adaptive  threshold,  and 
two  transmitters  performed  the  best. 

The  realistic  estimators  simulated  in  this  study  proved  to  improve  the  perfor¬ 
mance  in  all  cases.  The  performance  of  three  different  adaptive  threshold  systems 
are  compared  in  Fig.  34;  an  ideal  adaptive  threshold,  an  adaptive  threshold  with  an 
estimator  operating  at  fs  =  64 fa,  and  another  system  with  an  estimator  operating 
at  fs  =  16 fc-  For  a  single  transmitter  the  performance  for  the  fs  =  16 fc  estimator 
was  the  poorest  for  the  highest- Rytov  case  in  scenario  4  HHH.  For  this  and  all  other 
cases,  this  lower  sampling  rate  estimator  performance  greatly  improved  when  two 
transmitters  were  implemented.  The  single  transmitter  cases  have  more  variability 
in  the  received  irradiance  and  require  a  higher  fidelity  estimator  to  keep  up  with  the 
turbulence.  This  trickle-down  effect  indicates  multiple  transmitters  can  enable  the 
use  of  cheaper  lower-sampling-rate  estimators. 

Finally,  the  turbulence  effects  for  each  of  the  scenarios  are  compared  with  the 
BER  rate  performance  to  determine  causality.  The  only  difference  between  scenarios 
2  (HLL)  and  5  (HHL)  is  the  speed  of  the  turbulence  and  therefore,  as  expected,  their 
BERs  are  identical.  The  three  scenarios  with  high  Rytov  numbers  (scenarios  2,  4, 
and  5)  have  the  worst  performance,  but  the  improvement  provided  by  multiple  trans- 
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mitters  is  much  greater  for  these  scenarios.  The  combination  of  all  three  approaches 
provided  at  least  a  10  dB  gain  over  implementing  none  of  them  for  all  three  high  Ry- 
tov  number  scenarios.  Scenario  5  with  a  Rytov  number  lZsp  =  0.20  enjoyed  a  13  dB 
overall  improvement.  Due  to  the  improvements  afforded  by  these  multiple  techniques, 
the  high  Rytov  cases  of  7 Zsp  ~  0.20  were  on  par  with  Rytov  numbers  of  lZsp  ~  0.060 
without  these  techniques.  All  five  scenarios  tested  0.060  <  7 Zsph  <  0.20  were  within 
3  dB  of  each  other  when  all  of  the  improvement  techniques  were  implemented.  In  all 
scenarios,  the  coherence  diameter  r0  was  greater  than  the  diameter  of  the  receiver  D, 
therefore  changes  in  the  Rytov  number  had  a  much  larger  effect  than  r$.  If  D/tq  >  1 
the  phase  effects  due  to  the  turbulence  would  likely  have  had  a  larger  effect  on  the 
BER. 

Adaptive  thresholding  systems  provide  significant  improvement  over  optimal 
fixed  thresholds  for  both  single-Tx  and  double-Tx  systems,  providing  an  additional 
3-4  dB  over  both  systems.  As  long  as  the  estimator  kept  up  with  the  turbulence  the 
realistic  estimators  performed  well.  As  the  scintillation  effects  were  stronger  for  the 
single-Tx  high  Rytov  scenario  (HHH),  the  lower  bandwidth  estimator  performance 
lagged  behind  the  high  resolution  estimator  and  the  ideal  adaptive  threshold  system. 
The  primary  driver  to  the  adaptive  threshold  system  seemed  to  be  the  bandwidth, 
not  the  noisy  measurements.  This  is  most  likely  due  to  the  fact  that  the  SNR  in 
the  estimator  was  much  higher  than  for  the  digital  receiver  leg  (see  Fig.  28,  since  the 
required  bandwidth  was  so  much  lower  (kHz  vs.  GHz). 
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Figure  34:  BER  for  systems  with  an  optimal  fixed  threshold,  ideal  adaptive  threshold, 
fs  =  64 fG  estimator  adaptive  threshold,  and  fs  =  16  fc  estimator  adaptive 
threshold  for  single  Tx  (solid  lines)  and  double  Tx  (dashed  lines)  for  the 
tracker  case  of  Gj  =  crr/4.  Subplots  (a)-(d)  were  calculated  for  scenarios 
1,  2  &  5,  3,  and  4,  respectively.  The  data  rate  was  Rb  =  B  =  1  GHz. 
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VI.  Conclusion 


This  chapter  discusses  the  final  conclusions  of  this  research  including  the  chal¬ 
lenges  that  were  overcome  to  complete  the  project,  the  key  results,  first-time 
derivations,  and  recommendations  for  future  work. 

6.1  Challenges  met 

The  first  challenge  met  as  with  any  research  endeavor  dealt  with  defining  the 
problem  and  scoping  the  effort  to  a  reasonable  size  to  accomplish  the  research  in  the 
afforded  time.  There  has  been  quite  a  bit  of  research  into  free-space  optical  commu¬ 
nication  in  recent  years.  Fortunately,  FSOC  still  offers  many  different  fertile  areas  for 
research,  especially  for  the  air-to-air  link.  The  idea  for  using  multiple  transmitters 
came  from  the  realization  that  extended  objects  in  the  far  field  (like  the  moon  and  the 
planets)  do  not  twinkle  like  the  stars.  The  fluctuations  for  these  extended  objects  still 
occur,  but  they  are  averaged  over  the  extent  of  the  object.  Multiple  laser  transmitters 
incoherent  with  each  other  and  adequately  separated  act  in  much  the  same  way. 

Secondly,  for  all  the  literature  on  multiple  Tx-Rx  systems,  there  was  little  in¬ 
formation  about  how  far  apart  to  have  the  Txs  to  accomplish  good  averaging.  For 
the  multiple-Tx  case,  the  next  challenge  was  to  determine  analytic  equations  for  the 
requisite  anisoplanatic  separations  for  the  phase  and  amplitude  perturbations.  This 
analytic  research  built  on  developments  from  Roggemann  and  Welsh,  [66]  Sasiela,  [69] 
Andrews  and  Phillips,  [2]  and  Fried  [26]  to  determine  the  anisoplanatic  effects  in  tur¬ 
bulence.  Many  tacks  were  taken  to  determine  these  relations,  starting  with  a  study  of 
the  differential  tilt  variance  relations  in  Sasiela  to  determine  the  tilt  anisoplanatic  an¬ 
gle.  Determining  the  separation  required  for  the  tilt  anisoplanatic  angle  proved  to  be 
intractable.  A  new  tack  was  used  to  determine  tilt  anisoplanatic  angle  by  way  of  the 
phase  anisoplanatic  angle.  This  approach  started  with  the  phase  isoplanatic  effects 
and  determined  for  the  first  time  a  relation  for  the  phase  anisoplanatic  angle.  As  for 
the  amplitude  effects,  using  the  principle  of  reciprocity,  the  requisite  Tx  separation 
was  determined  by  using  the  correlation  width  pcw  previously  used  to  predict  how  Rx 
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size  effects  scintillation.  Similarly,  this  pcw  relation  can  be  determined  by  referring  to 
the  log-amplitude  structure  function.  The  problem  with  the  amplitude  effects  is  that 
they  saturate  for  strong  turbulence,  making  these  relations  only  valid  for  the  weak 
turbulence  regime  and  increasing  the  importance  of  simulation  and  experimentation. 

The  third  challenge  was  to  determine  the  practical  feasibility  of  a  multiple  trans¬ 
mitter  approach  for  the  airborne  scenario.  Since  these  analytic  relations  had  not  been 
studied  for  the  airborne  scenario,  they  needed  to  be  tested  to  determine  which  geome¬ 
tries  might  be  feasible.  Scenarios  for  the  air-to-air,  ground-to-air,  ground-to-ground, 
and  air-to-ground  were  tested  in  random  independent  realizations  of  turbulence.  The 
Rx  model  included  a  tracking  system  for  both  the  Rx  and  the  Tx  as  well  as  modeling 
the  coupling  into  a  single  mode  fiber. 

The  fourth  challenge  was  to  design  a  method  to  simulate  these  multiple  scenar¬ 
ios.  A  wave-optics  simulation  for  each  of  the  scenarios  was  developed  from  scratch, 
with  each  requiring  special  considerations  for  sampling.  These  considerations  were 
dealt  with  by  implementing  a  split-step  Fresnel  propagation  system  to  reduce  the 
sampling  requirements.  The  long  propagations  using  the  split-step  Fresnel  method 
still  require  large  arrays  to  satisfy  the  sampling  constraints  in  Eqs.  (163)  to  (166). 
In  this  research  2048  x  2048  arrays  were  used  to  calculate  the  propagations  through 
10  independent  phase  screens  for  the  propagations  in  Chapter  V.  Special  coding 
considerations  were  used  to  properly  handle  these  large  arrays  to  avoid  “out  of  mem¬ 
ory”  conditions.  The  analysis  that  followed  determined  that  an  air-to-air  multiple  Tx 
system  was  feasible  since  the  optimal  separation  was  only  about  31cm  apart  for  the 
100  km  long,  10  km  altitude  path. 

Next,  implementing  the  time-series  wave-optics  simulation  presented  new  chal¬ 
lenges  in  modeling.  A  new  relation  for  determining  the  temporal  sampling  of  a  wave- 
optics  simulation  was  developed  that  takes  into  account  the  temporal  aspects  of  the 
amplitude  perturbations.  The  Fourier-series-based  phase  screens  enabled  the  calcula¬ 
tion  of  time-correlated  phase  screens  without  requiring  the  screens  to  be  repeated. 
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Finally,  the  fixed  optimal  threshold  required  knowledge  of  turbulence  PDF,  but 
determining  the  PDF  of  these  multiple  turbulence  scenarios  with  multiple  Txs  was 
intractable.  Therefore,  numerical  PDFs  were  determined  from  the  simulations  and 
used  to  determine  the  fixed  optimal  threshold.  The  ideal  adaptive  threshold  system 
used  in  other  research  [13,18,72]  was  compared  to  a  new  linear  estimator  designed  for 
this  research.  For  all  calculations  the  measurement  noise  sources  included,  thermal, 
shot,  and  ASE  noise  from  the  EDFA. 

6.2  Key  results 

6.2.1  Anisoplanatic  effects.  This  research  designed  a  multiple- transmitter 
system  to  average  different  effects  of  turbulence  by  first  determining  the  requisite  sep¬ 
aration  distances.  Roggemann,  Fried,  Saseila  and  others  determined  isoplanatic  con¬ 
ditions  so  that  the  effects  of  two  paths  were  relatively  correlated  [26,51,66,69,70,84], 
This  separation  allows  AO  systems  to  measure  and  correct  for  these  effects  in  real 
time.  Here,  the  research  of  Roggemann,  Welsh,  and  Fried  was  extended  to  determine 
when  these  effects  become  uncorrelated  [26,66].  The  phase  anisoplanatic  angle  9^ind 
was  first  derived  in  support  of  this  research  [46,47,49]  and  defines  the  angle  over  which 
the  phase  effects  between  the  propagation  paths  of  two  point  sources  are  nearly  uncor¬ 
related.  This  relation  increases  with  outer  scale  and  turbulence  strength  and  defines 
the  separation  where  the  beams  wander  nearly  independently.  This  angle  describes 
the  minimum  separation  for  uncorrelated  beam  wandering  which  determines  the  opti¬ 
mal  separation  distance  for  untracked  multiple  Tx  systems.  The  likely  application  for 
an  untracked  system  is  a  ground-to-ground  last-mile  communication  system,  since  the 
system  is  stationary  and  the  outer  scale  is  small,  allowing  for  reasonable  separation 
distances  (i.e.  25-100  cm).  As  long  as  there  are  sufficient  number  of  Txs  with  wide 
enough  beams,  a  system  could  be  designed  to  maintain  at  least  one  beam  on  the  tar¬ 
get  without  the  need  for  tracking.  This  phase  anisoplanatic  angle  also  corresponds  to 
the  tilt  anisoplanatic  angle,  since  the  lowest-order  phase  effects  manifest  themselves 
as  tilt  and  beam  wander. 
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Since  this  research  is  concerned  with  air-to-air  laser  communication  a  separation 
distance  is  required  for  a  tracked  system.  In  this  case,  a  separation  angle  is  required 
to  average  the  higher-order  phase  and/or  scintillation  effects.  The  beams  must  be 
separated  far  enough  apart  that  the  scintillation  effects  are  uncorrelated.  This  occurs 
when  the  variance  of  the  difference  between  the  two  effects  is  equal  to  twice  the 
variance  of  a  single  path  effect.  In  previous  research  [2],  this  condition  was  determined 
as  the  correlation  width  pcw  for  amplitude  variations  of  two  paths.  This  previous  work 
applied  to  how  large  a  receiver  aperture  needs  to  be  averaging  out  scintillation  effects. 
Here,  the  principle  of  reciprocity  determines  how  far  apart  the  transmitters  need  to 
be  so  the  scintillation  effects  can  be  averaged  at  the  receiver.  At  this  separation, 
the  hot  spots  and  signal  drop  outs  due  to  strong  scintillation  across  the  aperture 
are  “averaged  out.”  Research  has  been  done  to  determine  this  requisite  separation 
distance  [4,5,59,60],  but  it  has  not  been  accomplished  for  air-to-air  scenarios  and  not 
put  in  terms  of  the  correlation  width  pcw.  The  optimal  Tx  separation  distance  for 
reducing  BER  was  1-3  times  pc  for  the  air-to-air  and  ground-to-air  scenarios. 

The  improvement  due  to  implementing  two  transmitters  can  be  scaled  to  a  small 
degree  for  multiple  transmitters.  A  limit  to  this  improvement  does  exist,  since  sepa¬ 
rating  beams  much  greater  than  2 pc  results  in  diminishing  improvement  as  shown  in 
Chapter  IV  [46,49].  The  mechanism  for  declining  performance  improvement  depends 
upon  how  the  beams  are  separated.  If  the  beams  are  separated  in  angle,  the  average 
AOA  of  the  beams  at  the  receiver  causes  the  beams  to  land  off-center  of  the  detector. 
If  the  beams  are  too  far  apart  angularly,  the  beams  on  average  will  land  off  of  the 
detector  in  the  focal  plane,  consequently  reducing  the  power  at  the  detector.  If  the 
beams  are  parallel,  separating  them  further  causes  each  beam  to  land  further  off  axis 
at  the  aperture,  reducing  the  individual  beam’s  received  power  and  increasing  the 
scintillation. 

Since  scintillation  effects  begin  to  decorrelate  for  fairly  small  separations,  the 
transmitters  only  need  to  be  separated  by  approximately  2 pc  for  air-to-air  scenarios 
of  100  km  in  length  at  4km  in  altitude  and  above. 
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Conventional  AO  systems  only  correct  for  phase  and  cannot  correct  for  strong 
scintillation,  but  here  the  different  transmitters  “average  out”  the  strong  scintillation 
effects  by  incoherently  summing  up  the  two  beams  in  the  receiver.  Additionally,  since 
the  Rx  aperture  diameter  was  smaller  than  r0,  an  AO  system  would  not  have  provided 
much  improvement. 

6.2.2  Temporal  aspects  and  fade  statistics.  Research  by  Haas,  Anguita, 
Lee,  Peleg,  and  others  into  multiple  transmitters  have  not  considered  the  temporal 
impacts  of  multiple  transmitters  [4,5,29,30,43,59].  This  research  determined  that  the 
multiple  transmitter  approach  not  only  reduces  the  length  of  fades,  but  also  reduces 
the  number  of  fades  by  simultaneously  reducing  the  variation  and  increasing  the 
received  power.  The  averaging  effect  of  two  transmitters  reduced  the  BER  and  the 
length  of  faded.  This  effect  had  not  been  studied  before.  If  further  error  reduction  is 
required,  an  inter leaver/FEC  receiver  could  be  implemented.  Multiple  Tx  techniques 
and  adaptive  thresholding  do  not  require  physical  elements  to  adaptively  compensate 
for  turbulence,  but  they  reduce  fade  lengths.  This  temporal  aspect  of  multiple-Tx 
systems  had  not  been  studied  before.  In  addition,  the  length  of  fades  for  higher  air 
velocities  are  shorter  than  the  lower  air  velocities.  Shorter  fade  lengths  require  shorter 
interleavers,  reducing  data  latency  and  making  shorter  interleavers  more  effective. 

6.2.3  Optimal  fixed  and  adaptive  thresholding.  Burris  used  estimates  of  sig¬ 
nal  mean  and  variance  to  determine  the  optimal  fixed  threshold  using  the  LRT  [13]. 
Crabtree  calculated  the  optimal  fixed  threshold  using  turbulence  PDF  estimates, 
which  is  better  since  it  takes  into  account  an  estimate  of  the  turbulence  PDF’s 
shape  [18].  This  research  and  Schmidt’s  research  used  the  actual  discretized  PDFs 
from  the  simulations  [72],  The  ideal  adaptive  threshold  system  in  other  research 
[13,18]  was  compared  to  a  new  linear  estimator  designed  for  this  research.  All  previ¬ 
ous  cases  of  adaptive  thresholding  systems  were  implemented  for  the  single-Tx  case. 
For  high-Rytov  scenarios,  the  low-sampling-rate  estimator  performance  greatly  im¬ 
proved  when  two  transmitters  were  implemented.  The  single  transmitter  cases  have 
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more  variability  in  the  received  irradiance  and  require  a  higher  fidelity  estimator  to 
keep  up  with  the  turbulence.  This  trickle-down  effect  indicates  multiple-Tx  systems 
enable  the  use  of  cheaper  lower-sampling-rate  estimators. 

These  and  other  trickle-down  effects,  like  the  reduction  in  fade  lengths  for  mul¬ 
tiple  transmitter  systems  indicate  the  techniques  implemented  in  this  research  were 
complimentary.  In  addition,  the  performance  improvements  for  these  techniques  were 
additive  and  pushed  the  performance  of  high-Rytov  (JZsp  ~  0.20)  turbulence  cases  to 
the  performance  of  much  more  benign  turbulence  regimes  ( 7Zsp  ~  0.06). 

This  research  enables  free-space  optical  communication  systems,  promising  to 
provide  covert,  difficult  to  jam  or  intercept,  high-speed,  broadband  connectivity  to 
airborne  platforms  [86] .  This  connectivity  will  provide  unparalleled  situational  aware¬ 
ness  to  the  warfighter.  The  next  steps  in  this  research  include  designing  hardware  to 
implement  this  system  in  a  hardware  in  the  loop  test  followed  by  flight  testing. 

6. 3  Recommendations  for  future  work 

Throughout  this  research  effort  multiple  areas  for  future  research  were  identified. 
These  areas  are  summarized  below. 

6.3.1  Diversity  techniques.  To  illustrate  the  optimal  separation  only  two 
beams  were  used,  so  the  full  impact  of  this  approach  has  yet  to  be  explored.  Others 
have  shown  that  four  beams,  for  instance,  can  be  very  effective  in  a  ground-to-ground 
link  scenario  [5].  Multiple-transmitter  configurations  could  be  studied  to  determine 
the  point  of  diminishing  return  and  what  configuration  is  optimal  for  both  parallel 
and  angular  separations.  Since  the  improvements  for  two  Txs  peaked  for  separations 
of  about  2-3  pc  and  diminishing  returns  occurred  for  separations  larger  than  about 
3  pCl  4-5  transmitters  is  potentially  the  limit  before  improvements  diminish,  but  this 
needs  to  be  demonstrated  conclusively.  In  addition,  a  comparison  of  parallel  versus 
angular  separated  beams  might  provide  some  interesting  design  considerations. 
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A  combination  of  multiple  transmitters  and  conventional  AO  might  provide 
significant  improvement  without  the  need  for  sophisticated  multi-conjugate  adaptive 
optics  systems.  The  multiple  transmitters  would  provide  the  improvement  in  terms 
of  reducing  scintillation  while  the  AO  system  would  handle  the  phase  perturbations, 
likely  improving  the  fiber  coupling.  A  comparison  of  conventional  AO  and  multiple 
transmitter  system  performance  for  varying  turbulence  conditions  could  also  provide 
some  design  considerations  for  future  air-to-air  FSOC  systems. 

6.3.2  Optical  receiver  design.  The  receiver  coupled  the  laser  light  into  a 
single  mode  EDFA.  The  optimal  focal  length  used  in  this  research  was  based  on  a 
plane  wave  incident  on  the  aperture.  If  the  AOA  variance  is  large  at  the  Rx,  longer 
focal  lengths  would  be  required  to  maintain  an  angle  less  than  the  NA  angle  to  couple 
more  light  into  the  fiber.  A  trade-off  would  have  to  be  made  as  well  since  increasing  the 
focal  length  also  increases  the  spot  size.  AOA  variance  due  to  turbulence  and  multiple- 
transmitter  systems  could  provide  a  rich  design  space  for  investigating  tradeoffs  in  the 
focal  length  of  the  receiving  telescope. 

Fiber  amplifiers  offer  advantages  in  simplicity,  efficiency,  and  high-speed  per¬ 
formance.  More  detailed  fiber  amplifier  modeling  could  distinguish  other  advantages 
including  automatic  gain  control  afforded  by  the  gain  saturation.  In  addition,  schemes 
to  improve  the  coupling  efficiency  could  ameliorate  their  performance. 

The  propagation  of  lasers  through  extended  turbulence  causes  significant  hot 
spots  and  dropouts.  Therefore,  an  optical  and/or  electrical  AGC  would  likely  provide 
significant  performance  improvement  with  inherently  dynamic-range-limited  detec¬ 
tors.  EDFAs  do  afford  some  optical  automatic  gain  control,  but  research  into  other 
techniques  could  provide  additional  improvements. 

6.3.3  Signal  processing  techniques.  Adapting  the  threshold  provided  signif¬ 
icant  improvement.  This  improvement  was  compared  to  the  optimal  fixed  threshold 
case,  which  is  difficult  to  implement  since  it  requires  knowledge  of  the  turbulence 
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PDF.  If  comparisons  were  made  with  more  realistic  thresholds,  the  performance  im¬ 
provement  would  have  been  more  significant.  Comparisons  with  multiple  methods 
including  a  threshold  determined  solely  from  the  average  power  level  over  a  time  pe¬ 
riod  might  provide  a  more  realistic  improvement  metric.  In  addition,  the  estimator 
used  in  this  work  could  be  extended  to  an  optimal  estimation  of  the  decision  threshold 
from  a  noisy  measurement. 

6.3.4  Tracking  systems.  A  point-ahead  angle  study  for  very  fast  moving 
transmitters  and  receivers  might  lead  to  a  determination  that  higher-order  phase 
compensation  using  conventional  adaptive  optics  systems  might  be  difficult  at  best, 
unfeasible  at  worst. 

The  centroid  trackers  used  in  this  research  performed  relatively  poorly.  This  was 
due  to  the  significant  scintillation  of  the  extended  turbulence  propagation.  Techniques 
to  track  the  multiple  hot  spots  or  select  the  optimal  tracking  path  could  improve  the 
performance.  (Many  times  the  centroid  was  located  in  a  low  intensity  location  on 
the  aperture)  In  addition,  new  tracking  and  pointing  algorithms  could  be  developed 
for  multiple-transmitter  systems,  since  the  separation  required  to  average  out  the 
scintillation  is  smaller  than  the  tilt  isoplanatic  angle.  Scintillation  does  adversely 
affect  centroid  determination  in  a  centroid-tracking  system  [48]. 

In  this  work,  each  transmitter  had  its  own  tracking  system.  Since  the  trans¬ 
mitter  separation  is  likely  to  be  within  the  tilt  isoplanatic  angle,  a  single  tracker  for 
this  configuration  might  be  adequate.  Multiple  transmitters  alone  might  be  able  to 
improve  the  tracking  system  since  it  would  reduce  the  hot  spots  in  the  receiver.  This 
reduction  in  hot  spots  could  enable  a  better  centroid  tracking  system  since  it  would 
not  jump  around  due  to  the  scintillation  effects. 

6.3.5  Large  communication  system  design.  On  a  system-level  basis  one  could 
study  the  design  consideration  of  a  high-altitude  optical  communication  hub.  This  hub 
would  provide  satellite  uplink/downlink  to  coordinate  and  route  information  from  the 
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future  Transformational  Satellite  Communications  System  (TSAT).  An  airborne  relay 
mirror  system  could  be  considered  to  provide  over-the-horizon  high-speed  line-of-sight 
communication  for  the  battle  space.  It  would  provide  much  faster  communication 
than  satellite  communication  systems  and  would  provide  a  high  bandwidth  approach 
to  long  haul  ground-to-ground  communication  systems. 
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